1  No.  FHWARD  75-52 


TE 

662 
.A3 
no. 

«jj-*"  RATION  OF  FLOOD  RISK  FACTORS  IN  THE 
^nrtsiGN  OF  HIGHWAY  STREAM  CROSSINGS 


Vol.  II  Analysis  of  Bridge  Backwater  Experiments 


Dept  of  Ti i 


il  1977 
Final  Report 


Document  is  available  to  the  public  through 
the  National  Technical  Information  Service, 
Springfield,  Virginia  22161 


Prepared  for 

FEDERAL  HIGHWAY  ADMINISTRATION 
Offices  of  Research  &  Development 
Washington,  D.  C.  20590 


Foreword 

This  report  describes  large  scale  laboratory  tests  of  bridge  backwater 
studies  and  will  be  of  interest  to  bridge  and  hydraulics  engineers 
interested  in  the  effects  of  bridges  on  water  surface  profiles. 

This  report  is  Volume  II  of  a  five  volume  series  presenting  the  results 
of  the  Water  Resources  Engineers,  Inc.,  research  on  "Evaluation  of 
Flood  Risk  Factors  in  the  Design  of  Highway  Stream  Crossings."  This 
is  the  last  of  the  five  volumes  to  be  distributed;  the  other  four 
volumes  were  distributed  as  they  were  completed  over  a  2-year  period 
beginning  in  1975.  The  research  was  conducted  for  the  Federal  Highway 
Administration,  Office  of  Research,  Washington,  D.C.,  under  Contract 
DOT-FH-1 1-7669. 

This  document  is  disseminated  under  the  sponsorship  of  the  Department 
of  Transportation  in  the  interest  of  information  exchange.  The  United 
States  Government  assumes  no  liability  for  its  contents  or  use  thereof. 

Sufficient  copies  of  the  report  are  being  distributed  by  FHWA  Bulletin 
to  provide  a  minimum  of  two  copies  to  each  FHWA  Regional  office,  one 
copy  to  each  FHWA  Division  office,  and  two  copies  to  each  State  highway 
agency.  Direct  distribution  is  being  made  to  the  Division  offices. 


•£or   Charles  F.  Scheffey 

Director,  Office  of  Research 
Federal  Highway  Administration 


TECHNICAL  REPORT  STANDARD  TITLE  PAGE 


1.     Report  No. 


FHWA-RD-75-52 


2.     Government  Accession  No. 


3.     Recipient's  Catalog  No. 


EVALUATION  OF  FLOOD  RISK  FACTORS  IN 
HIGHWAY  STREAM  CROSSINGS.  Vol.  II. 
Bridge  Backwater  Experiments 


THE  DESIGN  OF 
Analysis  of 


5.     Report  Dote 

April    1977 


6.     Performing  Organization   Code 


7.     Author's) 

A.  J.  Knepp 
R.  P.  Subinski 


M.  T.  Tseng 


-  -  "■  , 


B.     Performing  Organization   Report  No. 


WRE  20810 


9.     Performing  Organization  Name  and  Address 

Water  Resources  Engineers,  In 
8001  Forbes  Place 
Springfield,  Virginia  22151 


10.     Work   Unit  No. 


1 1 .     Contract  or  Grant  No. 

D0T-FH-1 1-7669 


12.     Sponsoring  Agency  Name  and  Address 

Offices  of  Research  and  Development 
Federal  Highway  Administration 
U.  S.  Department  of  Transportation  ,/ 
Washington,  D.C.  20590 


13.     Type  of  Report  ond  Period  Covered 


Final   Report 


£ 


14.     Sponsoring  Agency  Code 


15.     Supplementary  Notes 

FHWA  Contract  Administrator: 
FHWA  Contract  Manager: 
FHWA  Project  Manager: 


Gilbert  Trainer 
J.  Sterling  Jones 
Roy  E.  Trent 


FHWA  Implementation  Manager: 
Dan  O'Connor 


.  Abstrocr  . 

lis  report  presents  the  results  of  a  series  of  experiments  carried  out  using  a  large 
:ale  physical  model  of  backwater  caused  by  highway  stream  crossings  of  wide,  heavily- 


16.     Abstroct 

Th 

scale  phy_ 

vegetated  flood  plains.  Variables  studied  included  wingwall  and  spillthrough  abute- 
ments,  spur  dikes,  prototype  bridge  openings  between  100  and  1100  ft  and  prototype 
flows  between  8,000  and  36,000  cfs.  An  important  aspect  of  these  experiments  was  the 
use  of  large  scale  roughness  elements  protruding  the  free  surface  to  provide  the  nec- 
essary ranges  of  roughness  to  simulate  prototype  conditions.  Five  combinations  of 
roughness  densities  in  patterns  were  tested  in  the  course  of  experimentation.  Primary 
emphasis  was  on  crossings  normal  to  flow;  however,  eccentric  and  skewed  crossings  were 
also  tested. 

Results  of  the  study  included  an  analysis  of  the  sensitivity  of  water  surface 
elevation  to  the  principal  variables,  quantification  of  both  the  location  and  height 
of  the  maximum  backwater,  and  data  for  possible  future  study  and  modification  of 
existing  backwater  methods. 

The  other  volumes  of  this  study  are: 

I  Experimental  Determination  of  Channel  Resistance  for  Large  Scale  Roughness  75-51 

III  Finite  Element  Model  for  Bridge  Backwater  Computation  FHWA-RD-75-53 

IV  Economic  Risk  Analysis  for  Design  of  Bridge  Waterways  FHWA-RD-75-54 

V  Data  Report  for  Spur  Dike  Experiments  FHWA-RD-75-55 


17.   Key  Words 

Backwater  computation,  analysis  of 
variance,  sensitivity  analyses,  highway 
stream  crossing  design 


18.     Distribution   Statement 


No  Restrictions.  This  document  is 
available  to  the  public  through  the 
National  Technical  Information  Service, 
Springfield,  Virginia  22151 


19.     Security  Clossif.  (of  this  report) 

Unclassified 


20.     Security  Clossif.  (of  this  page) 

Unclassified 


21.   No.  of  Pages 

182 


22.     Price 


Form  DOT  F  1700.7  (8-69) 


TABLE  OF  CONTENTS 

Page  No, 


NOTICE 

ii 

TABLE  OF  CONTENTS 

i  i  i 

ACKNOWLEDGEMENTS 

V 

LIST  OF  FIGURES 

vi 

LIST  OF  TABLES 

X 

I. 

INTRODUCTION 

1 

STUDY  OBJECTIVES 

1 

BACKGROUND 

2 

II. 

LITERATURE  REVIEW 

7 

III. 

HYDRAULIC  MODEL  DESIGN 

22 

GENERAL  DESCRIPTION 

22 

DYNAMIC  SIMILITUDE 

23 

GOVERNING  FLOW  EQUATIONS  IN  BRIDGE  WATERWAYS 

26 

REQUIREMENTS  OF  MODEL  SCALES 

30 

CHANNEL  RESISTANCE 

36 

IV. 

EXPERIMENTAL  EQUIPMENT 

42 

FLUME  DESCRIPTION 

.42 

ROUGHNESS  ELEMENTS 

44 

DESCRIPTION  OF  CROSSING  CONFIGURATIONS 

51 

V. 

EXPERIMENTAL  PROCEDURE 

58 

GENERAL  DESCRIPTION 

58 

HYDRAULIC  FLUME  TESTS 

63 

DATA  COLLECTION 

64 

QUALITY  CONTROL 

80 

OPERATIONAL  PROBLEMS 

80 

VI. 

PRESENTATION  AND  ANALYSIS  OF  DATA 

83 

11 


VII. 
VIII. 


TABLE  OF  CONTENTS 
(Cont'd) 


INTRODUCTION 

DIMENSIONAL  ANALYSIS 

DEFINITION  OF  VARIABLES 

ANALYSIS 

TWO-DIMENSIONAL  ANALYSIS  OF  WATER  SURFACE 
ELEVATIONS 

ANALYSIS  OF  VARIANCE 

ANALYSIS  OF  DATA  FOR  ROUGHNESS  PATTERN  C 

ANALYSIS  OF  TOTAL  BACKWATER  COEFFICIENT,  K* 

ANALYSIS  OF  MAXIMUM  BACKWATER  LOCATION  U*/t>) 

ANALYSIS  OF  MAXIMUM  BACKWATER  HEIGHT  (h*/y  ) 

CONCLUSIONS 

RECOMMENDATIONS 
DATA  COLLECTION 
DATA  ANALYSIS  AND  INTERPRETATION 

LIST  OF  REFERENCES 

APPENDIX 


Page  No 


83 

84 

84 

85 

89 

90 

95 

103 

133 

153 

164 

166 

166 

168 

170 

A-l 

m 


ACKNOWLEDGEMENTS 


Numerous  personnel  from  both  the  Federal  Highway  Administration  and 
Water  Resources  Engineers,  Inc.,  have  made  valuable  contributions  to  the 
completion  of  this  project.  While  it  is  impossible  to  individually  acknowledge 
the  contributions  of  each  individual,  certain  do  deserve  special  recognition. 
Among  these  are  the  staff  of  the  Environmental  Control  Group,  Office  of 
Research  and  Development,  Federal  Highway  Administration,  particularly 
J.  Sterling  Jones  and  Roy  E.  Trent. 

Within  WRE,  overall  project  direction  has  been  the  responsibility, 
initially,  of  G.  Kenneth  Young  and  in  the  last  three  years  of  the  study  effort 
Robert  P.  Shubinski.  Ming  Te  Tseng  designed  the  experiment  and  directed  the 
data  collection  activities.  The  data  analysis  described  in  this  volume  was 
performed  and  documented  by  Anthony  J.  Knepp;  Mr.  Knepp  and  Dr.  Shubinski 
assume  all  responsibility  for  the  conclusions  and  recommendations  of  this 
volume. 

Of  the  many  others  at  WRE  who  worked  on  the  project  a  particular  note 
of  thanks  is  extended  to  Stanley  W.  Zison,  Richard  A.  Schmalz,  Larry  S.  Costello, 
John  L.  Matticks  and  Bert  Black.  For  the  final  format  of  the  document  a 
grateful  acknowledgement  is  extended  to  Phyllis  Wei ner,  who  drafted  Chapter  II, 
"Literature  Survey,"  and  provided  technical  editing;  Luanne  Wohler,  who  was 
responsible  for  the  technical  drafting;  and  to  Kathy  Settle  and  Joan  Lagasse, 
who  typed  the  manuscript. 


iv 


LIST  OF  FIGURES 
Figure  No.  Page  No, 

II-l      Definition  Sketch  of  Backwater  Reach  9 

1 1-2      Backwater  Ratio  for  Channels  Having  Vertical  - 

Faced  Constrictions  with  Square-Edged  Abutments  11 

I 1-3      Energy  Losses  for  Schneider  Method  13 

I 1-4      Backwater  Analysis  17 

I II-l      Types  of  Roughness  Elements  38 

IV-1      Basic  Flume  Characteristics  43 

IV-2      Definition  of  Roughness  Pattern  46 

IV-3      Manning's  "n"  for  Roughness  Patterns  as  a  Function 

of  Flow  48 

IV-4      Wingwall  Abutment  Design  53 

IV-5      Spill  through  Embankment  55 

IV-6      Spur  Dike  Design  56 

V-l      Crossing  Configurations  58 

V-2      Flow  Lines  for  Typical  Normal  Crossing  59 

V-3      Abutment  Shapes  61 

V-4      Velocity  Measurement  Locations  for  Centered 

Crossing  Configuration  74 

V-5      Skewed  Model  Velocity  Sampling  Configuration  76 

V-6      Eccentric  Model  Velocity  Sampling  Configuration  76 

V-7      Point  Gages  77 

V-8      Beam  Apparatus  Used  for  Locating  Maximum 

Backwater  82 

VI-1      Analysis  of  Variables  Affecting  Water  Surface 

Elevation  97 

VI-2      Analysis  of  Variables  Affecting  Water  Surface 

Elevation-Excluding  Spur  Dikes  101 

VI-3      Maximum  Backwater  Coefficient  (K*  vs  M)         109 


Figure  No. 

VI 

-4 

VI- 

-5 

VI- 

-6 

VI- 

-7 

VI- 

-8 

VI- 

-9 

VI- 

-10 

VI- 

-11 

VI- 

-12 

VI- 

-13 

VI- 

-14 

VI- 

-15 

VI- 

-16 

VI- 

-17 

VI- 

-18 

VI-19 


LIST  OF  FIGURES 
(Cont'd) 


Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spill  through  Roughness  Pattern  A 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spill  through  Roughness  Pattern  B 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spill  through  Roughness  Pattern  C 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spill  through  Roughness  Pattern  D 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spur  Dike  Roughness  Pattern  A 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spur  Dike  Roughness  Pattern  B 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spur  Dike  Roughness  Pattern  C 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Spur  Dike  Roughness  Pattern  D 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Wingwall  Roughness  Pattern  A 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Wingwall  Roughness  Pattern  B 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Wingwall  Roughness  Pattern  C 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Wingwall  Roughness  Pattern  D 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Wingwall  Roughness  Pattern  P 

Experimental  Eccentricity  Definition 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Skew  y  =  15° 

Maximum  Backwater  Coefficient  (K*  vs  M 
for  Skew  y  =  15°  (centered) 


Page  No. 

110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 
125 

129 

130 


VI 


LIST  OF  FIGURES 
(Cont'd) 


Figure  No.  Page  No 

VI-20     Maximum  Backwater  Coefficient  (K*  vs  M) 

for  Skew  y   45°  131 

VI-21     Maximum  Backwater  Location  [z*/b  vs  M) 

for  Wingwall  Roughness  Pattern  P  137 

VI-22     Maximum  Backwater  Location  [a*/b  vs  M) 

for  Wingwall  Roughness  Patterns  A  and  B  138 
VI-23     Maximum  Backwater  Location  U*/b  vs  M) 

for  Wingwall  Roughness  Patterns  C  and  D  139 
VI-24     Maximum  Backwater  Location  {i*/b   vs  M) 

for  Spill  through  Roughness  Patterns  A  and  B  140 
VI-25     Maximum  Backwater  Location  [l*/b   vs  M) 

for  Spillthrough  Roughness  Patterns  C  and  D  141 
VI-26     Maximum  Backwater  Location  (£*/b  vs  M) 

for  Spur  Dike  Roughness  Patterns  A  and  B  142 
VI-27     Maximum  Backwater  Location  (£*/b  vs  M) 

for  Spur  Dike  Roughness  Patterns  C  and  D  143 
VI-28     Plot  of  Residuals  vs  £*/b  148 

VI-29     Maximum  Backwater  Height  (h*/y  vs  M) 

for  Spillthrough  Roughness  Patterns  A  and  B  156 
VI-30     Maximum  Backwater  (h*/y0  vs  M) 

for  Spillthrough  Roughness  Patterns  C  and  D  157 
VI-31     Maximum  Backwater  Height  (h*/y0  vs  M) 

for  Wingwall  Roughness  Patterns  A  and  B  158 
VI-32     Maximum  Backwater  Height  (h*/y0  vs  M) 

for  Wingwall  Roughness  Patterns  C  and  D  159 
VI-33     Maximum  Backwater  Height  (h*/y0  vs  M) 

for  Wingwall  Roughness  Pattern  P  160 


vn 


LIST  OF  FIGURES 
(Cont'd) 

Figure  No.  Page  No, 

VI-34     Maximum  Backwater  Height  (h*/y  vs  M) 

for  Spur  Dike  Roughness  Patterns  A  and  B      161 

VI -35     Maximum  Backwater  Height  (h*/y0  vs  M) 

for  Spur  Dike  Roughness  Patterns  C  and  D      162 


vm 


Table  1 

to. 

III- 

-1 

IV- 

-1 

V- 

-1 

V- 

-2 

LIST  OF  TABLES 


Page  No 

Values  of  a,  3,  cq  and  Bi  40 

Roughness  Combinations  47 

Summary  Table  of  Experimental  Runs  63 

Experimental  Investigation  Chart  Centered 

Opening  Runs  65 

V-3     Experimental  Investigation  Chart  Eccentric 

Runs  68 

V-4     Experimental  Investigation  Chart  Skew  Runs       69 

V-5     Experimental  Investigation  Uniform  Flows         70 

V-6     Spillthrough  Runs  Not  Analyzed  (All  Roughness 

Pattern  B)  71 

V-7     Raw  Data  Collected  for  Each  Run  72 

V-8     Numbers  of  Moveable  Gages  78 

V-9     Stationary  Gage  Locations  78 

VI-1     Analysis  of  Variance  93 

VI-2     ANOVA  of  K*  With  and  Without  Spur  Dikes  106 

VI-3     Testing  of  K*  With  and  Without  Spur  Dikes  at 

95  Percent  Confidence  Level  107 

VI-4     ANOVA  of  K*  for  15°  and  45°  Skew  Crossings  127 

VI-5     Testing  of  K*  for  Skew  Crossing  at  95  Percent 

Confidence  Level  128 

VI-6     ANOVA  of  £/b  With  and  Without  Spur  Dikes  134 

VI-7     Testing  of  £*/b  With  and  Without  Spur  Dikes  135 

VI-8     Correlation  Coefficients  for  Independent 

Variables  145 

VI-9     Log-Linear  Regressions  for  £*/b  147 

VI-10    Beta  Weights  150 


IX 


LIST  OF  TABLES 
(Cont'd) 


Table  No.  Page  No, 

VI-11    ANOVA  of  h*/y0  for  Roughness  Pattern  C  With 

and  Without  Spur  Dike  Data  154 

VI-12    Testing  of  h*/y_  for  Roughness  Pattern  C  at 

95  Percent  Confidence  Level  155 


CO 

EC 

o 
►— 
u 
< 

z 
o 

CO 

cc 


E        -o 


.5        2-" 


<       S 


E    »    5-'  ° 

f  SIS 


I  I  g 

i 


*  s  I 


g  a 


6    E 


u.     O   -L    O  I 
I  I 


CJ 

zz 

U       103 

61 

8! 

i! 

91 

SI 

H 

n     \l\ 

II 

01 

6 

8 

L 

9 

s 

* 

c 

z 

1 

ll!l!!lll 

lllllllll 

III!  llll  llll  !!!i 

lilillill 

I!!!!lll! 

mini! 

lllllllll 

lllllllll 

lllllllll  lllllllll 

lllllllll 

llll 

llll 

llll 

llll 

llll 

illl 

lllllllll 

llll 

11!! 

lllllllll 

llll 

llll 

llll 

llll 

!!lllu;l 

1 
llll'llil 

o 
o 


9  8  1 


I      !  I   '  I  I  1   '      I  I  I   I  I  I  *      I  |  I  I  I  II      I  I  I  II  I  I      I  |  I  III  I      I  [  I    I   II  I 


I  |  <  I  III      III  llll      INI   III 


1  inches 


Is  I 


5   E    p  'J2 

O     O     ft)     v 


J    C    cr    it 


e 


O     O     <M    © 


O    O   O    O    (*)    o    o 


£  ».  E 


—     13      4> 


cr    c    c-    o 


C      3      O 

3      O     .C 

o    a    in 


c    9    E 


■g    <*  ill    ^  if    - 


=    E 


XI 


I.  INTRODUCTION 


This  report  is  the  second  in  a  series  of  five  volumes  comprising  the 
final  report  for  the  study  entitled  Evaluation  of  Flood  Risk  Factors  in  the 
Design  of  Highway  Stream  Crossings ,  authorized  by  the  Federal  Highway  Adminis- 
tration (FHWA)  under  Contract  No.  DOT-FH-1 1-7669.  The  overall  objective  of 
the  study  is  to  develop  an  engineering  systems  analysis  method  to  enhance  the 
decision-making  process  in  the  design  of  highway  stream  crossings.  This 
method  applies  economic  risk  techniques  as  well  as  standard  hydraulic  and 
hydrologic  factors  in  the  design  of  bridge  waterways. 

The  present  volume  describes  the  analysis  of  the  relationship  between 
combinations  of  bridge  openings,  abutment  shapes,  flows  and  roughness  densities 
and  the  resultant  backwater.  The  analysis  is  based  on  data  taken  during  a 
two-year  period  of  operation  of  a  large  hydraulic  flume.  The  relationships 
which  were  defined  in  the  analysis  can  be  extended  into  the  computerized 
Finite  Element  Model  simulation  of  the  flood  flow  and  backwater  as  described  in 
Volume  III. 


STUDY  OBJECTIVES 

The  principal  purposes  of  this  phase  of  the  study  were: 


1.  To  collect  raw  data  for  the  analysis  and  quantification  of  the 
effects  of  abutment  shapes,  channel  roughness,  and  crossing 
skew  and  eccentricity  on  the  magnitude  and  location  of 
maximum  backwater  under  varied  flow  conditions. 

2.  To  provide  additional  information  needed  to  update  and  possibly 
modify  the  FHWA  backwater  analysis  method. 


3.  To  evaluate  two-dimensional  aspects  of  backwater  and  the  effects 
of  main  channel  and  overbank  flows. 


BACKGROUND 

Highway  design  engineers  have  long  understood  the  problems  associated 
with  constricting  flow  in  natural  waterways.  The  introduction  of  embankments 
and  piers  onto  the  flood  plain  and  in  the  main  channel  alters  the  hydraulic 
behavior  of  the  waterway  and  can  appreciably  affect  the  flood  stage  of  the 
river  during  a  major  flooding  event  miles  upstream.  The  result  of  a  higher 
stage  may  be  property  and  personal  damages  over  what  would  have  occurred  had  the 
waterway  not  been  obstructed. 

Awareness  of  the  problem  has  led  state  highway  departments  to  include 
an  estimation  of  the  hydraulic  efficiency  of  the  highway  crossing  structure 
along  with  the  design  specifications  and  cost  of  the  structure.  Usually  the 
50-year,  100-year,  or  maximum  flood  of  record  is  used  as  the  design  flood  for 
the  hydraulic  analysis.  The  importance  and/or  cost  of  the  structure  predicates 
the  exact  choice  of  the  design  storm.  Once  this  choice  is  made,  the  stage 
and  velocity  reached  in  the  waterway  are  determined  for  both  the  unconstricted 
flow  and  the  flow  with  the  crossing  in  place.  This  additional  stage  height 
attributable  to  the  constriction  of  the  flow  is  known  as  backwater. 

Numerous  methods  exist  for  the  estimation  of  backwater;  some  of  them 
are  described  in  detail  in  the  literature  review.  Currently,  three  distinct 
methods  are  extensively  used:  the  Federal  Highway  Administration  Method 
(FHWA),  the  Geological  Survey  Method  (USGS),  and  the  Corps  of  Engineers 
Method  (HEC-2).  The  FHWA  method  [1]*  is  derived  from  the  work  of  Liu, 
Bradley  and  Plate  [2]  at  Colorado  State  University  in  1957.  The 

*Bracketed  numbers  refer  to  references  placed  at  the  end  of  the  text 
material . 


expression  derived  from  the  work  for  the  computation  of  backwater  is: 
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where 
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K* 

An2 

Vn2 
A- 
Ai 


=  total  backwater  (ft), 

=  total  backwater  coefficient, 

=  kinetic  energy  coefficient, 

=  gross  water  area  in  constriction  measured  below  normal 
stage  (sq  ft) , 

=  average  velocity  in  constriction  at  normal  depth  or 
Q/An2  (fps), 

=  water  area  at  section  4  where  normal  stage  is 
reestablished  (sq  ft),  and 

=  total  water  area  at  section  1,  (where  the  backwater  is 

greatest),  including  that  produced  by  the  backwater  (sq  ft). 


The  total  backwater  coefficient  (K*)  is  empirically  determined  and  accounts 
for  the  effect  of  piers,  angle  of  skew,  abutment  shape  and  eccentricity.  These 
effects  are  assumed  to  be  additive  and  are  functionally  related  to  the  dimen- 
sionless  parameter  M,  where  M  is  the  ratio  of  the  flow  that  can  pass 
unimpeded  through  the  bridge  opening  to  the  total  flow.  This  method  is 
restricted  since  it  provides  only  a  one-dimensional  analysis  along  the  centerline 
and  it  does  not  allow  for  natural  river  meandering,  variable  channel  conveyance 
capacities,  and  cross  flow  from  the  channel  to  the  overbank  and  vice  versa. 


The  USGS  method  experimentally  derived  in  1965  by  Kindsvater  and 
Carter  [3]  is  based  on  energy  and  continuity  relationships  developed  between  the 
upstream  reaches  and  the  location  of  the  vena  contracta.  Energy  loss  is 
calculated  by  use  of  a  conveyance  method  and  a  division  of  upstream  areas  into 
portions  where  energy  is  being  dissipated  due  to  both  friction  and  contraction. 


A  recent  report  written  by  Schneider  et  al .  [4]  presents  a  modified 
method  which  updates  the  old  USGS  method  and  reportedly  improves  its  accuracy 
to  a  considerable  extent.  Schneider  employs  a  more  accurate  estimate  of 
streamline  length  between  sections  using  potential  flow  theory.  He  also 
proposes  a  slightly  modified  estimation  of  energy  losses  in  the  flow  expansion 
reach.  However,  this  method  also  is  based  on  a  single-dimension  formulation 
of  energy  relationships  and  thus  is  limited  as  to  the  type  of  problem  which 
can  be  solved  and  the  accuracy  of  the  solution. 

The  third  widely  applied  backwater  method  is  that  developed  by  the 
Corps  of  Engineers  in  1973  and  called  HEC-2  [5].  This  computerized  steady 
state  model  computes  the  water  surface  profile  for  river  channels  of  any  cross 
section  for  either  subcritical  or  supercritical  flow  conditions.  The  effects 
of  various  hydraulic  structures  such  as  bridges,  culverts,  weirs,  embankments 
and  dams  are  considered  in  the  computation.  The  computational  procedure 
applies  Bernoulli's  Theorem  for  the  total  energy  at  each  cross  section  and 
Manning's  formula  for  the  friction  head  loss  between  cross  sections.  Total 
head  loss  is  calculated  using  an  iterative  process  and  includes  transition 
loss,  losses  through  bridges,  and  through  culvert  structures. 

A  fourth  method,  used  as  part  of  the  present  study,  uses  the  finite 
element  techniques  for  the  solution  of  the  equations  describing  flow  and  bulk 
transport  of  incompressible  fluids.  A  computer  model  solving  these  equations, 
i.e.,  continuity  and  momentum  relationships,  was  used  in  simulating  backwater 
profiles  in  a  two-dimensional  plane  at  a  bridge  site  located  in  Tallahalla, 
Mississippi,  with  good  results.  Among  the  advantages  offered  by  this 
approach  are:  (1)  a  two-dimensional  analysis;  (2)  irregular  cross  sections 
and  channel  bends  can  be  specified;  (3)  energy  losses  are  divided  into 
frictional  and  eddy  losses  allowing  the  effects  to  be  estimated  separately; 
(4)  roughness  can  be  specified  in  the  detail  to  which  it  is  known;  (5)  boun- 
dary effects  can  be  determined;  and  (6)  multiple  runs  are  possible  with  only 
slight  additional  effort.  The  principal  disadvantages  are:  (1)  the  technique 


presently  does  not  handle  transitional  flow  across  the  constriction;  (2)  the 
model  must  be  calibrated  with  the  field  data  to  ascertain  correct  values  of 
the  eddy  coefficients;  (3)  the  model  reliability  is  sensitive  to  user's  ability 
to  establish  a  grid  for  the  finite  elements;  and  (4)  the  model  requires  sub- 
stantial computer  capability. 

This  report  presents  the  analysis  of  information  collected  from  the 
operation  over  a  two-year  period  of  a  22.7  x  184  ft  flume  used  in  studying  the 
effects  on  backwater  of  various  bridge  openings,  abutment  shapes,  flows 
and  roughness  densities.  The  study  was  restricted  to  subcritical  flow  conditions 
The  range  of  bridge  openings  varied  from  three  to  eleven  feet  while  abutment 
shapes  included  45°  wingwall,  spillthrough  and  elliptical  spur  dikes.  Though 
most  model  openings  were  centered  in  the  flume  and  perpendicular  to  the  flow, 
both  eccentric  and  skew  crossings  were  also  investigated. 

The  interactive  effects  of  flood  plain  and  main  channel  flow  areas 
on  backwater  were  simulated  through  the  use  of  artificial  large  scale  roughness 
elements  placed  in  the  flume  to  represent  the  differing  effect  on  flow  of  both 
the  heavily  vegetated  overbank  and  the  main  channel  areas.  The  actual  degree 
of  roughness  produced  in  the  flume  was  even  higher  than  observed  because  the 
distortion  of  the  scale  model  requires  additional  roughness  to  satisfy  the 
law  of  dynamic  similitude  between  the  model  and  the  prototype.  Conceptually, 
the  energy  loss  due  to  the  viscosity  of  the  fluid  and  the  turbulence  generated 
by  the  roughness  elements  was  treated  as  energy  loss  attributable  to  frictional 
forces  alone. 

The  choice  of  the  density  and  location  of  the  roughness  elements  is 
the  subject  of  Volume  I  of  this  study  entitled  Experimental  Determination  of 
Channel  Resistance  for  Large  Scale  Roughness s    Report  No.  FHWA-RD-75-51  , 
available  from  the  National  Technical  Information  Service. 


Chapter  II  of  this  report  is  a  literature  review  of  current  one- 
dimensional  backwater  methods.  Chapter  III  presents  an  analysis  of  the 
hydraulic  model  design  which  forms  the  theoretical  basis  for  the  experimental 
configurations.  The  experimental  equipment  is  described  in  Chapter  IV 
and  the  important  procedural  methods,  both  in  the  collection  and  the  reduc- 
tion of  raw  data,  are  described  in  Chapter  V.  Operational  and  analysis 
problems  encountered  either  in  the  collection  or  the  analysis  of  data 
are  included  in  this  section.  Chapter  VI  describes  the  dimensional  analysis, 
which  to  a  large  extent  directs  the  data  analysis  presented  in  the  same 
chapter.  Conclusions  and  recommendations  follow  in  Chapters  VII  and  VIII, 
respectively. 

Listings  of  the  processed  and  sorted  data,  including  condensed  output 
summaries,  is  available  in  a  separately  bound  appendix  available  through  the 
Federal  Highway  Administration  Office  of  Research  (HRS-42),  Washington,  D.C. 
20590.  Condensed  output  summaries  for  each  flume  run  are  contained  in  an 
appendix  attached  to  this  report. 


II.  LITERATURE  REVIEW 


In  the  last  half  of  the  1800' s  and  during  the  early  part  of  this 
century  the  effects  of  stream  constriction  and  pier  design  on  backwater  were 
first  defined  by  D'Aubuisson,  Eytelwein,  Nagler,  Hutton,  Rehbock,  Lane,  Yarnell 
and  many  others.  The  earlier  investigators  used  mathematical  aralysis  while 
Nagler  and  his  contemporaries  developed  hydraulic  models  and  gathered  data 
from  real-world  observations.  Nagler  focused  primarily  on  the  backwater  effects 
of  piers,  comparing  various  shapes,  sizes  and  orientations  to  flow  [6]. 
Rehbock,  in  Germany,  also  used  hydraulic  models  to  determine  the  backwater 
effects  of  piers,  using  three  classifications  of  flow.  He  found  that  the 
maximum  backwater  caused  by  pier  obstruction  depends  on  the  channel -bridge 
waterway  contraction  ratio,  the  Froude  number  of  the  unobstructed  flow  and  the 
pier  geometry  [7] . 

In  the  1930's  D.L.  Yarnell  [8]  experimentally  verified  existing 
backwater  formulas  for  pier  effects  and  in  the  process  came  to  a  number  of 
important  conclusions: 

1.  The  height  of  the  backwater  due  to  bridge  piers  varies 
directly  as  the  depth  of  unobstructed  channel, 

2.  For  the  lower  velocities,  a  lense-shaped  nose  and  tail 
offer  minimum  resistance, 

3.  The  optimum  ratio  of  pier  length  to  width  probably  varies 
with  the  velocity  and  is  generally  between  4  and  7, 

4.  Placing  the  piers  at  an  angle  with  the  current  has  an 
insignificant  effect  on  the  amount  of  backwater  if  the 
angle  is  less  than  10°, 

5.  Placing  the  piers  at  an  angle  of  20°  or  more  with  the 
current  materially  increases  the  amount  of  backwater,  the 
increase  depending  upon  the  quantity  of  flow,  the  depth, 
and  the  channel  contractions. 


E.W.  Lane  [9]  was  the  first  to  study  the  effects  of  channel  constriction 
on  open-channel  flow,  laying  the  groundwork  for  backwater  formulas  by  developing 
an  empirical  discharge  coefficient  from  the  discharges  and  difference  of 
surface  elevation  upstream  and  downstream  from  the  constriction. 

A  working  but  generally  inapplicable  method  for  computing  the  maximum 
backwater  height  due  to  channel  constriction  was  developed  in  the  early  1950's 
by  H.J.  Tracy  and  R.W.  Carter  [10].  This  method  was  based  on  the  laboratory 
work  of  Kindsvater,  Tracy  and  Carter  [3]  to  estimate  discharge  through  a 
channel  constriction  (see  Figure  II-l).  The  relationship  established,  based  on 
one-dimensional  energy  and  continuity  relationship,  is  given  below: 


where 


Q   =  C  bh3/2g  [Ah  +  oi  V?/2g  -  Ef(l_3)l  (1) 

Q   =  discharge  through  the  opening  (cfs), 

C   =  Kindsvater  discharge  coefficient, 

b   =  width  of  the  opening  (ft), 

h3  =  water  surface  elevation  at  Section  3  (ft) 

(Sections  1,  2,  3  and  4  are  shown  in  Figure  II-l), 

Ah  =  difference  in  surface  elevations  at  Section  1  and  3  (ft), 

al     =  energy  correction  factor  to  account  for  unequal 
velocity  distribution  at  Section  1, 

Vf/2g  =  velocity  head  (ft),  and 

F_f(i-3)  =  head  loss  due  to  friction  between  Sections  1  and  3. 

The  relationship  is  a  convenient  one  for  computing  discharge  from  field  measure- 
ments at  an  existing  bridge.  It  is  not,  however,  in  terms  of  the  maximum 
backwater  h*,  which  is  of  interest  to  highway  designers. 

Parameters  shown  in  the  figure  but  not  included  in  Equation  1  are: 


B  =  flood  plain  width, 

„  =  contraction  coeffii 
e 

ht  =  maximum  backwater  above  normal  elevation  at  Section  1 


C   =  contraction  coefficient,  and 
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(c)   Elevation    adapted  to  assumption  of  zero  friction  loss 


FIGURE  II-l.  DEFINITION  SKETCH  OF  BACKWATER  REACH 
(Source:  Reference  10) 


Five  cross  sections  are  shown  in  Figure  Il-lb.  Section  0  refers  to 

the  area  upstream  of  the  constriction  unaffected  by  backwater  effects 

emanating  from  the  constriction,  section  1  defines  the  cross  section  at  the 

location  of  the  maximum  backwater,  section  2  is  located  in  the  center  of  the 

constriction,  section  3  is  defined  by  the  location  of  the  vena  contracta 

which  is  the  point  at  which  the  live  stream  boundaries  are  closest,  and  finally, 

section  4  is  the  cross  section  where  normal  flow  depth  returns.  Through 

dimensional  analysis  it  was  shown  that  the  Kindsvater  discharge  coefficient,  C, 

is  a  function  of  approach  channel  shape  and  roughness,  Froude  number 

(F  = ) ,  the  geometry  of  the  constriction  and  the  contraction  ratio  (m). 

bh3/gh3 


The  contraction  ratio,  m,  is  defined  as 

m  =  1  -  y~ 

KB 

K   =  conveyance 

K,   =  conveyance  of  section  of  length  b 

KD  =  conveyance  of  section  of  length  B. 

D 

Conveyance  is  defined  from  Manning's  equation  as: 
K  =L«AR2/3 


where 


where 


n  =  Manning's  roughness  factor 
A  =  cross-sectional  area 
R  =  hydraulic  radius. 
Note  that  if  channel  slopes  are  constant 


m  =  1  -  — 

% 
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where 


=  flow  carried  in  width  b,  and 


qp  =  flow  carried  in  width  B, 


% 
>B 


Tracy  and  Carter  [10]  constructed  a  family  of  curves  showing  the 
relationship  between  C  and  the  dimensionless  parameters  which  exhibited 
a  functional  response.  These  curves  are  not  included  here. 

Tracy  and  Carter  later  found  that  -r-  was  primarily  a  function  of 
the  percentage  of  contraction.  The  quantity  Ah  was  defined  as 


Ah  = 


4 


2gb2h32C2 


Vi2 

"  *1      2g"+   E 


f(l-3) 


(2) 


The  backwater  h*  was  then  found  by  multiplying  Ah  by  the  backwater  ratio 

h* 

—r-     ,  selected  from  measured  laboratory  data  summarized  in  Figure  II-2. 

These  calculations  assume  a  standard  condition,  i.e.,  vertical  constrictions 
with  square-edged  abutments  in  a  rectangular  (and  horizontal)  flume  for 
three  variations  of  channel  roughness. 
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FIGURE  II-2.  BACKWATER  RATIO  FOR  CHANNELS  HAVING  VERTICAL-FACED  CONSTRICTIONS 
WITH  SQUARE-EDGED  ABUTMENTS.  Source:   Reference  10. 
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The  application  of  this  method  is  limited  for  several  reasons: 

1.  Since  Ah  is  an  independent  variable,  the  estimation  of  h*  through 

J.  is  a  trial  and  error  process  not  practicable  for  the  direct 

Ah 

use  of  h*  in  bridge  design  calculations. 

2.  The  small  scale  of  the  hydraulic  model  limits  the  application  of 
relationships  between  the  contraction  ratio,  channel  roughness 
and  constriction  geometry. 

3.  The  hydraulic  model  used  in  the  experiments  was  designed  with 
a  level-bottom  channel,  in  which  for  a  given  discharge  the 
velocity,  depth  and  energy  gradient  of  unobstructed  flow  vary 
from  section  to  section;  that  is,  the  flow  is  nonuniform  and 
there  is  no  standard  flow  condition. 

4.  The  scope  of  the  methodology  was  limited  to  steady  tranquil  flows 
through  single-opening  constrictions  and  does  not  include  effects 
of  scour  around  bridge  piers  and  abutments  or  boundary  friction 
due  to  heavy  vegetal  growth  along  the  stream  banks. 

The  method  proposed  by  Tracy  and  Carter  [10]  has  most  recently  been 
updated  and  modified  by  field  data  collected  by  the  USGS  in  conjunction  with 
the  FHWA  and  the  State  Highway  Departments  of  Mississippi,  Alabama  and 
Louisiana.  In  their  report,  not  yet  published,  Schneider  et  al .  [4]  propose 

refinements  to  both  the  contraction  reach  and  expansion  reach  energy  loss 
formulations  as  given  in  Equation  1,  which  is  repeated  here  for  convenience 
of  referral : 


Q  =  C  bh3/2g[Ah  +  ai  Vf/2g  -  Ef(l_3)] 

The  method  of  solution  for  Equation  1  for  backwater  is  an  iterative  approach 
using  a  "step-backwater"  technique.  Normal  water  surface  elevation  is  first 
calculated  followed  by  constricted  flow  calculations.  The  distinction  between 
the  Carter  method  and  the  Schneider  method  is  that  modifications  for  head 
loss  computation  proposed  by  Schneider  are  incorporated  into  the  solution 
algorithm. 

The  energy  loss  term  due  to  friction,  Ef(-   \,  represents  the  head  loss 
occurring  in  moving  from  cross  section  1  to  cross  section  3  as  shown  in 
Figure  II-3.  The  "L"  values  shown  on  the  horizontal  axis  represent  distances 
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between  cross  sections.  Thus  L  ,  is  the  center  line  distance  from  cross  section 
"d"  to  cross  section  2.  Expanding  Ew   \  according  to  Figure  1 1-3 : 


:f(l-3)=^  +  ^^ 


KtK 


1-n     -d-3     K  2 


(3) 


where 


Lw  = 

Ld  = 
L  = 

Ki  *■ 
V 

K3  = 


distance  from  the  location  of  the  maximum  backwater  to 
the  toe  of  the  spur  dikes  (ft), 

length  of  the  spur  dike  (ft),  and 

length  of  the  bridge  abutment. 

approach  conveyance  at  section  1  (ft3/sec), 

that  portion  of  the  approach  conveyance,  Kls  corresponding 
to  the  bridge  width  b  (ft3/sec), 

conveyance  of  the  cross  section  at  the  toe  of  the  spur 
dikes  (ft3/sec),  and 

conveyance  of  the  cross  section  at  the  abutments  (ft3/sec). 
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FIGURE   1 1-3.      ENERGY  LOSSES  FOR  SCHNEIDER  METHOD 
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Each  of  the  three  components  of  Equation  3  represents  the  energy  loss 
due  to  friction  over  the  stated  length.  The  components  of  Equation  3  represent 
the  geometric  mean  of  the  energy  slopes  at  the  end  cross  sections  of  each 
reach  times  the  distance  between  the  sections.  Schneider  assumed  the  uniform 
flow  equation  for  some  cross  section  i  having  conveyance  K-  as: 


Q  =K.Sf 


Vz 


which  implies 


l/2 

Sf    =  Q/K. 
The  geometric  mean  friction  slope  can  then  be  represented  by 

5f(i-j)^ 


K.K. 


where  subscripts  i  and  j  refer  to  upstream  and  downstream  conveyances, 
respectively. 

Schneider  et  al  [4]  propose  two  modifications  to  Equation  3  which 

represent  changes  in  approach  and  expansion  reach  energy  losses.  First,  in 

the  modification  for  the  contraction  reach,  L,  Schneider  assumes 

w 

that  for  large  contraction  ratios  the  average  streamline  length  L  is  actually 

much  greater  than  the  straight  line  distance  L.  For  the  determination  of 

3  w 

the  average  streamline  length,  L  ,  he  studied  an  idealized  constriction  using 
potential  flow  theory.  Average  streamline  length  as  a  function  of  the  length 
of  the  approach  reach  divided  by  the  width  of  the  opening  and  the  contraction 
ratio  were  determined.  Tabular  values  are  not  presented  in  this  report,  but 
can  be  found  in  reference  4. 

Second,  Schneider  calculates  head  losses  in  the  flow  expansion  reach  as 
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where 


and 


\(3-k) 


'f  C3-f) 


hf  (3-4)    +  h 

bQ2 
KcK.n 

Q2 
2gA42 

(u, 

) 


ah)     +   2| 


3  A^ 


+  a: 


fe) 


(4) 
(5) 

(6) 


Equation  5  represents  the  frictional  head  loss  between  sections  3  and 

4.  Note  that  normal  depth  is  assumed  to  occur  one  bridge  length  downstream 

from  section  3.  K   is  the  conveyance  at  normal  depth  for  section  4  and  K  is 
hX\  c 

the  smaller  of  the  conveyances  K  or  K3 .  Equation  6  represents  flow  expansion 
losses  for  an  idealized  expansion  in  open  channel  flow  according  to  Henderson  [11] 
Here  beta  and  alpha  are  momentum  and  energy  coefficients  while  A  is  cross 
sectional  area.  Subscripts  refer  to  cross  section  locations. 

Analyzing  the  data  collected  for  31  floods,  Schneider  reports  that  the 
unmodified  Geological  Survey  and  FHWA  methods  underpredict  backwater  by  45  and 
47  percent,  respectively.  Using  the  proposed  modification  the  mean  error  for 
computing  backwater  was  six  percent. 

A  simple  formula  for  computing  backwater  was  proposed  by  C.F.  Izzard 

in  1954  [12]  in  his  discussion  of  the  Kindsvater  and  Carter  paper  [3].  This 

provided  the  basis  for  the  beginnings  of  the  method  which  became  known  gener- 
ally as  the  FHWA  method.  It  is 

V  2 
hi  -  hn  =  hf  =  K. 


n2 


vb  2g 


where  V 


_  Q 


,  ,  —  is  a  hypothetical  velocity  and  K,  is  an  empirical  backwater 
coefficient.  The  data  compiled  by  Tracy  and  Carter  were  used  by  Izzard  [13] 
to  develop  a  simple  correlation  of  the  maximum  backwater/contraction  ratio  and 
Froude  number.  The  correlation  was  useable  for  rough  estimates  but  did  not 
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consider  important  variables  such  as  channel  slope  and  roughness,  skew  and 
eccentricity  of  the  crossing  and  piers. 

Liu,  Bradley  and  Plate  [2]  at  the  Colorado  State  University  used 
Izzard's  concept  for  a  series  of  more  sophisticated  hydraulic  model  experiments. 
They  defined  the  backwater  expression  by  applying  the  principle  of  conservation 
of  energy  between  the  point  of  maximum  backwater  upstream  from  the  contraction 
and  the  point  downstream  at  which  the  flow  has  reached  normal  stage.  This 
method  assumes  that  the  channel  in  the  vicinity  of  the  bridge  is  essentially 
straight,  the  cross  sectional  area  of  the  stream  is  reasonably  uniform  and  the 
gradient  of  the  bottom  is  constant  between  cross  sections  1  and  4  (Figure  1 1-4) . 
The  analysis  applies  only  to  steady  tranquil  flow. 

The  theoretical  development  by  Liu  et  al  [2]  is.  essentially  an 
application  of  one-dimensional  energy  relationships  between  cross  sections  1  and 
4,  respectively.  The  analysis,  condensed  from  reference  14,  follows.  Equating 
the  two  sections: 

2 

so^-o+  yi   +^g~=^  +  ^"+  hT 


where 


S   =  channel  slope  (ft/ft), 

L(j-it)   =  distance  between  sections  1  and  4, 

y\     =   depth  of  water  at  section  1  , 

aiVi2 

=  average  velocity  head  at  section  1 


2g 

yk     =   normal  depth  of  water  at  section  4,  and 
hy  =  total  energy  loss  between  sections  1  and  4. 

Since  the  frictional  loss  term  in  feet/foot  of  channel  is  contained 
in  the  total  energy  loss  term  hT,  these  quantities  can  be  eliminated  from  both 
sides  of  the  equation  and  hT  replaced  by  h,.  Rearranging  thus  yields: 
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a4V        aiV22 


2g~+hb 


The  additional  loss  h.  can  be  expressed  as  a  product  of  velocity  head  in  the 
constriction  at  normal  flow  and  an  energy  loss  coefficient,  K*. 


hb  =K*Tg 


ni 


Replacing  h,  with  the  above  expression  and  y1   -   y^  with  h*  yields 


a9V  2 


h*  =  k* 

ni     K   2g 


n2 


ai+Vt^2   a1V1: 


2g     2g 


The  analysis  assumes  that  the  cross  sectional  areas  at  1  and  4  are 
essentially  the  same  and  from  continuity--  A^  =  Ai+V^  =  A  V   --allows: 

A  V 
w     n2  n2   „  ■ 
V4  =   A    ,  and 


An2Vnz 


Replacing  Vi  and  V4  with  the  above  and  factoring  V   leaves 


h* 


vn22 

=  K*«2  -2g~  +  ai 


A  2 
n2 


A  2 
n2 

A!2 


V   2 

n2 
2g 


(7) 


Through  the  aid  of  dimensional  analysis  a  backwater  coefficient  K*  was 
determined  experimentally  and  was  found  to  be  a  function  primarily  of 


1.  Stream  constriction  as  measured  by  the  opening  ratio  M,  where 


2.  Type  and  shape  of  bridge  abutment--wingwall ,  spj 11th rough,  etc., 

3.  Number,  size,  shape  and  orientation  of  piers  in  the  constriction, 

4.  Eccentricity  or  asymmetric  position  of  bridge  on  the  flood  plain, 

5.  Skew  (bridge  crosses  river  of  flood  plain  at  other  than  90  degrees), 
and 

6.  Froude  number. 

The  hydraulic  model  used  for  these  experiments  had  a  variable  flume 
slope,  two  flume  widths  and  two  types  of  bed  roughness.  The  experiments 
conducted  with  this  flume  used  unobstructed  flow  data  (normal  depth  and  Froude 
number)  taken  for  each  slope,  discharge  rate  and  roughness  type  in  order  to 
make  a  realistic  judgment  of  the  effects  of  the  constriction  on  the  flow 
through  the  constriction  and  the  resulting  backwater.  The  flume  constriction 
was  formed  with  models  of  three  types  of  bridge  abutments  and/or  piers.  The 
specific  purpose  of  the  study  was  to  determine  the  maximum  height  of  backwater 
for  a  given  constriction.  Crossing  conditions  studied  were: 

1.  Simple  normal  crossing, 

2.  Abnormal  stage-discharge  condition, 

3.  Dual -bridges  contraction, 

4.  Bridge  girders  partially  submerged, 

5.  Skew  crossing, 

6.  Eccentric  crossing,  and 

7.  Piers  with  and  without  abutments. 

U.S.  Geological  Survey  field  surveys  in  Mississippi  conducted  to 
verify  the  laboratory  tests  showed  that  the  model  application  is  limited  to 
medium  size  bridges  (i.e.,  structures  up  to  220  feet  in  length  with  flood 
plains  up  to  .5  mile  in  width)  since  it  underpredicts  backwater  for  bridges 
with  a  larger  channel/opening  ratio  [14]. 

The  calculations  derived  from  the  Colorado  State  University  study 
and  improvements  made  since  are  incorporated  in  the  1970  edition  of  the  Bureau 
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of  Public  Roads  Design  Series,  "Hydraulics  of  Bridge  Waterways"  [14].  This 
publication  presents  the  most  advanced  method  in  wide  use  for  estimating 
backwater  in  a  noncritical,  uniform  flow  situation.  The  manual  also  presents 
means  for  making  preliminary  estimates  of  the  effects  of  spur  dikes,  critical 
stage-discharge  conditions,  and  Manning  roughness  coefficients  for  stream 
channels  and  flood  plains. 

The  backwater  effects  of  flood  plain  flow  were  the  focus  of  an  article 
by  E.M.  Laursen  [15]  which  followed  the  Colorado  State  University  investigation, 
Laursen  noted  the  two-dimensional  aspects  of  flood  flow  in  a  wide  valley  where 
the  bridge  embankment  forces  overbank  flow  to  move  laterally  to  the  bridge 
opening  and,  on  the  downstream  side  of  the  constriction,  to  move  laterally 
back  to  the  flow  boundary.  Laursen  postulated  that  the  downstream  diffusion 
in  section  4  is  possibly  a  significant  factor  in  calculating  upstream  backwater 
and  based  his  analysis  for  sections  1  and  4  on  momentum  principles.  He 
pointed  out  that  previous  experiments  with  hydraulic  models  actually  simulated 
cross  sections  2  and  3  because  of  the  relatively  small  dimensions  involved, 
leaving  unconsidered  the  flow  behavior  of  cross  sections  1  and  4  such  as  would 
be  characteristic  of  the  prototype. 

C.L.  Yen  has  contributed  to  the  description  of  the  characteristics  of 
flood  plain  flow  versus  channel  flow  and  in  a  1971  paper  with  D.R.  Overton  [16] 
presented  a  method  for  subdividing  the  flood  plain  cross  section  in  order  to 
evaluate  turbulent  flow  resistance  for  each  subsection.  Yen  and  Overton 
suggested  that  since  Manning's  n-values  for  both  main  channel  and  flood  plain 
are  lower  than  those  for  wide  rectangular  channels  of  the  same  materials,  the 
n-values  reflect  not  only  the  roughness  but  also  the  effects  of  channel 
geometry.  While  these  ideas  have  not  been  approached  in  the  current  study, 
they  offer  an  avenue  for  greater  accuracy  in  flood  plain  discharge  calculations 
and  thus  backwater  calculations. 

A  review  of  the  literature  demonstrates  there  are  two  basic  methods 
in  which  the  estimation  of  backwater  has  advanced.  The  first,  and  by  far  the 
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most  rapid,  is  the  physical  modeling  approach.  The  second  method  is  an  empirical 
approach  in  which  laboratory  analysis  techniques  are  applied  to  field  studies 
and  the  results  are  used  to  modify  the  flume-based  theories  of  Kindsvater  and 
Liu.  As  the  effect  of  flooding  on  ever-urbanizing  areas  and  the  increased  cost 
of  construction  and  maintenance  of  crossing  structures  has  escalated,  so  has  the 
interest  in  more  reliable  methods  of  estimating  backwater.  The  present  study 
recognizes  this  and  provides  a  large  data  base  for  the  current  needs  of  this 
study  as  stated  in  the  purposes  for  the  study  and  for  the  needs  of  both  present 
and  future  proposed  backwater  techniques. 
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III.  HYDRAULIC  MODEL  DESIGN 


GENERAL  DESCRIPTION 

The  difficulties  involved  in  measuring  flows  and  velocities  during 
major  flood  events  make  the  study  of  hydraulic  behavior  in  the  field  an 
expensive  if  not  impossible  task.  For  this  reason,  hydraulic  engineers  have 
historically  used  physical  models  of  flood  plain  behavior  in  order  to  exercise 
control  over  the  events  taking  place  and  to  be  able  to  vary  parameters  in  order 
to  determine  important  variables  and  establish  functional  relationships.  If 
dynamic  similitude  is  maintained  between  the  model  and  the  .prototype,  results 
are  directly  applicable  to  the  field  situation.  It  must  be  recognized  of  course 
that  the  wide  variation  in  width-depth  ratio  observed  in  many  streams  makes  it 
physically  impossible  to  model  the  entire  stream  or  even  the  entire  area 
affected  by  backwater.  However,  the  physical  model  can  be  used  to  study  flow 
patterns  adjacent  to  the  bridge  embankments  and  opening. 

This  chapter  presents  the  analysis  demonstrating  dynamic  similitude 
between  an  average  set  of  prototype  conditions  and  the  model.  Model  distortion, 
which  is  practically  unavoidable  in  studying  any  flood  plain  situation,  is 
developed  and  its  effect  on  model  roughness  quantified.  The  chapter  concludes 
with  a  development  of  the  theory  applicable  to  the  artificial  roughness 
elements  and  configurations  used  in  the  model  studies. 


22 


DYNAMIC  SIMILITUDE 

The  water  surface  elevation  in  an  open-channel  flow  is  governed  by  the 
mutual  interaction  of  forces  of  inertia,  gravity  and  viscosity.  Due  to  the 
presence  of  a  free  surface  in  open-channel  flow,  the  gravity  forces  are  always 
important.   In  some  stretches  of  river  where  the  water  surface  slope  is  mild, 
the  viscous  effect  may  become  important  and  both  gravity  and  viscous  forces 
control  the  flow. 

The  effect  of  gravity  is  characterized  by  the  Froude  number: 

IF  -% 

gR 

where 

V  =  average  velocity 

R  =  hydraulic  radius 

g  =  gravitational  constant 

and  the  viscous  effect  by  the  Reynolds  number: 

R     VRp  _  VR 
e      y      v 

where      R   =  Reynolds  number, 
e 

y  =  dynamic  viscosity, 

p  =  density,  and 

v  =  kinematic  viscosity. 

When  a  scale  model  is  used  to  study  the  flow  characteristics  of  its 
prototype,  dynamic  similarity  requires  that  the  Froude  number  and  Reynolds 
number  be  the  same  in  both  model  and  prototype.  That  is, 

V  R     V  R              V  2     V n2 
m  m     p  p  _m _  _p_ 

v       v  q  R      g  R 

m       p  ym  m     ap  p 
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Under  ordinary  conditions,  g  =  g  ,  hence  to  satisfy  both  the  above  conditions 


where  the  subscripts  m  and  p  refer  to  the  model  and  prototype,  respectively. 

Under  ordinary  conditions,  g  =  g  ,  her 
j  3m   3p 

simultaneously  it  is  necessary  to  have 

3/2 


(d 


That  is,  in  order  to  achieve  the  dynamic  similarity  a  scale  model  will  be  such 

that  the  kinematic  viscosity  of  the  fluid  in  the  model  must  be  less  than  that 

in  the  prototype.  Unfortunately,  there  are  few  liquids  less  viscous  than 

water,  and  none  of  those  is  inexpensive  and  safe  to  use.  If  water  is  used 

in  the  model,  i.e.,  v  =  v  ,  then  R  must  equal  R  .  Hence,  the  model  is 

m    p       m      ^     p 

identical  to  the  prototype.  In  other  words,  the  realization  of  a  dynamically 
similar  small  scale  model  is  practically  impossible. 

In  open  channel  flow  the  gravity  forces  and  thus  the  Froude  number  are 
considered  to  be  the  predominant  factors.  Consequently,  the  Froude  number  of 
the  model  and  prototype  should  be  held  identical.  In  prototype  open  channel 
flow,  the  Reynolds  number  is  high,  thus  the  viscous  forces  are  usually  only 
secondary.  In  order  to  account  for  the  viscous  effect  a  scale  model  is  normally 
calibrated  for  its  roughness  by  trial  and  error  process.  Such  a  process 
involves  successive  adjustment  of  model  roughness  until  the  appropriate  flows 
and  depths  in  the  prototype  are  reproduced  in  the  model.  The  model  was  not 
calibrated  in  this  study  because  the  study  is  intended  for  a  spectrum  of 
riverine  situations  rather  than  a  specific  stream  crossing  site.  The  roughness 
field  selected  in  the  study  represents  an  average  condition  over  the  spectrum. 

A  further  complication  frequently  ensues  in  dealing  with  physical  models 
of  river  and  tidal  conditions.  Since  the  laboratory  space  is  usually  limited, 
the  model  scale  in  depth  is  necessarily  small.  The  actual  water  depth  in  the 
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model  may  be  only  a  fraction  of  an  inch,  producing  a  laminar  flow.  Since  the 
flow  in  the  prototype  is  invariably  turbulent,  the  viscous  forces  in  the  laminar 
flow  of  the  model  would  be  of  an  entirely  different  type  compared  with  the 
prototype.  As  a  result,  the  two  systems  would  be  virtually  unrelated  to  one 
another,  and  the  model  results  could  not  be  used  to  represent  the  prototype 
conditions. 

This  difficulty  may  be  overcome  by  the  use  of  different  scales  for  ver- 
tical and  horizontal  dimensions  in  the  model .  The  vertical  (depth)  scale  is 
exaggerated  in  relation  to  the  horizontal  (distance)  scale  to  increase  the 
velocity  scale  and  produce  turbulent  flow  in  the  model.  As  described  by 
Henderson  [11],  the  scale  distortion  does  not  seriously  affect  the  flow  pattern, 
and  in  fixed-bed  models  it  gives  satisfactory  results.  Such  exaggerated  scale 
models  have  the  advantage  of  saving  laboratory  space.  For  example,  if  M  is  the 
ratio  of  the  vertical  scale  to  the  horizontal  scale,  then  there  is  an  M-fold 
economy  for  the  model  space  requirements  compared  to  an  undistorted  model. 

The  background  analysis  leading  to  the  experimental  setup  of  this  study 
is  presented  in  the  following  sections.  It  has  been  assumed  that  the  flow  is 
subcritical,  the  width  of  flood  plan  is  large  in  comparison  with  the  main 
channel,  and  the  channel  geometry  is  nonprismatic.  The  flow  has  been  assumed 
to  be  dominated  by  gravity  and  channel  resistance  except  in  the  immediate 
vicinity  of  the  bridge. 

GOVERNING  FLOW  EQUATIONS  IN  BRIDGE  WATERWAYS 

Consider  a  gradually-varied  flow  in  an  open  channel.  The  total  energy 
at  a  section  x  is  given  as 

H  =y  +  |^  (8) 

2g 
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where 

y  =  water  surface  elevation, 

V  =  average  velocity, 

a  =  energy  coefficient  and 

g  =  gravitational  acceleration. 

Differentiating  Equation  8  with  respect  to  x  yields 


d_H  =  dy  +  o_  d_V2 
dx    dx   2g  dx 


(9) 


Since      y  =  h  +  z 


dy  =  dh  dz  =  dh  nox 

dx    dx   dx   dx   b  [W) 


where 


h  =  flow  depth, 

z  =  bottom  elevation,  and 

s  =  slope  of  the  channel  bed. 

If  the  flow  rate  of  the  channel  is  Q,  and  the  cross-sectional  area  at  x  is  A, 

v  -jf  (id 

The  second  term  in  Equation  9  can  be  expressed  as 


o_  dV2=  oQ2dA22=   oQ2  dA 
2g  dx  "  2g  dx    "  gA3  dx 


(12) 
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In  the  nonprismatic  channel,  A  is  a  function  of  x  and  h.  In  functional  form 
A  =  cf>  (x,  h).  Thus, 


dA    9A   9A  dh   9A  ,  R  dh  M-» 

dx  ~  3x   9h  dx   ax   b  dx  [U) 


where 


B  =  width  of  water  surface  at  depth  h. 
Substituting  Equation  13  into  Equation  12,  we  obtain 


dV2   =   _aQi/8A  +  B 


dh 


2g-dx-    -  -  g-^^^dx";  (14) 

By  substituting  Equations  10  and  14  into  Equation  9,  we  obtain 


\Jt+  B  d7J 


dH     _  dh       _       a  Q2 

dx"    "  H5T-  s  "  g^7  l^'       ^ 


Since  the  total  energy,  H,  decreases  when  x  increases,  the  rate  of  change  in 

dH 

-t—  is  always  negative.  Hence  the  energy  slope  (gradient)  is 


At  this  time  the  exact  determination  of  the  energy  gradient  of  a  nonuniform 
flow  is  unknown.  It  is  customary  to  assume  that  the  uniform  flow  formula  and 
roughness  coefficients  may  be  used  to  evaluate  the  energy  gradient  of  gradually 
varied  flow  (where  the  energy  loss  may  be  regarded  as  entirely  due  to  surface 
resistance).  At  a  given  section,  when  the  Manning  formula  is  used, 

n2V2 

S/r  = 


f   2.22R^3  (17) 
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and,  when  the  Chezy  formula  is  used, 

v2 
s^   =  V  (18) 


>f     C2R 


where, 


n  =  Manning's  roughness  coefficient, 
C  =  Chezy  roughness  coefficient,  and 
R  =  hydraulic  radius. 

Determination  of  the  energy  gradient  is  further  complicated  if  the  flow 
is  rapidly  varied.  In  this  case  the  overall  flow  resistance  consists  of  boun- 
dary friction  and  local  energy  loss  caused  by  the  abrupt  change  in  flow  geometry. 
In  the  case  of  bridge  waterways  the  geometric  change  involves  the  flow  con- 
traction and  expansion  at  the  bridge.  The  local  energy  loss  is  usually  treated 
as  proportional  to  the  velocity  head  at  a  given  section.  For  a  particular 
geometric  variation,  i  (e.g.,  channel  expansion),  the  energy  loss  is  expressed 
as: 

(HL).   =  K.  £  (19) 

where  K.  is  the  loss  coefficient.  The  total  energy  loss  due  to  N  geometric 
variations  existing  in  the  river  reach  L  is: 


a  ''j/ill^l-  (20) 


The  average  energy  gradient  of  reach  L  is 


^.J<Y£  (21) 

L    L  2g  '  {dl} 
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Hence,  for  rapidly-varied  flow,  the  value  of  sf  is  given  by  the  sum  of 
frictional  losses  and  the  loss  incurred  due  to  the  changes  in  the  stream 
geometry.  Thus 


or 


dH 
dx 


dH 
dx 


n2V2      K  V2 

S,-  =  n—  + 


2.22R'3 


c   ■  V2  ,  K  V2 
Sf  "  ,,„  *  L  2g 


L  2g 


C2R 


(22) 


(23) 


Equation  23  can  be  rearranged  to  yield: 


dH 
dx 


{***f)S 


(24) 


Substituting  Equation  24  into  Equation  15,  we  obtain 


dh 

dx 


(g         Kr\  V2   a  Q2  (  9A   R  dh  \ 
\F      K)    ^+gF  \  9^  +  B  d7  / 


(25) 


Applying  A  =  RP  and  f  =   V,  where  P  is  the  wetted  perimeter,  to  Equation  25 
we  have 


where 


1  dh 

s  dx 


1 


IF2 


g    KR 
X?+   2L 


)■ 


1  -  #  HF2 


uF  = 


gR 


a  3A_ 
P  9X 


(26) 


Equation  26  is  the  differential  equation  for  one-dimensional  non- 
uniform flow.  The  flow  could  be  either  gradually  varied  or  rapidly  varied. 
The  purpose  of  deriving  Equation  26  is  to  use  the  equation  for  establishing 
the  similarity  requirements  for  the  scale  model,  as  will  be  described  in  the 
following  section. 
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REQUIREMENTS  OF  MODEL  SCALES 


Equation  26  may  be  written  for  the  prototype  as 


1  dhp 

1 

..V 
sp 

L\  p           P  / 

9An" 

a           p 

P      8X 
P        P. 

spdxp 

1 

aB„ 

(27) 


The  model  version  of  Equation  27  is: 


,  dh 
1 m 

s  dx 

m   m 


1TF  2  T/     K  R  \ 

m   |  g  ,  mm  I 

m  L  \  m     m  / 


A    1 

a 

m 

P 

X 

in 

m  J 

1  - 


aB 

n-^F  2 
P_   m 


(28) 


Substituting  a  a^  for  a  in  Equation  28,  we  obtain 
r  p     m 


i  dh  /  h  \ 
] P-[-J^-l  = 

s  dx  \s  x  / 

p   p  \  r  r/ 


SP   [VW7)    2LP   V     srLr  /] 


aB  /B  HF   z 

J2.  „F  2(    r     r 

p  p   \     p 

p  \       r 


m 


i  - 


A     flF  ?A  \ 
a p      r     r] 

PP  sp  srPrxr/ 


(29) 


— 2-  uF  2 


If  the  model    is  expected  to  be  dynamically  similar  to  the  prototype  with 
respect  to  the  section  properties  only,  then  the  prototype  and  model   versions 
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of  Equations  27  and  28  of  the  di mens ion less  equation  of  motion  must  necessarily 
be  identical.  From  a  comparison  of  Equations  27  and  29  it  follows  that  these 
equations  can  be  identical  only  if  each  of  the  expressions  in  the  square 
brackets  of  Equation  29  is  equal  to  unity,  that  is: 


s  x 
r  r 


HF  2K  R 
r  r  r 

srLr 


IF  2A 
r  r 

s  P  x 
r  r  r 


B  uF  2 
r  r 


=  1 


1FF  2 

s~C?  "  ] 
r  r 


=  1 


=  1 


=  1 


(30) 


where  the  subscript  r  represents  the  ratio  of  model  to  prototype. 

If  the  flow  cross  section  is  sufficiently  wide,  it  is  possible  to  express  the 
scale  of  the  wetted  perimeter  P  in  terms  of  the  horizontal  scale  alone.  We 
may  arrive  at  the  following  set  of  scales  of  the  geometric  properties  of  the 
flow  cross  section: 


hr  =  yf  ,  Ar  =  xpyr  ,  Br  =  xr  ,  Pf  =  xr 


(31) 


Substituting   Equation  31    in  Equation  30  yields 


31 


br    — 

A  y* 

r  (32) 

Kr  =  1 

(32) 

C  2   =  — 

Equation  32  gives  the  conditions  of  dynamic  similarity  for  nonuniform 
flow,  including  gradually  varied  as  well  as  rapidly  varied,  in  distorted  models. 
According  to  Equation  32,  a  vertically  exaggerated  model  must  be  designed  to 
have: 

1.  The  same  Froude  number  as  the  prototype, 

2.  A  larger  slope  than  the  prototype, 

3.  A  larger  roughness  than  the  prototype,  and 

4.  Similar  coefficients  of  local  energy  loss  to  that  of  the 
prototype. 

In  the  design  of  a  scale  model  for  studying  bridge  backwater,  problems 
exist  with  the  third  condition  in  Equation  32  since  loss  coefficients  associated 
with  the  bridge  constriction  have  not  been  established  yet.  Indeed,  these 
coefficients  are  part  of  the  information  sought  in  the  study  itself.  Fortunately, 
coefficients  of  the  local  energy  losses  are  usually  a  function  of  certain 
dimensionless  ratios  describing  the  geometry  responsible  for  the  local  energy 
losses.  It  appears  that  these  ratios  are  often  identical  for  both  model  and 
prototype,  in  spite  of  the  length  distortion.  Consequently,  the  third  condition 
(i.e.,  K  =  1)  in  Equation  32  can  frequently  be  regarded  as  satisfied 
automatically. 

xr 
In  specifying  the  roughness  requirement,  C  2  =  w—  in  Equation  32, 

roughness  factors  may  be  expressed  by  Manning's  n  as  well  as  Chezy's  C. 

The  expression  for  the  n  ratio  in  the  model  and  the  prototype  is  derived  as 

follows: 
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Since 
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» 
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Equation  33  then  yields 


1/     yr 

n  2  =  r  A   .  JL 
r  r  x„ 


or 


r/6  y/2 
'r   ="VW 


(33) 


(34) 


Application  of  Equation  31   to  Equation  34  yields 

2A 


y, 


n 


J? 


(35) 


Equation  35  can  also  be  derived  by  the  direct  application  of  the 
Manning  equation. 

For  the  practical  purposes  of  this  study,  the  design  of  the  scale  model 
is  based  upon  the  conditions  of: 


nF  =  1  ,  s  =  —  ,  C  = 
r 


y. 


K 


y2/3 

r 
,  or  n  =  — r 


% 


(36 
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Field  data  collected  by  the  U.S.  Geological  Survey  [1]  for  over  one 
hundred  streams  in  the  states  of  Alabama,  Louisiana  and  Mississippi  have 
yielded  the  following  information: 

1.  Width  of  flood  plain:  200  to  9800  ft,  average  2000  ft, 

2.  Depth:  3  to  1 5  ft, 

3.  Bridge  length:  45  to  5900  ft,  average  550  ft, 

4.  Manning  n:  0.03  to  0.2, 

5.  Slope:  0.001  to  0.0018, 

6.  Discharge:  1400  to  205,000  cfs;  average  25,000  cfs. 

From  these  data  the  following  scales  and  physical  model  requirements  were  derived: 

yr  =  =jy        thus,    ym  =  0.25  ft  to  1.25. ft 

\=m  xm  =  20  ft>  bm  =  5-5  ft 

Distortion  =8.33 

s  =  8.33  s  =  0.008  to  0.015 

r  m 

C  -  0.346 
r 

n  =  1.9  n  =  0.057  to  0.38 

r  m 

The  actual  model  was  22.7  ft  wide,  184  ft  long  and  sloped  at  0.0022 
ft/ft.  Discharges  in  the  model  ranged  from  1.0  to  9.6  cfs,  Manning's  "n" 
values  ranged  from  0.05  to  0.36,  measured  depths  from  0.5  to  1.5  ft,  and 
bridge  lengths  from  1.0  to  11  ft  (predominant  range  was  3  to  8  ft).  Experi- 
mental results  should  be  limited  to  those  field  conditions  actually  modeled: 

1.  Flood  plain  width  approximately  2,000  ft. 

2.  Bridge  opening  100  to  1,100  ft. 

=  bridge  opening   = 
flood  plain  width 

It  is  suggested  that  the  flood  plain  width  can  vary  from 
400  to  16,000  ft  if  the  above  ratio  is  within  the  specified 
range. 
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4.  Manning's  "n"  0.026  to  0.189. 

5.  Depth  6  to  18  ft. 

6.  Slope  around  0.0003. 

7.  Discharge  8,000  to  40,000  cfs.  The  discharge  range  is  computed 
from  the  Froude  criterion  where: 

'Fm  "  ,Fp 


*Fr  = 

1   or 

V2 
9Rm" 

V2 

9R  n 

p 

Qi_   - 

Q2 

A2R 
m 

A2Rp 

Qm 

or     Q  =  -  m 
P   V 


The  model  adequately  scaled  all  the  parameters  for  Alabama,  Louisiana  and 
Mississippi  bridge  sites  except  for  slope.  While  it  would  have  been  desirable 
to  run  the  test  at  other  slopes,  such  an  investigation  was  beyond  the  scope  of 
this  study. 

The  applicability  of  the  model  results  can  be  viewed  from  a  different 
perspective  by  using  the  scale  ratios  of  Equation  36.  Since  the  slope  is 
the  primary  limiting  factor,  it  is  possible  to  use  the  slope  ratio  as  a 
controlling  factor.  Assume  that  a  stream  has  a  slope  of  0.0014  and  a  flood 
plain  width  of  1,000  ft.  Then  the  critical  scale  ratios  and  prototype 
limitations  can  be  deduced  as  follows: 


s   : 

s 
m 

0.0022 

=    1 

r 

SP 

SP 

xm 

22 

22 

Ar 

XP 

XP 

1000 

=0.022  -  1_ 


y„  =  s  x,  =  0.0345  =  ~ 

J  r      r  r  29 
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y  % 

nr  -  -~  =  0.715 


Bridge  opening  67  to  495  ft, 
Manning's  "n"  0.036  to  0.25, 
Depth  14.5  to  43  ft. 
Discharge  860  to  3,800  cfs . 


CHANNEL  RESISTANCE 

The  resistance  of  the  prototype  flood  plain  channel  was  simulated  for 
the  purposes  of  this  study  by  an  arrangement  of  artificial  roughness  elements. 
The  determination  of  size,  shape  and  location  of  these  roughness  elements 
required  for  an  accurate  simulation  is  described  in  Volume  I  of  this  study, 

Experimental  Determination  of  Channel  Resistance  for  Large  Scale  Roughness . 
The  process  of  determining  the  characteristics  of  the  roughness  elements  is 
summarized  briefly  below. 

In  the  prototype  condition  of  densely  forested  flood  plains,  the 
energy  losses  of  the  flow  are  due  to  bed  roughness,  bank  roughness  and 
the  resistance  of  bushes,  plants  and  trees  in  the  flood  plains.  These 
roughness  elements  either  are  submerged  or  protrude  through  the  free  surface 
during  floods.  Their  distribution  is  invariably  random,  making  it  impracticable 
to  scale  size  and  distribution  patterns  in  the  model  flume. 

Traditionally,  roughness  element  studies  have  been  conducted  with 
bottom  roughness  elements  that  are  completely  submerged.  The  bridge  backwater 
problem,  on  the  other  hand,  is  influenced  by  trees  and  brush  that  penetrate  the 
water  surface  and  are  spaced  randomly.  It  was  not  considered  feasible  to  use 
model  trees  for  the  experiment,  so  attention  concentrated  on  achieving  various 
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levels  of  channel  resistance  and  relating  that  resistance  to  statistical  repre- 
sentations of  spacing  parameters  where  roughness  elements  are  spaced  randomly 
as  well  as  on  a  regular  pattern. 

The  hydraulic  experiments  for  this  study  were  conducted  in  a  22.7-foot 
wide  flume  with  large  scale  roughness  to  simulate  the  densely  vegetated  flood 
plains  in  the  prototype.  The  roughness  fields  to  be  installed  in  the  large 
flume  were  determined  by  performing  preliminary  testing  and  screening  experiments 
of  a  9-inch  wide  flume.  These  small  flume  experiments  isolated  the  effects  of 
various  shapes  and  densities  of  roughness  patterns. 

Altogether  there  were  seven  different  shapes  of  roughness  element 
tested  in  the  preliminary  study.  Three  types  of  element  distribution,  random,  ' 
rectangular,  and  staggered  (or  diamond),  were  tested  for  flow  rates  ranging  from 
0.1  to  0.6  cfs.  For  each  roughness  element  shape  and  pattern,  the  elements  were 
attached  to  the  channel  bed  and  were  of  sufficient  length  to  protrude  through 
the  water  surface.  The  density  (number  of  elements  per  square  foot  of  the 
channel  bed)  of  the  roughness  elements  was  determined  for  each  configuration. 
All  tests  were  performed  for  a  steady,  nonuniform  flow  condition.  Resistance 
coefficients  for  each  roughness  configuration  were  determined  from  the  test 
data. 

The  experiments  were  performed  using  eight  types  (see  Figure  III-l) 
of  roughness  elements.  All  the  elements  were  approximately  18  inches  in  length. 
Their  projected  widths  varied  from  0.25  inche  to  1.06  inches  and  their 
cross-sectional  areas  varied  from  approximately  0.05  square  inch  to  0.8  square 
inch. 
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The  small  flume  study  established  the  relationships  between  roughness 
and  the  Froude  and  Reynolds  numbers,  depth  to  width  ratio,  element  concentration, 
and  element  pattern.  The  functional  form  is  given  below. 

f  =  *i(F,  R.  £-,  a,  ?) 

where  f  is  the  function  factor. 

The  study  concluded  that: 

f     =   cf>2    (a)  K  =   K 

Essentially  that  for  a  given  roughness  element  pattern,  £,  the  channel 
resistance,  f,  is  a  unique  function  of  the  roughness  concentration,  a- 
Relationships  were  developed  of  the  form: 

t   -  ao 
or 

n   =  al   c    l 

to  quantify  both  the  friction  factor,  f,  and/or  Manning's  n  as  a  function  of 
o.  Values  of  a,  3>  als  and  3X  are  given  in  Table  III-l. 
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Table  III-l.  Values  of  a,  3,  a2  and  gx 


where 


Roughness 
Pattern 

a 

6 

ai 

01 

Random 

5.60 

0.887 

0.208 

0.480 

Rectangular 

3.80 

0.870 

0.183 

0.520 

Diamond 

5.10 

1.045 

0.210 

0.602 

Diamond 

11.40 

0.955 

0.357 

0.706 

Diamond 

7.60 

0.941 

0.275 

0.635 

Diamond 

10.08 

1.333 

0.345 

0.862 

Rectangular 

5.90 

1.360 

0.215 

0.824 

Diamond 

6.03 

0.980 

0.238 

0.603 

Diamond 

7.40 

0.735 

0.293 

0.497 

Diamond 

4.80 

0.897 

0.247 

0.598 

Diamond 

4.40 

0.837 

0.226 

0.565 

Diamond 

11.50 

1.007 

0.374 

0.673 

Diamond 

13.30 

0.992 

0.374 

0.589 

f  =  a  a  ,  and 


pi 
n  =  a!  a  1 


As  stated  previously  the  objective  of  this  study  was  to  determine  the 
roughness  patterns  to  be  placed  in  a  large  test  flume  to  produce  sufficiently 
high  resistance  to  characterize  the  flow  field  in  heavily  vegetated  flood  plains 
The  selection  of  such  roughness  patterns  was  governed  by  (1)  ease  of  installing 
the  roughness  elements,  (2)  degree  of  roughness  in  the  prototype  flood  plain, 
and  (3)  scale  of  the  model. 

The  bridge  backwater  experiments  conducted  in  the  large  flume  were 
not  intended  to  represent  the  hydrualics  of  any  site-specific  case;  rather,  they 


40 


cover  a  wide  range  of  typical  hydraulic  characteristics  of  bridge  crossing 
sites.  The  kind  of  roughness  field  to  be  installed  in  the  large  flume  thus 
must  provide  sufficient  range  of  variation  to  characterize  the  field  conditions. 
As  stated  previously,  roughness  coefficients  show  a  range  of  between  0.03  to 
0.2  for  the  flood  plains. 

The  range  of  n  values  (0.057  to  0.38)  required  for  the  large  test  flume 
was  used  to  find  a  roughness  pattern  frum  the  small  flume  test  data  which 
satisfied  the  selection  factors  given  at  the  beginning  of  this  section.  For 
a  given  value  of  n,  that  pattern  is  selected  which  gives  the  minimum  number  of 
roughness  elements  required  to  produce  the  specified  roughness  field. 

In  this  manner  it  was  found  that  the  v-shaped  metal  joists  roughness 
elements  in  the  diamond  pattern  yield  the  Manning's  n  values  needed  for  the 
model.  Moreover,  the  metal  roughness  elements  were  considered  the  easiest  to 
install  in  the  large  flume. 
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IV.   EXPERIMENTAL  EQUIPMENT 


FLUME  DESCRIPTION 

All  experimental  runs  were  conducted  in  a  constant  slope-rectangular 
flume  consisting  of  a  concrete  substructure  with  baffles,  weirs,  pumps,  and 
tailgate  located  as  shown  in  Figure  IV-1 .  Dimensions  of  the  flume  were  22.7 
ft  wide,  184  ft  long,  and  a  maximum  depth  of  3.0  ft  with  a  bottom  slope  of 
0.0022.  At  the  upstream  end,  pumps  were  installed  to  provide  a  variable  flow 
of  water  from  a  reservoir  situated  on  the  level  below  the  flume.  The  water  was 
pumped  from  the  reservoir  and  discharged  through  the  palette  diffusers  into  a 
stilling  pond  bounded  upstream  by  the  flume  wall  and  downstream  by  a  rock  and 
wire  baffle.  The  baffle  was  designed  for  through-flow  and  was  never  used  as 
a  weir.  This  arrangement  was  not  adequate  to  distribute  the  flow,  so  additional 
baffles  were  added  after  the  initial  runs  were  made. 

The  stone  and  wire  baffle  formed  the  upstream  boundary  for  a  second 
stilling  pond,  which  in  turn  was  bounded  downstream  by  a  wood  and  metal -pi ate 
weir.  Downstream  from  the  weir,  the  water  flowed  over  and  through  a  gravel 
and  wire  baffle,  and  thence  into  the  main  test  section  of  the  flume. 

At  the  downstream  end  of  the  flume,  a  large  wooden  plank  baffle  was 
used  to  regulate  flow  out  of  the  flume.  This  wooden  baffle  or  tailgate  was 
calibrated  to  produce  uniform  flow  for  each  flow  and  roughness  configuration 
without  the  bridge  model  in  place.  During  test  runs  with  the  bridge  models  in 
place,  the  downstream  baffle  was  set  at  proper  calibrated  openings. 

An  arbitrary  reference  station  was  established  at  the  longtudinal  28  ft 
11  in.  mark  of  the  flume.  This  was  designated  as  90  ft.  Relative  to  this  value, 
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the  crossing  centerline  (centerline  of  the  simulated  crossing)  was  at  150  ft 
and  the  tailgate  at  225  ft.  The  model  embankments  had  a  top  width  of  1  ft  and 
a  bottom  width  of  2  ft  with  sloping  embankments  extending  6  in.  at  the  bottom 
from  each  side.  Embankments  were  constructed  in  4-ft  modular  sections  to  allow 
the  opening  width  to  be  easily  varied.  The  bottom  of  the  flume  was  lined  with 
4  ft  x  8  ft  x  3/4  in.  marine  plywood  to  provide  a  surface  for  attaching  rough- 
ness elements.  Extending  1  ft  into  the  flume  on  north  and  south  walls  were 
energy  absorbing  mats  of  rubberized  horsehair  matting.  This  allowed  a  data 
gathering  area  of  20.7  ft  x  125  ft.  Both  the  weir  and  subsequently  a  manometer 
were  calibrated  with  a  volumetric  tank  (not  shown  in  Figure  IV-1)  located  on 
the  level  below  the  flume  and  connected  to  the  flume  by  a  flow  bypass. 

A  rectangle  grid  system  was  defined  and  marked  on  the  flume.  A  point  just 
beyond  the  gravel  wire  baffle  against  the  energy  absorbing  matting  on  the  north 
wall  was  given  the  coordinate  (90.0,  0.0).  The  first  coordinate  of  the  pair,  x, 
is  in  the  direction  of  flow  and  ranged  between  90.0  and  215  ft.  The  center  of 
the  crossing  structures  (embankments)  was  located  at  x  =  150.0  ft.  The  second, 
or  Y  coordinate,  ranged  between  0.0  and  20.7  ft.  The  centerline  was  located  at 
Y  =  10.3  marker.  All  data  collected  is  referenced  using  this  coordinate  system. 
In  some  instances  the  X  coordinate  is  translated  by  90  ft  to  give  a  (0.0,  0.0) 
reference  point.  This  causes  the  X  coordinate  to  range  between  0.0  and  125.0. 
When  this  occurs  it  is  either  stated  in  the  text  or  obvious  from  the  context  of 
the  material . 


ROUGHNESS  ELEMENTS 

A  primary  objective  of  this  study  is  to  quantify  the  relationship 
between  backwater  and  roughness.  Roughness  in  the  large  flume  was  simulated 
by  bolting  V-shaped  metal  joists  to  the  3/4  in.  marine  plywood  sheets  which 
were  anchored  to  the  concrete  floor  of  the  flume.  Each  roughness  element  was 
18  in.  high  and  extended  through  the  free  surface.  This  was  purposely  done  to 
achieve  the  conditions  noted  on  heavily  vegetated  flood  plans  where  often 
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trees  are  not  inundated  by  the  flood  waters.  No  attempt  was  made  to 
approach  geometric  similitude  in  the  flood  plain  with  the  arrangement  of 
roughness  elements.  However,  since  previous  flume  studies  had  underpredicted 
backwater,  possibly  by  not  including  an  accurate  representation  of  roughness, 
part  of  the  goals  of  this  study  was  to  estimate  their  effect. 

It  is  important  to  note  that  the  bridge  backwater  experiments 
conducted  in  the  large  flume  were  not  intended  to  represent  the  hydraulics  of 
any  site-specific  case.  Instead,  they  cover  a  wide  range  of  typical  hydraulic 
characteristics  of  bridge  crossing  sites.  As  reported  in  Chapter  III,  the 
prototype  range  of  Manning's  n  was  between  0.03  and  0.2  for  flood  plains. 
The  exaggeration  in  vertical  scale  (distortion)  of  the  model ,  with  the 
dynamic  similitude  derivation  given  in  Chapter  III,  requires  roughness 
in  the  large  flume  model  to  range  between  0.06  and  0.4.  The  development  of 
the  shape  and  location  of  the  elements,  briefly  summarized  in  the  preceding 
chapter,  is  the  subject  of  Volume  I  of  this  report,  Experimental  Determination 
of  Channel  Resistance  for  Large  Seale  Roughness,   Report  No.  FHWA-RD-75-71 . 
Figure  IV-2  shows  actual  roughness  patterns  that  were  used  in  the  large 
flume.  The  number  in  percent  for  the  various  roughness  patterns  represents 
the  percentage  density  compared  to  the  reference  pattern  shown  in  Figure  IV-2, 
The  various  patterns  will  be  identified  hereafter  in  this  report  simply  by 
these  percentages  and/or  a  corresponding  letter  code. 

The  combination  of  two  separate  densities  was  used  to  simulate  flood 
plain  and  main  channel  conditions  for  centered  crossings  only.  The  main 
channel  was  defined  as  that  area  extending  1.5  ft  on  either  side  of  the  model 
centerline  (10.3)  and  extending  the  length  of  the  flume. 

Table  IV-1  lists  the  roughness  combinations  used  in  this  study. 
Figure  IV-3  shows  the  measured  Manning's  n  as  observed  in  the  flume  for  each 
of  the  roughness  combinations  as  a  function  of  flow.  These  calculations 
were  made  using  a  flume  slope  of  0.0022,  normal  runs,  and  the  overall  average 
normal  depth. 
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TABLE  IV-1 
ROUGHNESS  COMBINATIONS 


Roughness  Density 
Expressed  as  Percent  of  Density 


Letter 

Designation 

Fl 

ood  PI 

ain 

Ma 

in  Channel 

P 

100 

100 

A 

100 

50 

B 

75 

75 

C 

75 

25 

D 

50 

50 

*See  Figure  IV-2 
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.35- 
.25- 

.15- 
.05 


ROUGHNESS  P 


1 1 1 1— 

1.0       2.0       3.0      4.0 


-i 1 1— 

7.0       8.0      9.0 


5.0      6.0 

Q  (cfs) 


ROUGHNESS  A 


.25- 


i r 

8.0       9.0 


Q  (CFS) 


<o    "       " 


MANNING'S   n     FOR    ROUGHNESS    PATTERNS  AS  A  FUNCTION    OF 

FLOW 


FIGURE  IV-3 
(cont'd.) 


48 


.35- 


.25- 


15- 


.05. 


ROUGHNESS  B 


~1 1 1 r— 

1.0       2<0       3.0      4.0 


~i 1 r- 

5.0      6.0       7.0 
Q   (CFS) 


~l r~ 

8.0      9.0 


.35- 


.25- 


5- 


.05- 


ROUGHNESS  C 


-i 1 1 -|- 

1.0       2.0      3.0      4.0 


-i 1- 

5.0      4.0 

Q   (CFS) 


-I 1 1— 

70      8.0      9.0 
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DESCRIPTION  OF  CROSSING  CONFIGURATIONS 

Each  crossing  configuration  consisted  of  an  interior  4  ft  x  1  ft  x 
1-2/3  ft  metal  frame  constructed  of  steel  angles  covered  with  1/2  in.  marine 
plywood.  Wood  surfaces,  both  exterior  and  interior,  were  painted  with  two 
or  more  coats  of  high  quality  white  enamel  paint  to  prevent  their  deterioration 
in  water.  The  wooden  embankments  and  abutment  shapes  were  connected  to  the 
metal  frame  by  means  of  8  in.  toggle  bolts,  which  did  not  penetrate  the 
exterior  surface  of  the  marine  plywood.  Marine  plywood  extensions  were  appended 
to  the  metal  frame  on  both  the  upstream  and  downstream  faces  to  allow  easy 
adjustment  of  the  opening  size.  The  width  of  each  model  was  2  ft  measured  at 
the  toe  of  the  upstream  and  downstream  embankments  to  1  ft  at  the  top.  All  of 
the  models  had  the  same  height  of  20  in. 

The  center  of  each  crossing  configuration,  with  the  exception  of 
skewed,  was  located  at  the  150  ft  marker  in  the  flume.  The  toe  of  the  upstream 
and  downstream  embankments  were  located  at  149.0  and  151.0  respectively.  This 
provided  approximately  60  ft  upstream  and  75  ft  downstream  for  measurement. 
Skewed  configurations  had  their  centers  at  150.0  and  formed  an  angle  at  either 
15  degrees  or  45  degrees  with  respect  to  the  flow  direction. 

All  joints  of  the  crossing  configuration  were  well  caulked.  Care  was 
taken  to  maintain  the  roughness  characteristics  of  the  exterior  surface. 
Caulking  was  applied  around  the  upstream  base  of  the  embankments  at  the 
junction  of  the  wooden  appendages  to  the  wood  covered  metal  frame,  and  at  the 
abutment  face  to  the  frame.  The  hollow  crossing  configuration  was  filled  with 
water- 
Each  crossing  configuration  was  inserted  into  a  1  ft  wide  by  3/4  in. 
deep  cut  in  the  marine  plywood  sheets  covering  the  floor  of  the  flume.  This 
cut  extended  completely  across  the  flume.  The  roughness  elements  were  mounted 
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on  these  4  ft  x  8  ft  sheets.  After  the  crossings  were  put  in  place,  a  piece  of 
the  plywood  was  laid  across  the  cut  in  the  opening  to  keep  the  floor  even.  No 
roughness  element  was  placed  on  this  2-ft  strip  across  the  opening. 

WINGWALL 

The  abutment  shape  used  most  extensively  in  the  testing  was  the  45  degree 
wingwall  shown  in  Figure  IV-4.  Figure  IV-4  shows  the  as-constructed  dimensions 
of  the  embankments.  Embankment  sides  were  attached  to  the  upstream  and  down- 
stream faces  of  the  crossing  extending  from  the  edge  of  the  flume  to  the  abut- 
ment face. 

SPILLTHROUGH 

A  second  abutment  face,  the  spillthrough,  was  also  tested  in  the  flume. 
The  face  of  the  abutment  is  tapered  from  the  top  to  bottom  and  thus  makes  a 
trapezoidal  opening  shape.  Each  face  of  the  abutment  is  sloped  so  that  the  toe 
protrudes  6  in.  into  the  constriction,  or,  for  two  abutments,  1  ft.  Opening 
size  is  measured  from  the  toe  of  the  embankment  but  of  course  varies  with  depth. 
Embankments  attached  to  the  upstream  and  downstream  faces  of  the  crossing  were 
the  same  as  for  the  wingwall.  Figure  IV-5  provides  a  sketch  of  the  spillthrough 
abutment  and  pertinent  dimensions. 

SPUR  DIKES 

An  upstream  elliptical  spur  dike  was  the  third  shape  tested  in  the  large 
flume.  The  elliptic  spur  dikes  consisted  of  40  layers  of  1/2  in.  marine  plywood, 
18  in.  wide  at  the  bottom  and  6  in.  at  the  top,  for  an  embankment  ratio  of 
1.6:1  (horizontal/vertical).  The  elliptic  dikes  constructed  have  a  major  axis 
of  4.85  ft  with  a  ratio  of  major  axis  to  minor  axis  of  2.5:1.  Figure  IV-6 
shows  the  elliptical  spur  dike  used  in  all  experiments. 
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The  equation  describing  the  construction  is 


Y2   X2     , 

a2   b2  [      } 


where 

Y  =  coordinate  in  the  direction  of  flow, 

X  =  coordinate  perpendicular  to  the  direction  of  flow, 

a  =  major  axis, 

b  =  minor  axis,  and 

a/b  =  2.5. 

In  a  late  phase  of  the  study,  logarithmic  spur  dikes  were  tested. 
In  these  tests,  the  opening  was  moved  to  one  side  of  the  flume  so  that  the 
flume  served  essentially  as  a  symmetrical  half  of  a  larger  channel  which, 
in  effect,  doubled  the  discharge  and  horizontal  dimensions  of  the  model. 
Detailed  velocities  and  depths  were  recorded  for  flow  around  the  spur  dikes. 
These  data  are  summarized  in  Volume  V,  Data  Report  for  Spur  Dike  Experiments, 
but  their  analysis  was  beyond  the  scope  of  the  present  contract. 
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V.  EXPERIMENTAL  PROCEDURE 


GENERAL  DESCRIPTION 

The  flume  studies  analyzed  in  this  report  were  conducted  from  June 
1972  to  July  1973.  The  site  of  the  flume  was  the  old  National  Bureau  of  Standards 
hydraulics  laboratory  located  at  4200  Connecticut  Avenue,  N.W.,  Washington,  D.C. 
The  flume,  constructed  by  the  Federal  Highway  Administration,  was  operated  under 
the  supervision  of  Water  Resources  Engineers. 

During  the  course  of  the  experiment  three  types*of  crossings  were 

studied  (Figure  V-l).  These  were  centered,  Skewed,  and  eccentric  cross- 
ings. 

Eccentricity  is  defined  as: 

°c  Qa 

1  -  ---   ,  where  Q  <  Q,   or   1  -  jf-   ,  where  Q  <  Qr 
Qa        c   a  0C        a    c 

Flow,  or  Q,  is  defined  as  the  flow  occurring  over  the  cross  section  of  the 
flood  plain  obstructed  by  the  roadway  embankments  (Figure  V-2) . 
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Three  abutment  shapes  were  studied  in  conjunction  with  the  centered 
crossing.  These  were  45  degree  wingwall,  spillthrough,  and  upstream  elliptical 
spur  dikes  as  shown  on  Figure  V-3.  All  tests  conducted  for  skewed  and  eccentric 
crossings  used  wingwall  abutments. 

An  important  aspect  of  this  experiment,  previously  described  in  Chapter 
III,  was  the  simulation  of  roughness  by  the  use  of  artificial  roughness  elements 
penetrating  the  free  surface.  Roughness  elements  were  installed  and  checked  in 
both  the  flood  plain  area  and  main  channel.  Ten  separate  patterns  were  used  in 
the  flume  which  gave  an  effective  Manning's  n  ranging  from  0.06  to  0.4.  These 
values  are  necessarily  higher  than  those  observed  in  nature  in  order  to  compen- 
sate for  model  distortion.  All  eccentric  and  skewed  crossing  configurations 
used  the  same  roughness  pattern  (type  D  -  50  percent  flood  plain  and  main 
channel).  Centered  crossings  accounted  for  all  other  roughness  configurations. 
A  description  of  the  roughness  patterns  related  to  the  model  runs  is  presented 
in  following  sections. 

In  general,  bridge  openings  for  all  flume  runs  were  set  at  1,  3,  5,  8, 
or  11  ft,  respectively.  Few  flume  runs  were  made  with  the  11 -ft  opening  due  to 
the  difficulty  in  locating  and  measuring  the  maximum  backwater.  A  5-ft  opening 
modeled  the  average  prototype  opening  on  a  1:100  scale. 

Flows  for  model  runs  varied  between  1.0  and  9.6  cfs.  Typically  flows 
in  the  neighborhood  of  2.0,  5.0,  7.0  and  9.0  cfs  were  used  for  any  set  of 
roughness  configurations,  opening  size,  abutment  type,  and  crossing  type. 
Calculations  for  dynamic  similitude  indicated  a  flow  of  approximately  6.0  cfs 
is  needed  to  satisfy  average  prototype  flow  conditions  of  25,000  cfs. 

Operating  procedures  for  ewery   series' of  runs  were  standardized  and 
carefully  controlled.  Flows  that  were  to  be  used  for  the  entire  series  of 
openings  and  models  were  tested  in  the  flume  with  roughness  patterns  in  both 
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the  flood  plain  and  main  channel  but  without  any  embankment  construction.  This 
permitted  the  calculation  of  normal  depth  for  each  flow  at  all  points  in  the 
flume.  All  the  different  flows  used  in  the  experiment  were  measured  in  a  volu- 
metric tank  of  known  cross-sectional  area  and  measurable  depth.  Simultaneously, 
all  flow  entering  the  flume  passed  over  an  annubar  connected  to  a  manometer 
which  was  then  calibrated  against  the  flow  measured  in  the  volumetric  tank. 
This  allowed  the  reproduction  of  flows  and  eliminated  the  need  for  the  constant 
use  of  the  volumetric  tank. 

After  the  normal  depth  of  each  flow  was  determined,  an  abutment  shape 
and  model  embankments  were  installed  in  the  flume  and  the  flows  to  be  tested 
were  run.  Typically,  after  flow  was  set  by  use  of  the  manometer,  a  period  of 
one  to  two  hours  was  required  by  recording  centerline  elevations  along  the 
entire  flume  every   15  minutes  until  the  same  values  were  recorded  in  successive 
measurements.  When  this  occurred,  the  entire  flume  was  considered  to  be  at 
equilibrium  and  data  were  collected  over  the  flume  for  that  run. 

The  determination  of  water  surface  elevation  required  the  measurement 
of  the  elevations  of  the  floor  of  the  flume  for  each  of  the  approximately  one 
hundred  depth  recording  stations.  The  floor  elevation  plus  depth  gives  the 
corresponding  water  surface  elevation.  Since  this  fairly  simple  procedure  was 
complicated  by  the  3/4  in.  marine  plywood  installed  over  the  reinforced  concrete 
floor,  precautions  were  taken  to  minimize  error  in  the  measurements.  First,  at 
the  start  of  daily  operations  the  flume  was  allowed  to  become  thoroughly  satu- 
rated before  any  data  were  taken.  Second,  bottom  elevation  measurements  were 
taken  on  a  weekly  basis  with  a  flow  representing  an  average  depth  in  the  flume 
to  check  for  changes  in  the  floor  due  to  warping  or  other  causes.  As  the 
experiment  progressed,  the  measurements  of  bottom  elevations  were  needed  less 
frequently. 
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Each  run  required  four  technicians  and  a  supervising  engineer.  Actual 
data  collection  usually  required  two  hours  after  the  steady  state  condition  was 
established. 


HYDRAULIC  MODEL  TESTS 

During  the  period  in  which  the  hydraulic  model  was  in  operation,  a  total 
of  285  roughness,  flow,  abutment,  and  opening  combination  runs  were  tested. 
Table  V-l  summarizes  the  test  matrix  for  the  roughness  pattern  types  and 
crossing  configurations. 


TABLE 

V-l 

SUMMARY 

OF  EXPERIMENTAL 

RUNS 

Crossing  Con 

figurations 

Centered 

Skew 

Eccen. 

Normal 

Roughness 
Type 

Wing- 
wall 

Spill- 
through 

Spur 
Dike 

Wing- 

wal  1 

Wing- 
wall 

Total 
Runs 

P 

12 

4 

-- 

-- 

-- 

6 

22 

A 

8 

12 

7 

-- 

-- 

6 

33 

B 

9 

27 

12 

-- 

-- 

6 

54 

C 

12 

12 

12 

-- 

-- 

4 

40 

D 

9 

12 

11 

38 

12 

4 

86 

Total 

50 

67 

42 

38 

12 

26 

285 

Both  skew  and  eccentric  crossing  configurations  used  wingwall  abutment 
shapes.  The  elliptical  spur  dikes  were  used  only  on  the  upstream  side  of  the 
embankment.  Spur  dike  data  collected  using  logarithmic  spur  dikes  are  not 
included  in  the  table. 
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Tables  V-2  through  V-5  present  a  complete  tabulation  of  all  the  flume 
runs  accomplished.  Table  V-2  contains  a  listing  of  the  centered  opening  runs 
subdivided  by  roughness  element  pattern.  Flow,  further  subdivided  by  either 
abutment  shape  or  spur  dike,  is  condensed  to  the  nearest  0.1  cfs  for  the  purpose 
of  the  presentation.  Tables  V-3  and  V-4  give  eccentric  and  skew  runs,  respec- 
tively, while  Table  V-5  lists  the  flows  used  for  normal  runs. 

All  the  runs  with  the  exception  of  the  fourteen  shown  in  Table  V-6  had 

roughness  elements  excluded  from  the  opening.  These  fourteen  runs  were  not 

used  in  any  analysis  presented  in  Chapter  VI  and  are  included  here  only  for 
completeness. 


DATA  COLLECTION 

Data  collected  included  primarily  water  surface  and  bottom  readings, 
and  velocity  and  direction  of  flow.  Additional  depth  measurements  were  taken 
at  centerline  positions  and  in  the  area  between  the  location  of  maximum  backwater 
and  the  embankments  of  the  crossing  model.  Depth  measurements  were  taken  both 
60  ft  upstream  and  65  ft  downstream  from  the  crossing  model  section.  Velocity 
measurements  were  primarily  confined  to  that  area  between  the  cross  sections 
defined  by  the  location  of  the  maximum  backwater  and  the  reattachment  points 
as  recorded  on  the  flume  walls.  The  reattachment  points  are  the  two  locations 
where  the  flow  lines  first  strike  the  flume  walls  downstream  of  the  constriction. 

Data  gathered  during  each  run  included  those  parameters  listed  in 
Table  V-7.  All  the  data  collected  have  been  computerized  and  will  become 
property  of  the  FHWA.  The  analysis,  as  presented  in  this  report,  required  the 
development  of  special  in-house  software,  which  is  also  the  property  of  the  FHWA. 
The  following  sections  describe  the  procedure  for  collecting  velocity,  depth  and 
flow  measurements.  Analysis  of  these  data  and  pertinent  calculations  using  raw 
data  are  presented  in  Chapter  VI. 
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TABLE  V-2 

EXPERIMENTAL  INVESTIGATION  CHART 

CENTERED  OPENING  RUNS 


Opening    Flow  Rate    Abutment 

Wingwall   Spill  through   Spur  Dike 


Size 
1ft) (cfs) 


ROUGHNESS  PATTERN  P 

( 1 00% 

FP 

,  100%  MC) 

5.2 

X 

2.9 

X 

1.9 

X 

0.9 

X 

7.4 

X 

5.2 

X 

2.9 

X 

1.9 

X 

7.9 

X 

7.4 

X 

5.2 

X 

2.9 

X 

ROUGHNESS  PATTERN  A 
(100%  FP,  50%  MC) 

7.4  x 

5.4  x 

2.7  x           x 

1.9  xx 

1.1  x 

7.4  xxx 

5.4  xxx 

2.7  x           x 

1.9  x 

9.6  x 
7.4  x  x 
5.3  x          x 

2.7  x  x 


NOTE:  FP  =  flood  plain       1  ft  =  0.3048  meters 
MC  =  main  channel      1  cfs  =  0.0283  m^/sec 
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TABLE  V-2  (Cont'd) 
EXPERIMENTAL  INVESTIGATION  CHART 
CENTERED  OPENING  RUNS 


n         rn   n  ,  Abutment 

Opening    Flow  Rate 


Size  (ft)     (cfs)      Wingwall   Spill  through   Spur  Dike 


11 


ROUGHNESS  PATTERN  A  (Cont'd) 


9.8 

X 

7.4 

X 

5.6 

X 

1 

ROUGHNESS 

PATTERN  B 

(75% 

FP, 

75%  MC) 

9.7 

X 

7.4 

X 

X 

5.3 

X 

X 

2.7 

X 

X 

1.9 

X 

9.7 

X 

7.6 

X 

X 

X 

5.4 

X 

X 

X 

2.7 

X 

X 

X 

9.7 

X 

7.4 

X 

X 

X 

5.4 

X 

X 

X 

2.8 

X 

X 

X 

9.7 

X 

X 

7.4 

X 

X 

5.3 

X 

X 

11 


ROUGHNESS  PATTERN  C 
(75%  FP,  25%  MC)  . 

9.7  xxx 

^         7.4  xxx 

6  5.4  x  x  x 

2.7  x  x  x 

NOTE:  FP  =  flood  plain        1  ft  =  0.3048  meters 
MC  =  main  channel       1  cfs  =  0.0283  m3/sec 
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TABLE  V-2  (Cont'd) 
EXPERIMENTAL  INVESTIGATION  CHART 
CENTERED  OPENING  RUNS 


n         ri   n  .  Abutment 

Opening    Flow  Rate 


Size(ft)     (cfs)      Wingwall   Spillthrough   Spur  Dike 


ROUGHNESS  PATTERN  C  (Cont'd) 


9.7 

X 

X 

X 

5 

7.4 

X 

X 

X 

5.3 

X 

X 

X 

2.7 

X 

X 

X 

9.7 

X 

X 

X 

8 

7.4 

X 

X 

X 

5.4 

X 

X 

X 

2.7 

X 

X 

X 

ROUGHNESS 

PATTERN 

D 

(50%  FP, 

50%  MC) 

7.1 

X 

1 

5.3 
3.9 

9.2 
7.1 

X 

X 
X 

X 
X 

X 
X 

3 

6.5 
5.3 

3.8 

9.3 

X 
X 
X 

X 
X 

X 

X 
X 

X 

5 

7.2 

X 

X 

X 

5.3 

X 

X 

3.8 

X 

X 

X 

9.2 

X 

X 

8 

7.2 
5.3 
3.9 

X 

X 

X 
X 
X 

NOTE: 

FP  = 

=  flood 

Pi 

ain 

1 

ft 

= 

0, 

.3048 

meters 

MC  = 

:  main  c 

ha 

nnel 

1 

cf 

s  = 

0, 

,0283 

IT|3 

/sec 
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TABLE  V-3 

EXPERIMENTAL  INVESTIGATION  CHART 

ECCENTRIC  RUNS 
(Wingwall  Abutment) 

Opening     Flow  Rate 

Size(ft)      (cfs)       Eccentricity 


ROUGHNESS  PATTERN  D 
(50%  FP,  50%  MC) 


9.7 

1.0 

6.9 

1.0 

3.6 

1.0 

9.7 

0.87 

6.9 

0.87 

3.6 

0.87 

9.0 

1.0 

6.9 

1.0 

3.6 

1.0 

9.7 

0.81 

6.9 

0.81 

3.5 

0.81 

NOTE:  FP  =  flood  plain        1  ft  =  0.3048  meters 
MC  =  main  channel       1  cfs  =  0.2083  m3/sec 
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TABLE  V-6 

SPILLTHROUGH  RUNS  NOT  ANALYZED 
(All  Roughness  Pattern  B) 


Run  No. 


Opening 
(ft)- 


Flow 
(cfs) 


Comment 


ST20B 

5 

2.8 

ST23B 

5 

5.3 

ST26B 

5 

7.4 

STUB 

8 

2.8 

ST14B 

8 

9.6 

ST27B 

5 

7.4 

ST24B 

5 

5.3 

ST21B 

5 

2.8 

ST18B 

8 

7.4 

ST28B 

5 

7.4 

ST25B 

5 

5.3 

ST22B 

5 

2.8 

ST19B 

8 

7.4 

ST29B 


75%  element  in  opening 
75%  element  in  opening 
75%  element  in  opening 

2  rows  of  elements  in  opening 
2  rows  of  elements  in  opening 

50%  element  in  opening 
50%  element  in  opening 
50%  element  in  opening 

25%  element  in  opening 
25%  element  in  opening 
25%  element  in  opening 


7.4    100%  element  in  opening 


NOTE:  FP  =  flood  plain 
MC  =  main  channel 


1  ft  =  0.3048  meters 
1  cfs  -  0.2083  m3/sec 
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TABLE  V-7 
RAW  DATA  COLLECTED  FOR  EACH  RUN 


Run  number 

Flood  plain  and  main  channel  roughness 

Flow 

Opening  width 

Crossing  type  (angle  of  skew  or  degree  of  eccentricity) 

Abutment  shape 

Vena  contracta  location 

Vena  contracta  width 

Location  of  maximum  backwater 

Flow  reattachment  points 

Water  surface  elevations  (approximately  100/run) 

Velocity  and  angle  of  flow  (approximately  15/run) 
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VELOCITY  MEASUREMENTS 

As  a  standard  procedure,  velocity  measurements  were  taken  between  the 
location  of  the  cross  section  containing  the  maximum  backwater  and  that  contain- 
ing the  reattachment  points.  Measurements  were  also  taken  in  the  opening  at  both 
edges  of  the  vena  contracta.  Measurements  were  taken  using  an  Ott  small  current 
meter.  Each  shaft  was  equipped  with  a  swinging  vane  to  measure  the  angle  of 
flow.  Each  meter  and  shaft  were  calibrated  at  the  beginning  of  the  experi- 
ment by  the  use  of  a  large  towing  tank  designed  for  this  purpose.  Actual 
measurements  were  taken  by  electronically  recording  the  number  of  "clicks" 
caused  by  flow  turning  the  propel  lor  over  a  50-sec  period.  This  was  done  twice 
for  each  measurement  and  the  average  value  recorded.  This  average  value  was 
then  used  to  determine  the  velocity  from  the  calibration  curves  for  the  particu- 
lar shaft  and  current  meter. 

Centered  Models 

The  location  of  velocity  measurements  for  centered  model  runs  were 
a  function  of  both  the  opening  size  and  the  location  of  the  maximum  backwater. 
Figure  V-4  indicates  the  location  of  velocity  measurements  for  all  the  runs 
having  the  crossing  configuration  centered  in  the  flume.  Those  points  Vll, 
V12,  and  V13,  directly  in  the  opening  were  usually  recorded  on  both  edges  of 
the  vena  contracta  and  at  the  centerline.  Point  V3  was  located  on  the  centerline 
at  the  point  of  maximum  backwater.  The  four  points,  VI,  V2,  V4,  and  V5  flanking 
the  centerline  point  on  both  sides  were  taken  to  fully  characterize  the  cross 
section.  This  cross  section  corresponds  to  the  beginning  of  section  1,  zone  of 
contraction  in  Figure  II-l.  Points  numbered  V6  through  V10  were  usually  taken 
at  the  cross  section  determined  by  the  reattachment  points.  This  cross  section 
corresponds  to  section  4,  Figure  II-l,  the  end  of  the  flow  expansion  zone. 
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SKEW  MODELS 

Five  to  seven  velocity  measurements  were  recorded  for  each  skew  run 
depending  on  opening  size.  A  typical  configuration  is  shown  in  Figure  V-5. 
Velocity  measurements  were  recorded  on  centerline  and  directly  behind  each  edge 
of  the  embankment. 

Eccentric  Model 

Velocities  were  recorded  for  the  eccentric  model  runs  across  the  flume 
at  the  maximum  bakcwater,  in  the  opening  at  the  vena  contracta,  and  across  the 
flume  at  the  reattachment  location.  A  typical  configuration  is  given  in 
Figure  V-6.  Locations  extending  across  the  flume  are  equidistant  from  each 
other  with  the  middle  point  on  the  centerline.  These  cross  sections  refer  to 
sections  1,  2  and  4  respectively,  Figure  II-l. 

DEPTH  MEASUREMENT 

By  far  the  most  significant  data  collected  from  the  large  flume  were 
depth  measurements.  A  typical  model  run  recorded  depth  measurements  at  approxi- 
mately 100  separate  locations.  These  locations  can  be  classified  as  either 
stationary  or  moveable.  Moveable  gage  (Table  V-8)  locations  were  a  function  of 
crossing  configuration,  opening  size  and  the  radius  of  maximum  backwater.  These 
gages  were  located  upstream  of  the  constriction  and  usually  within  the  cross 
section  bounded  by  the  radius  of  maximum  backwater.  The  stationary  gages 
numeric  code  and  their  location  are  given  in  Table  V-9.  The  stationary  gages 
were  gages  locked  in  place  and  located  both  upstream  and  downstream  of  the 
crossing.  Figure  V-7  is  a  description  of  the  typical  gage  employed  at  the  flume 
The  smallest  scale  gradation  was  0.001  ft. 

Water  surface  elevation  was  calculated  by  adding  the  bottom  elevation 
to  the  depth  of  the  water  in  the  flume.  Depth  was  calculated  as  the  difference 
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SKEWED  MODEL  VELOCITY    MEASUREMENT  LOCATIONS 

FIGURE  V-5 
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ECCENTRIC  MODEL  VELOCITY    MEASUREMENT  LOCATIONS 

FIGURE  V-6 


76 


up 
o 


t> 

'' 

1    1      1 

*=- 

\ 1 

-ft    1 

<^3 

x 

o 
o 


o 

I- 

GO 


o 

CM 


t> 

1 

1 

1             1 

t=i 

1             1 

H'Pt" 

fe               1 

up 

[Jy. 

ctt. 

1. 

Z 

UJ 


> 


h- 


jllll|llll|,l,l|llll|llll||llljMII|IUIJIU 

.      n      io      *      in 


7^ 


•111 


O 

Itlfl'llifll'.l,! 


y 

"I 

liJ  ^ 


i  7 


Cl< 
<  X 


"l""|" l|iiii|ini|iii I n n|.„,|i,i.| 

(7)   Q   —   cm   ro   <j   m 


o 

CL 


Q 

o 


JFR 


o 


o 


co 


CD 
CM 


UJ 
CD 


2 


< 

CM 


en 

o 
< 


O 
Q_ 


LU 

o 


77 


TABLE  V-8 
NUMBERS  OF  MOVEABLE  GAGES 


42 

58 

56 

63 

77 

55 

74 

60 

69 

82 

80 

75 

62 

76 

88 

67 

57 

81 

59 

68 

61 

70 

TABLE  V-9 
STATIONARY  GAGE  LOCATIONS 


Gage 

X 

Y 

Gage 

X 

Y 

Gage 

X 

Y 

1 

90.0 

10.3 

33 

180.0 

10.3 

87 

140.0 

20.7 

2 

100.0 

10.3 

34 

190.0 

10.3 

89 

149.0 

20.7 

3 

110.0 

10.3 

35 

200.0 

10.3 

90 

151.0 

20.7 

4 

120.0 

10.3 

36 

215.0 

10.3 

94 

170.0 

20.7 

5 

131.5 

10.3 

37 

90.1 

0.0 

95 

18-0.0 

20.7 

6 

135.0 

10.3 

38 

100.0 

1.0 

96 

190.0 

20.7 

7 

137.5 

10.3 

39 

120.0 

0.0 

97 

200.0 

20.7 

8 

140.5 

10.3 

40 

131.5 

11.6 

98 

215.0 

20.7 

9 

142.5 

10.3 

41 

140.0 

0.0 

153 

153.0 

10.3 

10 

143.5 

10.3 

43 

100.0 

0.0 

154 

154.0 

10.3 

11 

144.0 

10.3 

44 

151.0 

0.0 

155 

155.0 

10.3 

12 

144.5 

10.3 

48 

170.0 

0.0 

156 

156.0 

10.3 

13 

145.0 

10.3 

49 

180.0 

0.0 

157 

157.0 

10.3 

14 

145.5 

10.3 

50 

190.0 

0.0 

158 

158.0 

10.3 

15 

146.0 

10.3 

51 

200.0 

0.0 

159 

159.0 

10.3 

16 

146.5 

10.3 

52 

215.0 

0.0 

160 

160.0 

10.3 

17 

147.0 

10.3 

53 

135.0 

5.3 

161 

161.0 

10.3 

18 

147.25 

10.3 

54 

140.0 

5.3 

162 

162.0 

10.3 

19 

148.0 

10.3 

78 

135.0 

15.7 

163 

163.0 

10.3 

20 

149.0 

10.3 

79 

140.0 

15.5 

164 

164.0 

10.3 

21 

150.0 

10.3 

83 

90.0 

20.7 

165 

160.0 

20.3 

22 

151.0 

10.3 

84 

105.0 

18.7 

166 

160.0 

16.4 

23 

152.0 

10.3 

85 

120.0 

20.7 

167 

160.0 

4.3 

32 

170.0 

10.3 

86 

131.5 

19.1 

168 

160.0 

0.3 
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between  the  surface  reading  and  bottom  reading,  plus  the  gage  constant,  which 
was  the  distance  between  the  tip  of  the  point  and  the  tip  of  the  rod  shown  in 
Figure  V-7.  This  depth,  when  added  to  bottom  elevation,  known  previously, 
yielded  surface  elevation. 

FLOW  MEASUREMENTS 

Initially,  the  flow  was  measured  using  annubar  flow  meters  installed 
in  the  discharge  pipes  near  the  control  valves.  These  flow  measurement  devices 
were  found  to  be  insensitive.  The  flow  was  then  measured  using  a  stilling  well 
connected  to  the  upstream  side  of  the  overflow  weir  located  at  the  head  of  the 
flume  (Figure  IV-1 ) .  Depth  in  the  standpipe  was  measured  using  a  hook  gage,  and 
then  calibrated  using  the  volumetric  tank  at  the  base  of  the  flume  to  determine 
flow.  The  volumetric  tank  was  a  large  tank  of  known  constant  cross  sectional 
area  in  which  the  depth  could  be  accurately  determined.  This  permitted  the 
calculation  of  flow  rate  to  two  significant  decimal  places. 

After  three  months  of  operation,  it  was  clear  that  the  stilling  well 
method  did  not  yield  the  required  accuracy  as  flows  could  not  be  duplicated 
exactly.  This,  of  course,  was  due  to  the  large  variation  in  flow  caused  by  a 
slight  change  in  head  over  a  long  weir.  A  light  liquid  manometer  was  subse- 
quently connected  to  the  annubar  flow  meters  placed  earlier  and  an  additional 
valve  installed  near  the  discharge  pipe  outlet.  The  manometer  was  calibrated 
against  the  volumetric  tank  readings.  The  manometer  was  connected  to  a  manifold 
valve  system  attaching  to  the  annubar  meter  in  each  discharge  pipe.  Some 
fluctuation  in  manometer  readings  was  observed  during  the  course  of  the  experi- 
ment. 

This  was  taken  to  mean  that  the  flow  being  observed  did  not  correspond 
to  the  expected  flow,  i.e.,  the  flow  at  which  the  normal  depths  were  measured. 
This  was  partially  solved  by  regressing  flow  on  depth  for  the  normal  runs 
corresponding  to  a  roughness  configuration  and  thus  adjusting  all  the  normal 
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depths  accordingly.  A  detailed  description  of  this  problem  is  found  in 
Chapter  VI  and  the  Appendix. 

QUALITY  CONTROL 

As  described  earlier,  the  bridge  backwater  is  the  incremental  rise 
in  stage  over  the  natural  water  surface  elevation  due  to  constriction  of  flow 
passage  for  a  given  flow  rate.  Thus,  laboratory  determination  of  backwater 
requires  that  a  constant  flow  rate  in  the  test  flume  be  maintained  before  and 
after  the  crossing  models  are  placed.  Since  the  placement  of  crossing  models 
in  the  flowing  water  was  not  feasible,  model  tests  for  the  unconstricted  and 
constricted  conditions  could  be  performed  only  on  separate  occasions.  Thus, 
the  accurate  duplication  of  flow  rates  in  the  flume  during  test  conditions  was 
the  paramount  criterion  for  determining  quality  of  test  data.  Other  factors 
such  as  the  fluctuation  of  voltage  in  the  power  supply  to  the  pumps,  human 
error,  instrumentation  malfunction,  and  change  of  floor  elevations  also  affected 
the  quality  of  data.  Since  some  of  these  errors  are  local  in  nature,  they 
cannot  be  controlled  effectively.  An  inconsistency  in  discharge,  however,  can 
cause  the  stage  to  vary  throughout  the  test  flume,  representing  dissimilar  flow 
regimes.  The  degree  of  duplication,  therefore,  was  the  criterion  used  as  the 
basis  for  assessing  data  quality. 

Unexpected  shifts  in  elevation,  changes  in  expected  point  scatter,  and 
unusual  shifts  in  the  location  or  height  of  the  maximum  backwaters  were  indi- 
cated by  this  method  quickly  enough  to  allow  corrective  action  to  be  taken 
immediately  before  moving  on  to  other  areas  of  the  testing.  This  assured 
reasonable  quality  in  the  large  amounts  of  data  to  be  collected  in  the  short 
period  of  time  the  flume  was  in  operation. 
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OPERATIONAL  PROBLEMS 

A  number  of  operational  problems  and  decisions  significantly  affecting 
the  experiment  are  listed  here  in  an  effort  to  provide  insight  into  the  data 
reliability,  particularly  in  any  aspects  which  might  materially  affect  conclu- 
sions drawn  from  the  data. 

All  elevations  measured  in  the  flume  were  based  on  a  bench  mark  which 
was  assigned  an  elevation  of  10.000-ft.  A  flow  of  approximately  4-6  cfs  was 
maintained  when  bottom  elevations  were  taken  in  the  flume.  A  problem  which 
affected  depth  measurement  was  the  presence  of  standing  waves  in  the  test  flume. 
This  was  alleviated  by  the  placement  of  a  layer  of  6-in.  rubberized  horsehair 
mats  along  both  sides  of  the  test  flume.  Despite  the  installation  of  the  energy- 
absorbing  mats  slight  fluctuations  in  water  surface  elevation  still  existed. 
Minimum  fluctuations  of  0.002  ft  to  maximum  of  about  0.022  ft  in  double  ampli- 
tudes were  recorded.  The  wave  was  probably  generated  by  the  placing  of  the 
roughness  elements  aggravated  by  a  slight  slope  of  the  flume  in  the  north  to 
south  axis.  The  regular  spacing  of  the  elements  allowed  amplification  of  cortex 
strengths  and  resonance  of  these  vortices  interacted  with  vertical  flume  walls. 

One  of  the  difficulties  involved  in  the  experiment  was  the  accurate 
estimation  of  the  location  of  the  maximum  backwater.  Previous  studies  at 
Colorado  State  University  suggested  that  its  location  is  at  a  distance  of  one 
bridge  length  away  from  the  crossing.  The  search  for  its  location  was  made 
from  a  beam  (Figure  V-8)  having  an  adjustable  slope  and  a  moveable  point  gage. 
The  slope  was  set  equal  to  the  slope  of  the  flume's  bottom  surface.  Thus  by 
sliding  the  point  gage  along  the  beam  the  location  of  the  greatest  depth  was 
determined.  This  method  was  used  for  all  centered  and  eccentric  crossing 
configurations. 
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VI.   PRESENTATION  AND  ANALYSIS  OF  DATA 


INTRODUCTION 

This  chapter,  which  is  divided  into  five  major  subsections,  describes 
the  data  analysis  relating  to  the  objectives  of  the  experiment  listed  in 
Chapter  I.  Among  these  are  the  determination  of  the  effect  of  Froude  number, 
opening  ratio,  abutment  shape,  roughness  pattern,  skewness  and  eccentricity  on 
the  magnitude  and  location  of  maximum  backwater.  Due  to  the  unusual  nature  of 
the  roughness  elements  (i.e.,  protruding  through  the  free  surface),  an  important 
aspect  of  this  analysis  was  to  determine  their  relative  effectiveness  in  relation 
to  the  modeling  of  prototype  behavior. 

The  first  section  of  the  chapter,  "Dimensional  Analysis,"  presents  the 
theoretical  basis  for  the  data  analysis.  The  dimensionless  dependent  variables 
£*/b  and  h*/y  are  expressed  as  a  function  of  the  major  independent  variables 
varied  in  the  experiment.  This  analysis  is  extended  to  the  backwater  coeffi- 
cient, K*,  presented  in  the  third  section. 

The  "Analysis  of  Water  Surface  Elevation"  presented  in  the  second 
section  has  a  two-fold  purpose:  to  describe  the  sensitivity  of  the  most  basic 
parameter  recorded,  water  surface  elevation,  to  the  independent  variables  in  the 
experiment;  and  to  introduce  the  statistical  tool,  Analysis  of  Variance  (ANOVA), 
used  in  this  and  later  analyses  presented  in  the  chapter. 

The  ANOVA  statistically  estimates  the  presence  of  effects  in  the 
measured  data  caused  by  changes  in  the  independent  variables,  assuming  certain 
specific  conditions  of  the  data  are  met.  ANOVA  demonstrates  the  effect  caused 
by  changing  the  independent  variables  on  the  dependent  variable  being  measured 
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or  calculated,  i.e.,  the  effect  of  Froude  number  on  maximum  backwater  location. 
It  also  provides  an  estimate  of  the  reliability  of  the  conclusions  drawn  and  the 
ability  to  easily  handle  discontinuous  variables  such  as  abutment  shape  and 
roughness  patterns.  However,  due  to  time  exigencies  in  the  collection  of  data, 
only  a  small  portion  (12  percent)  of  the  data  can  be  analyzed  by  this  method. 
In  sections  four  and  five,  ANOVA  is  used  as  a  preliminary  screening  tool  to 
identify  trends  and  effects  in  the  analysis. 

The  fourth  section  presents  an  analysis  of  the  maximum  backwater 

location  and  the  fifth  concludes  with  an  analysis  of  the  height  of  maximum 

backwater.  Both  are  conducted  in  accordance  with  the  concepts  developed  in  the 
"Dimensional  Analysis"  section. 


DIMENSIONAL  ANALYSIS 

Any  given  physical  phenomenon  is  governed  by  a  set  of  independent 
quantities  or  characteristic  parameters.  The  laws  relating  physical  phenomena 
generally  are  expressed  in  the  form  of  mathematical  relations  among  the 
quantities  involved.  The  principles  of  dimensional  analysis  serve  as  a  vehicle 
for  deriving  the  functional  relationships  for  natural  phenomena.  This  analysis 
attempts  to  establish  these  relationships  required  to  describe  the  bridge 
backwater  phenomenon  among  the  pertinent  parameters  through  the  application  of 
dimensional  theory. 


DEFINITION  OF  VARIABLES 

Independent  Variables 

1 .  Geometry 
(a)  Channel 

Longitudinal  Slope  (S  ) 
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Shape  Factor  (?) 
Channel  Roughnes: 
Channel  Width  (B) 


Channel  Roughness  and  Distribution  (Cf) 


(b)    Bridge 

Abutment  Shape  (A) 
Opening  Size  (b) 
Skewness  (y) 
Eccentricity  (E) 
Bridge  Width  (w) 

.2.  Flow  Characteristics 
Discharge  (Q) 
Depth  (y) 

3.  Fluid  Characteristics 

Density  (p) 
Viscosity  (y) 

4.  Constant 

Gravitational  Constant  (g) 

De-pendent  Variables 

Maximum  Backwater  (h*) 

Maximum  Backwater  Location  U*) 


ANALYSIS 

From  the  variables  listed  above,  the  backwater  of  a  highway  stream 
crossing  can  generally  be  derived  from  the  following  characteristic  parameters 
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S  ,  g,  Cf,  B,  A,  b,  y>  E,  w,  Q,  y,  p,  y,  and  g. 

And  the  functional  relationship  betwetn  the  dependent  and  independent  variables 
may  be  expressed  by: 


h*,  i*     =  <h(S0,  5.  Cf,  B,  A,  b,  Yf  E,  w,  Q,  y,  p,  y,  g) 


(38) 


Replacing  Q  with  V  yields 


h*,  £*  =  <j>2(S0,  ?,  Cf,  B,  A,  b,  Y>  E,  w,  V,  y,  p,  y,  g) 


(39) 


Consider  the  flow  conditions  at  the  bridge  site  prior  to  the  construction 
idge.  A  given  dischan 
Equation  39  can  be  rewritten  as 


of  the  bridge.  A  given  discharge,  Q  ,  can  be  expressed  by  B,  y  and  V  .  Then 


h*,  £*  =  4>3(S0,  K,   Cf,  B,  A,  b,  Y,  E,  w,  VQ,  yQ,  p,  y,  g) 


(40) 


Selecting  p,  y  and  v  as  basic  quantities  (i.e.,  repeating  variables) 
we  obtain  the  following  dimensionless  equations: 


and 


b 


V  y  V  2 

Qrr^_/\t)    ,-W   oJo  o 

V  *'  f  y  '  '  B'  Y'  L'  B'  u  '  gy 

•'o  OJ 


V  y   V  2 

s   e  C   ^-  A  -  y  E  *-    -51° '  -°- 

o'  5'  f  y0'  '  B'  Y'  b'  B'  u   gyQ_ 


(41) 


(42) 


In  the  present  study,  the  bottom  slope  and  shape  of  the  flume  and  the 
bridge  width  are  fixed.  Therefore,  the  terms  S  ,  e,  and  =■  may  be  dropped  from 
Equations  41  and  42,  resulting  in  the  following  equations: 


h^ 

yn 


V  2  V  y 

o  oJo  ]d  B_  r     f  A 

gy.'  v)  '  B'  y  '  Lf  Y'  L'  rt 

_  0  0 


(43) 


and 
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b 


V   2    V   V      .     n 

o    cro  £  B_  r      r  fl 

Lgy0'  o  ■  b'  yo'  Lf  Y'  b'  A 


v„y. 


(44) 


The  lower  limit  of  the  Reynolds  Number  (i.e.,  °  °)  tested  was  approximately  in 
the  order  of  5,000,  which  is  in  the  turbulent  flow  region;  hence  the  effect  of  the 
Reynolds  Number  can  be  ignored.  Furthermore,  a  previous  study  at  Colorado  State 
University  has  shown  that  the  effect  of  the  Reynolds  Number  on  the  backwater 
is  negligible.  Equations  43  and  44  can  be  further  reduced  to: 


and 


Let 


h* 
yo 

"    *8 

"v  2 

0      b    B      r           fa 
gy   '  B'  y   '   Lf   Y'   L'  rt 

_     0               J  0 

_ 

b 

=    <j>9 

V  2 
o       b     B       r             F     . 

gy   '  B'  y   '  Lf  Y'   L»  AJ 

L     0                  0 

*Fr 

V 
_  _o_ 

b 

(45) 


(46) 


•'o 

hi 

yo 


and 


Then  Equations  45  and  46  become 


and 


5   =  *io(^r,  M,  x,  Cf,  0,  E,  A) 
e   =  *ll-(lFFr,  M,  t,  Cf,   ,  E,  A) 


(47) 
(48) 
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In  Equations  47  and  48,  nF  is  the  Froude  Number  of  the  flow,  M  is  the 
ratio  of  the  bridge  opening  to  the  valley  width,  t  is  the  width-depth  ratio  of 
the  stream,  C.  is  the  roughness  parameter  of  the  channel,  y   is  skewness,  E  is 
eccentricity,  and  A  is  the  shape  factor  of  the  abutments,  as  defined  previously. 
Equations  47  and  48  are  used  as  the  basis  for  the  data  analysis. 

Since  the  test  flume  provides  only  a  single  slope,  the  interrelation 
of  nF  and  t  cannot  be  examined  explicitly.  For  example,  it  is  impossible  to 
set  a  value  of  t  and  vary  values  of  IFF  without  changing  the  slope  of  the  flume 
and  the  flow  rate.  As  a  result,  a  change  in  flow  rate  will  cause  the  depth  and, 
hence,  the  value  of  B/y  to  change  accordingly,  since  B  is  a  constant. 

In  the  present  study,  the  value  of  B/y  ranges  from  13  to  82  in  the 
model  flume.  Based  upon  a  distortion  of  8.3  (i.e.,  100/12)  the  corresponding 
width-depth  ratio  for  the  prototype  is  110  to  680.  (Existing  field  data 
collected  at  approximately  100  bridge  sites  have  shown  that  the  values  of  the 
width-depth  ratio  vary  from  6  to  770.)  In  view  of  this  analysis,  the  effect  of 
the  width-depth  ratio  of  the  model  data  may  be  eliminated.  Equations  47  and  48 
can  then  be  rewritten  as 

6,e  =  <f,12  (*Fr>  M'  Cf  A'  Y'  E)  (49) 

Equation  49  has  six  independent  dimensionless  variables  useful  for 
the  development  of  a  functional  relationship  describing  6  and  e.   Of  course, 
y  and  E  are  only  important  for  skew  and  eccentric  model  openings.  The  elimina- 
tion of  these  two  leaves  four  independent  parameters.  The  development  of  the 
relationship  describing  6  and  e  are  developed  from  plotting  and  verified  by 
statistics.  Typically  in  this  type  of  analysis  heavy  reliance  is  placed  on 
Analysis  of  Variance  (AN0VA)  [17,18]  and  Regression  Analysis. 
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TWO-DIMENSIONAL  ANALYSIS  OF  WATER  SURFACE  ELEVATIONS 

This  section  presents  the  results  of  an  analysis  of  the  effects  of 
Froude  number  (F),  crossing  opening  ratio  (M)  and  abutment  shape  (A)  on  water 
surface  elevations  for  each  gage  position  used  in  collecting  data  over  the 
entire  flume  for  roughness  pattern  "C,"  i.e.,  main  channel  25  percent  and  flood 
plain  75  percent.  This  includes  36  runs  with  the  crossing  configuration 
centered  in  the  flume  at  the  150  ft  location.  Roughness  pattern  C  was  selected 
because  it  offers  a  complete  set  of  the  results  of  varying  Froude  number 
through  four  levels,  opening  ratio  through  three  levels  and  abutment  shape 
through  three  levels.  Complete  set  means  every  combination  of  three 
variables  was  tested  yielding  36  runs  (3x3x4).  Roughness  pattern  C  is  the  only 
set  of  centered  runs  meeting  this  definition.  Skew  and  eccentric  crossings 
are  also  complete  sets  but  with  a  different  set  of  variables,  i.e.,  they 
include  angle  of  skewness  and  eccentricity. 

The  purpose  of  this  analysis  is  to: 

1.  Determine  and/or  confirm  that  the  major  variables 
affected  backwater  in  a  measurable  way. 

2.  Determine  the  spatial  range  of  the  effects  of  these  variables 
in  the  flume. 

3.  Determine  where  normal  depth  was  re-established. 

4.  Introduce  the  concept  of  Analysis  of  Variance  (ANOVA) 
used  in  this  and  later  sections. 

This  analysis  demonstrates  that  the  opening  ratio,  Froude  number,  and 
abutment  shape  are  statistically  related  to  the  height  of  the  water  surface 
elevation,  each  independently  of  the  other.  This  is,  of  course,  by  no  means 
a  startling  conclusion.  However,  the  analysis  extends  this  over  the  entire 
flume  for  each  gage  where  data  were  taken.  This  permits  an  estimation  not 
only  of  where  opening  and  abutment  shape  affect  elevation  but  also  of  where 
they  do  not  affect  water  surface  elevation.  This  by  definition  represents  a 
return  to  normal  flow  conditions. 
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Implicitly  this  type  of  data  analysis  serves  as  a  very   good  indicator 
of  the  amount  and  magnitude  of  data  error  in  the  flume  along  with  an  indication 
of  the  sensitivity  of  the  dependent  variable,  water  surface  elevation,  to 
changes  in  the  independent  variables  (Froude  number,  opening  ratio,  and 
abutment  shape).  It  also  provides  an  indicator,  in  the  F  distribution,  as  to 
nonrandom  error  introduced  by  either  the  experimenter  or  model  design.  The  F 
distribution  is  a  statistic  used  in  testing  for  effects  and  specifically 
represents  a  ratio  of  mean  squares.  This  is  not  to  be  confused  with  Froude 
number  (F).  In  all  cases  the  meaning  should  be  clear  from  the  text. 

Since  the  ANOVA  is  used  in  following  sections,  we  will  digress  with  a 
brief  nontechnical  description  of  the  method,  its  advantages  and  disadvantages 
and  how  it  was  employed  here.  For  a  theoretical  presentation  the  reader  is 
referred  to  Reference  17.  The  purpose  of  the  following  section  is  merely  to 
acquaint  the  reader  with  the  principles  of  ANOVA. 


ANALYSIS  OF  VARIANCE 

An  analysis  of  variance  is  a  statistical  analytical  method  which  can 
be  used  on  certain  kinds  of  experimental  data  (i.e.,  data  meeting  certain  well 
defined  criteria)  in  order  to  reveal  the  importance  of  effects  of  individual 
independent  variables  on  some  dependent  variable.  In  this  problem  is  answered 
the  question  "Does  the  data  demonstrate  at  a  given  confidence  level  whether 
Froude  number,  opening  ratio  (M),  or  abutment  shape  have  a  direct  effect  on 
water  surface  elevation  or  does  any  combination  of  these  variables  have  an 
effect  on  water  surface  elevation?" 

Since  Froude  number  is  directly  linked  to  flow,  an  obvious  expected 
result  is  that  the  differing  Froude  numbers  have  an  effect  on  surface  elevation. 
It  is  also  expected  that  differing  M  ratios  should  show  significant  effects  at 
least  upstream  of  the  constriction.  It  is  not  so  obvious  whether  the  difference 
in  abutment  shapes,  i.e.,  wingwall  and  spill  through,  is  great  enough  to  cause 
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a  statistically  significant  effect  on  water  surface  elevation.  It  is  also  not 
obvious  what  the  spatial  distribution  of  these  effects  is  in  the  flume.  For 
example,  the  ANOVA  can  demonstrate  how  far  upstream  from  the  flow  constriction 
the  opening  ratio  affects  depth.  If  the  opening  ratio  affected  depth  over  an 
area,  that  area  would  be  classified  as  a  Zone  2  flow  region  (Figure  II-3). 

The  same  can  be  demonstrated  for  the  flume  below  the  constriction,  i.e., 
"Does  Zone  3  extend  the  length  of  the  flume  or  is  Zone  4  (normal  flow  and  de'pth) 
attained?"  The  answer  to  this  question  determines  the  location  where  normal 
depth  is  attained.  Throughout  this  report  any  question  answered  by  ANOVA  is 
assumed  to  be  prefaced  with  "According  to  the  statistical  model  used  does  the 
data  collected  demonstrate...?" 

The  ANOVA  used  in  this  section  and  generally  throughout  the  data 
analysis  is  a  three-variable,  factorial  design  with  one  replication.  "Three- 
variable"  refers  to  the  independent  variables,  i.e.,  Froude  number,  opening 
ratio,  and  abutment  shape.  Roughness  element  pattern  was  not  included  as  a 
variable  in  the  analysis  due  to  incompleteness  of  the  data  matrix.  Each  variable 
is  divided  into  levels.   In  the  analysis  there  are  four  levels  for  Froude  number, 
each  representing  one  of  the  four  flows  used,  three  levels  for  opening  ratio, 
corresponding  to  3,  5  and  8ft  openings  and  three  levels  representing  abutment 
shape  (wingwall,  spill  through  and  spur  dike).  For  ease  of  discussion  the  two 
abutment  shapes,  wingwall,  spill  through,  and  spur  dikes  are  collectively  referred 
to  as  abutment  shapes.  Each  represents  a  specified  level  of  the  variable  "A". 
Factorial  design  refers  to  the  combinations  of  the  levels,  in  this  case  the 
(4x3x3)  36  flume  runs.  This  means  there  are  36  water  surface  elevations  for 
each  point.  Ideally,  the  ability  to  discern  effects  is  greatly  increased  by 
replications;  however,  due  to  the  expense  of  the  data  this  was  not  done. 

An  important  difference  in  independent  variables  is  whether  they  are 
fixed  or  random.  Fixed  variables  are  variables  whose  levels  are  arbitrarily 
chosen  by  the  experimenter.  In  this  example,  abutment  shape  and  opening  ratio 
are  fixed  variables.  As  originally  planned,  flow  would  also  have  fallen  into 
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this  category  but  due  to  the  difficulty  of  generating  the  same  exact  flow,  as 
previously  described,  flow  was  classified  as  a  random  variable.  As  will  be  seen 

when  the  model  is  discussed,  this  is  generally  a  conservative  bias.  An  ANOVA 

which  uses  random  and  fixed  variables  is  referred  to  as  a  Mixed  Model. 

The  statistical  model  that  tests  effects  and  interactions  for  the  mix 
of  random  and  fixed  variables  is  developed  using  standard  procedures  described 
in  most  advanced  texts  [20].  A  list  of  the  variables  and  the  test  for  each 
is  shown  in  Table  VI-1 .  The  components  of  variance  represent  the  theoretically 
possible  variance  attributable  to  each  source  of  variance.  This  is  shown  under 
the  heading  "Components  of  Variance"  by  the  random  error  component,  sigma,  and 
a  second  and/or  third  component  attributable  to  the  effect  being  tested. 

The  effects  are  tested  by  calculating  the  various  ratios  shown  under 
the  column  "Test."  The  parameter  MS  represents  a  mean  square,  which  in  turn  is 
calculated  by  dividing  the  sums  of  squares  by  the  degrees  of  freedom.  Each 
mean  square  theoretically  can  contain  the  components  of  variance  attributable 
to  it  by  the  design,  as  shown  in  Table  VI-1.  Sums  of  squares  and  mean  squares 
are  the  most  mechanistic  part  of  an  analysis  of  variance.  The  calculating 
formulas  for  each  sum  of  squares  representing  each  component  can  also  be  found 
in  an  advanced  text  describing  an  ANOVA  (see  Myers,  Ref.  17).  The  calculating 
procedure  is  exactly  the  same  regardless  of  the  classification  of  variables  as 
fixed  or  random.  Note  that  in  the  development  of  the  testing  ratios  each  major 
variable  is  isolated  from  the  effect  of  the  others.  In  essence,  the  effect  of 
opening  ratio  on  surface  elevation  can  be  tested  with  the  effects  of  Froude 
number  and  abutment  shape  removed. 

For  example,  the  mean  square  for  Froude  number,  MSp,  contains  a 
random  error,  sigma,  and  possibly  a  variance  component  due  to  the  differences 
in  the  effects  of  Froude  number  on  the  dependent  variable,  water  surface 
elevation.  The  determination  of  whether  there  is  an  effect  of  Froude  number  is 
demonstrated  by  the  ratios  shown  under  the  heading  "Test"  in  Table  VI-1.  If 
the  conditions  of  normality,  independence  and  homogeneity  hold,  the  ratio  of 
two  mean  squares  can  be  shown  to  have  an  F  distribution.  Consequently,  the 
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TABLE  VI-1 
ANALYSIS  OF  VARIANCE 


Source  of  Variance 


Type 


df    Components  of  Variance 


Test 


Froude  (F) 

Random 

3 

Abutment  Shape  (A) 

Fixed 

2 

Opening  Ratio  (M) 

Fixed 

2 

FxA  Interaction 

Mixed 

6 

FxM 

Mixed 

6 

MxA 

Mixed 

4 

FxAxM 

Mixed 

12 

a2  +  BMne2, 


a2  +  QBne2A  +  Bne2FA 

a2  +  QMne2M  +  Mne2FM 

a2  +  BNe2 

a2  +  Mne2 

a2  +  Qne2  A  +  ne2 


FA 


QB 

HA   n9  FAM 


a2  +  ne2 


FAM 


MSF/MSFA 

MSA/MSFA 
MSM/MSFM 

No  Test 

MSFM/MSFfl 

MSMA/MSFAM 
MSFAH/MSFfl 


Where: 


x 

df 

sigma 

B 

M 

Q 

e 

MS 

n 


=  interaction, 

=  degrees  of  freedom, 

=  random  error  component  of  variance, 

=  number  of  levels  of  openings, 

=  number  of  levels  of  abutment  shapes, 

=  number  of  levels  of  Froude  number, 

=  variance  (either  fixed  or  random), 

=  mean  square, 

=  number  of  replicates 
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probability  of  obtaining  various  F  values  for  a  given  confidence  level  and  thus 
the  presence  of  an  effect  can  be  determined.  For  example,  if  the  effect  of 
Froude  number  is  negligible  on  water  surface  elevation  (which  it  is  not),  the 
hypothesis  that  BMno^O  would  be  true.  The  F  test  would  be: 

F(.05)     o2  +  BMn  o2 

(3,6)  "  1-  (51) 

a2  +   BMn  Q2^ 

A  prior  assumption  in  this  design  is  that  BMne§„=0.  The  reasons  for  this 
assumption  are  described  in  a  later  part  of  this  section.  Rewriting  Equation  51, 

c(.05)  _  a2  +  BMn  02  ,„>> 

h(3,6)  = [bl) 

a2 

If  the  hypothesis  is  true,  i.e.,  <C=0,  Equation  52  reduces  to 

^  =  i 

a2 

Of  course,  in  reality  the  above  ratio  is  never  exactly  1,  whether  the  hypothesis 
is  true  or  false.  However,  if  the  hypothesis  is  false  and  BMnej^O,  the  ratio 
is  greater  than  1.  How  much  greater  than  1  determines  the  level  of  confidence 
that  can  be  placed  in  the  result.  The  95  percent  level  of  confidence  for  the 
above  F  test  is  4.76.  For  the  hypothesis  that  differing  Froude  numbers  do  not 
have  an  effect  on  water  surface  elevation  to  be  rejected  at  the  95  percent 
confidence  level,  the  ratio  of  the  two  calculated  mean  squares  must  be  less 
than  4.76.  The  95  percent  confidence  level  should  be  interpreted  as  there  is 
a  95  percent  chance  that  if  the  ratio  equals  4.76  that  the  hypothesis  is  false 

(H  :BMne9=0)  and  Froude  number  does  have  an  effect  on  water  surface  elevation. 

v  o     z 

A  similar  process  is  used  to  test  the  other  effects  as  shown  in  Table  VI-1 
under  the  heading  "Test." 
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The  term  o2   in  each  component  of  variance  requires  special  mention. 
This  term  ideally  contains  only  random  error  and  as  such  has  an  expected  mean 
of  0  and  variance  of  1 .  In  the  F  test  presented  above,  note  that  a2  appears 
in  the  denominator.  Intuitively  this  would  indicate  that  the  larger  the  random 
error  term  becomes,  the  more  difficult  it  is  for  a  hypothesis  to  test  false 
(i.e.,  an  effect  exists).  In  experimental  problems  the  estimate  of  o2  is 
probably  one  of  the  most  important  and  difficult  problems  that  must  be 
answered.  If  random  measuring  error  or  nonrandom  error  is  too  large  in  comparison 
to  the  effect  being  measured,  positive  results  cannot  be  drawn  from  the  data. 
The  efficacy  of  the  F  test  lies  in  the  fact  that  this  problem  is  recognized  in 
the  analysis. 

In  the  analysis  presented  in  this  chapter  the  ANOVA  is  often  used  as 
a  screening  mechanism  to  determine  whether  sufficient  reliability  can  be  placed 
in  the  data  to  justify  further  analysis  using  regression  techniques.  It  is 
recognized  that  regression  analysis  alone  can  provide  results  similar  to  those  of 
the  ANOVA  but  since  ANOVA  more  easily  handles  the  noncontiguous  type  of  variables 
required  for  this  experiment,  that  analysis  was  performed  initially. 


ANALYSIS  OF  DATA  FOR  ROUGHNESS  PATTERN  C 

The  following  conditions  were  tested  for  roughness  pattern  C  (75  per- 
cent in  flood  plain,  25  percent  in  main  channel): 

Froude  Number       Opening  Ratio  (M)        Abutment  Shape  (A) 


.107 

.387 

.082 

.242 

.072 

.145 

.059 

Wingwal 1 
Spill  through 
Spur  Dike 


The  analysis  of  variance  took  the  form  previously  shown  in  Table  VI-1 
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It  is  standard  procedure  to  use  wi thin-cell  variance  as  an  estimate  of 
random  error,  sigma,  if  replicates  of  the  experiment  are  possible.  If  this 
proves  impossible,  as  in  this  case,  either  a  two-  or  three-way  interaction  can 
be  used  as  a  conservative  estimate  of  sigma. 

The  data  base  for  this  analysis  consisted  of  36  runs,  each  consisting 
of  approximately  100  gage  points.  Thus  at  each  point  36  measurements  of  water 
surface  elevations  were  expected.  An  examination  of  the  data  base  revealed 
that  for  most  gage  points,  the  data  matrix  (4x3x3)  was  complete.  However, 
occasionally  between  one  and  three  points  were  missing.  The  missing  values  were 
estimated  since  for  the  most  part  there  was  only  one  datum  lacking  and  although 
it  is  possible  to  analyze  the  matrix  of  data  with  a  limited  amount  of  missing 
data,  substantial  changes  in  technique  must  be  made  or  the  analysis  will  be 
biased.  The  data  were  highly  linear  for  opening  size  and  flow,  so  that  a 
linear  regression  estimate  of  missing  data  would  not  seriously  bias  results  and 
the  error  introduced  would  be  minimal.  A  "front  end"  was  constructed  and  added 
to  the  program  which  performed  a  linear  regression  over  the  three  directions  of 
the  matrix  to  generate  three  estimates  of  a  missing  value.  The  correlation 
coefficient  (R2)  was  computed  for  n>3,  and  the  best  value  was  selected  from 
three  estimates.  Generally,  R2  (for  the  regression  on  flow)  was  very   close  to 
unity,  and  estimates  were  based  on  that  fitted  line. 

Figure  VI-1  presents  the  analytical  results  in  a  gage  map  of  the 
flume.  As  may  be  seen  from  the  first  gage  map,  the  water  surface  elevation  at 
any  point  in  the  flume  is  affected  by  flow  at  the  95  percent  confidence  level. 
Since  flow  is  a  significant  variable  everywhere,  even  at  the  99.9  percent  level, 
it  is  not  shown  on  the  gage  map. 

Opening  ratio  (M)  was  found  to  significantly  affect  surface  elevation 
everywhere  above  the  constriction,  demonstrating  that  the  backwater  clearly 
extended  back  beyond  the  upper  limit  of  the  flume.  Note  also  that  the  effects 
of  opening  ratio  essentially  are  limited  to  above  the  constriction.  This  is, 
of  course,  to  be  expected. 
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The  variable  abutment  shape  was  found  to  affect  stage  in  a  region 
bounded  for  practical  purposes  by  the  120-ft  line  (30  ft  above  the  constriction) 
although  the  demarcation  is  not  clear  cut  and  some  other  gages  above  this  line 
appeared  to  have  been  affected. 

The  only  interactive  effect  which  was  significant  for  a  large  number 
of  gages  was  AxM.  This  is  the  interactive  effect  on  water  surface  elevation  of 
abutment  shape  and  opening,  independent  of  both  of  their  already  noted  main 
effects.  Interestingly,  in  general  terms  these  appear  to  congregate  within  the 
constriction;  there  appears  to  be  a  region  on  the  center! ine  above  the 
constriction  where  the  interaction  is  not  significant;  and  finally,  above  the 
120-ft  line  and  on  the  north  side  of  the  flume,  all  gage  depths  appear  to  be 
affected. 

Except  for  a  few  anomalous  two-way  interactions,  no  effects  of 
opening  or  abutment  shape  are  observed  below  the  constriction. 

This  leads  to  the  conclusion  that  normal  depth  has  been  reestablished 
below  the  constriction  well  within  the  flume.  The  line  of  reasoning  follows: 
Assume  normal  depth  has  not  been  established.  If  this  were  the  case  the 
various  stages  below  the  constriction  would  show  an  effect  of  either  opening  or 
abutment  shape  since  these  are  the  principal  changing  variables  affecting  depth 
in  this  experiment  other  than  Froude  number.  Certainly  the  reattainment  of 
normal  depth  would  be  affected  by  opening  ratio  if  this  were  the  case.  However, 
the  analysis  demonstrates  that  below  the  constriction  only  Froude  number  has  an 
effect  on  elevation.  This  would  indicate  normal  depth  has  been  reattained. 

Establishing  that  normal  depth  has  been  attained  allows  the  method  to 
be  applied  for  establishing  flow  based  on  specific  gages  as  described  in  the 
Appendix.  Only  those  gages  on  the  centerline  where  normal  depth  was  shown  to 
occur  were  used  as  estimators  of  average  depth  for  a  particular  run. 
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It  should  also  be  noted  that  the  analysis  of  variance,  like  every 
other  probability-based  statistical  analysis  using  a  confidence  level,  provides 
a  finite  probability  of  not  only  making  a  correct  decision,  but  a  wrong 
decision  as  well.  In  the  case  of  a  95  percent  confidence  level,  it  should  be 
noted  that  there  is  a  5  percent  chance  that  the  significant  ratio  is  coincidental, 
that  is,  anomalous,  and  the  effect  thus  indicated  as  being  significant  actually 
is  not.  This  will  occasionally  cause  the  analyst  to  assume  that  an  effect  has 
been  manifested  when  it  has  not.  This  may  be  partially  responsible  for  the 
occasionally  apparently  anomalous  significant  interactive  effects. 

Another  possible  explanation,  and  one  which  accounts  for  gage 
locations  where  several  interactive  effects  are  anomalously  significant,  is  a 
lack  of  validity  of  one  or  more  of  the  assumptions,  especially  that  error  is 
independent  of  the  variables,  i.e.,  the  error  introduced  in  the  surface  elevation 
measurement  for  changes  in  flow  is  the  same  as  for  changes  in  abutment  shape. 
This  is  especially  true  under  the  present  circumstance  of  using  one  of  the 
interactive  effects  to  estimate  error.   If  error  was  indeed  small  for  the 
interaction  but  large  over  other  variables,  this  would  be  manifested  as  anomalous 
groups  of  significant  variables  and  interactions.  Such  a  situation  may  well  have 
given  rise  to  the  condition  at  gages  38  and  1  and  those  downstream  showing  a 
significant  effect  other  than  flow. 

The  strong  effect  of  the  abutment  variable  on  water  surface  was 
selected  for  special  duty,  in  view  of  the  very   close  geometric  similarity  of  two 
of  the  abutment  shapes  (spillthrough  and  wingwall)  and  the  very   much  different 
spur  dike.  It  was  hypothesized  that  the  strong  effect  was  caused  by  the 
disparity  between  spur  dike  and  the  other  two  shapes. 

There  are  two  straightforward  ways  to  test  this  hypothesis  using 
analysis  of  variance.  One,  commonly  called  "orthogonal  contrasts,"  allows 
the  analyst  to  extract  the  maximum  amount  of  information  from  the  data  as 
dictated  by  the  number  of  degrees  of  freedom.  The  other  method,  which  was 
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considered  more  practicable  under  the  circumstances,  was  to  remove  from  the 
statistical  analysis  that  part  of  the  data  base  relating  to  spur  dike  and  to 
proceed  to  determine  whether  abutment  shape  was  or  was  not  a  significant 
independent  variable  relative  to  surface  elevation  at  the  various  gages.  Since 
the  analysis  had  already  been  performed  with  spur  dike  data  included,  such  an 
analysis  would  provide  for  two  related  pieces  of  information.  First,  were  the 
strong  effects  of  the  abutment  variable  due  to  the  large  structural  differences 
between  the  spur  dike  and  the  other  two  abutment  shapes,  and  second,  if  so,  just 
how  much  impact  did  the  difference  between  spill  through  and  wingwall  have  on 
surface  readings. 

The  results  of  this  analysis  are  presented  in  Figure  VI-2  in  the  form 
of  a  second  gage  map  of  the  flume.  Comparing  the  results. of  this  analysis  with 
those  for  two  abutment  shapes  and  spur  dikes  combined,  certain  interesting 
conclusions  may  be  drawn.  First,  it  appears  that  the  spur  dike  did  in  fact 
cause  some  of  the  significant  effects  of  the  variable  abutment  on  surface 
elevation,  particularly  between  the  130-ft  line  and  the  constriction  and  off  the 
centerline.  Interestingly,  however,  the  gages  situated  on  the  centerline  and 
between  the  130-ft  line  and  the  constriction  appear  to  be  affected  by  the 
difference  between  spillthrough  and  wingwall.  Also  interesting  is  the  apparently 
substantial  reduction  of  the  number  of  gages  exhibiting  the  MxA  interactive 
effect,  particularly  in  and  near  the  constriction.  It  would  appear  that  the 
strong  interactive  effect  observed  in  the  analysis  described  earlier  was  due 
to  a  considerable  extent  to  the  presence  of  spur  dike  data.  Thus,  interpreting 
this  observation,  we  may  conclude  that  while  spur  dike  alters  the  impact  of 
opening  size  on  depth  the  direct  effect  of  either  abutment  shapes  or  spur  dike 
is  strong  for  even  minor  differences  in  geometry. 

These  results  must  be  interpreted  carefully,  however,  Because  of  the 
approach  taken  for  investigating  the  impact  of  spur  dike,  the  total  degrees 
of  freedom  and  the  degress  of  freedom  for  the  abutment  variable  and  all  abutment- 
incorporating  interactions  are  reduced.  This  reduces  the  sensitivity  of  the 
analysis  to  the  presence  of  effects.  This  means  that  we  may  accept  without 
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reservation  those  gages  which  continue  to  show  some  effect,  but  have  reservations 
with  those  observations  where  effects  have  been  minimized  or  are  not  observed.  In 
terms  of  the  results  of  these  two  analyses,  it  is  reasonable  to  accept  the  premise 
that  spur  dike  was  not  the  sole  reasons  for  the  impact  of  the  abutment  variable 
on  surface  readings  near  the  constriction.  However,  the  observation  that  the 
interactive  effect  (AxM)  is  reduced  may  be  due  in  part  to  a  reduced  data  base. 


The  conclusions  from  an  analysis  of  roughness  pattern  C  (FP-75,  MC-25) 


are: 


1.  Opening  ratio  (M)  affected  surface  elevation  the  entire 
length  of  the  flume  above  the  constriction.  No  effects 
were  noted  more  than  5  ft  below  the  constriction. 

2.  Normal  depth  was  reattained  within  5  ft  below  the 
constriction  for  the  flows  and  openings  studied. 
Individual  effects  of  flow  with  a  given  opening 
size  in  regard  to  the  attainment  of  normal  depth 
were  not  investigated. 

3.  Water  surface  measurement  was  sensitive  enough  to 
detect  differences  in  wingwall  and  spill  through 
abutment  shapes  (A). 

4.  Froude  number  (F)  was  a  significant  factor  in 
water  surface  elevation  measurement. 

5.  Interactive  effects  between  abutment  shape  and 
opening  ratio,  though  significant  in  some  cases, 
are  generally  inconclusive  due  to  their  scatter. 
This  interaction  might  be  interpreted  as  an 
effect  of  turbulence  due  to  abutment  shape  and 
constriction. 

6.  .  The  use  of  spur  dikes  was  responsible  for  a  large 

portion  of  the  water  surface  elevation  effects 
attributed  to  the  variable  abutment  shapes. 

In  general,  the  ANOVA  analysis  demonstrates  that  the  primary  parameter 
recorded,  water  surface  elevation,  was  sufficiently  sensitive  as  measured  to 
demonstrate  the  effects  of  each  of  the  independent  variables,  Froude  number, 
opening  ratio,  and  abutment  shape  in  the  data  collected. 
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ANALYSIS  OF  TOTAL  BACKWATER  COEFFICIENT,  K* 


The  backwater  equation  as  developed  by  Liu  and  Bradley  [1],  solved 
for  K*  is  given  as: 


h^ 


K* 


'ri2 


a  V2/2g    a2 
2  n2'  3 


(W-M 


(53) 


Cross-sectional  area  of  the  opening  at 
normal  flow 


V   =  Velocity  in  the  opening  at  normal  flow 


n2 


In  terms  of  application  where  the  upstream  and  downstream  areas,  A2 
and  A^,  vary  little,  the  calculation  of  K*  is  essentially  a  function  of  the 
velocity  head,  V22/2g.  In  the  flume  experiments  of  the  present  study,  Aj  and 
A4  are  the  same  since  the  flume  was  of  a  regular  shape  with  both  Ax  and  A4 
reflecting  normal  depth.  Alpha  2  (a2),  which  gives  a  weighted  average  of  the 
kinetic  energy  distribution  across  the  opening,  was  essentially  equal  to  one  for 
the  flume  runs  reported.  However,  Equation  53  was  used  for  the  calculation  of 
all  backwater  coefficients  presented  in  this  section.1 

Previous  studies  [1]  have  established  the  experimental  relationship 
that  K*  varies  primarily  with 


1.  Stream  constriction  ratio  (M), 

2.  Type  and  shape  of  bridge  abutment  (A), 

3.  Number,  size  and  shape  of  piers  in  the  constriction  (P) 

4.  Asymmetric  position  of  the  crossing  in  relation  to  the 
flow  (E), 

5.  Angle  of  skewness  (y),  and 

6.  Froude  number  (F) . 


'Data  for  the  analysis  described  in  this  section  are  in  the  Appendix  under 
the  heading  "KB"  in  the  Consolidated  Data  Table. 
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K*   =  f(M,  A,  P,  E,  Y,  F,  Cf) 

In  the  flume  experiments  the  effect  of  piers  was  not  a  variable  and  the  effect 
of  skewness  and  eccentricity  was  measured  using  only  one  abutment  shape  and 
roughness.  These  conditions  allow  the  estimation  of  K*  for  centered  models: 

K*   =  fT(M,  A,  F,  Cf)  and 

and  for  skewed  and  eccentric  models: 

K*   =  f2(M,  F,  E,  Y). 

It  is  important  to  note  that  these  independent  variables  are  not 
theoretically  related  to  K*  with  the  exception  of  Froude  number,  F.  The 
backwater  coefficient  can  thus  be  rewritten  as: 


2h* 
n2  n 


ignoring  a2,  where  F  „  and  h  are  theoretical  values  related  to  normal  flow, 
a     a  z>       n2     n 

In  most  cases  the  degree  of  interdependence  between  independent  variables  does 
negate  the  usefulness  of  each  independent  variable  in  predicting  the  dependent 
variable,  K*. 

BACKWATER  COEFFICIENT    (K*)   FOR  CENTERED  CROSSING 

Analysis  of  Variance    (AN OVA) 

The  purpose  of  the  analysis  of  variance  as  used  in  this  and  following 
sections  is  to  serve  as  a  preliminary  screening  device  using  a  small  portion 
of  the  data  base  (approximately  12  percent).  Since  the  data  base  encompasses 
a  reasonable  distribution  of  the  abutment  shapes,  openings,  and  flows,  results 
derived  from  this  portion  of  the  analysis  can  be  considered  to  hold  for  the 
entire  data  base  within  limits.  Reasonably  strong  positive  or  negative 
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results  can  be  extrapolated  with  a  degree  of  certainty  while  marginal  results 
should,  of  course,  be  viewed  with  caution. 

In  any  case,  the  analysis  of  variance  will  complement  the  graphical 
exposition,  and  these  in  turn  will  complement  the  regression  analysis. 

The  K*  was  analyzed  by  ANOVA  using  the  same  data  as  used  in  the  section 
on  "Analysis  of  Water  Surface  Elevations."  The  dependent  variable,  water  surface 
elevation,  was  replaced  with  K*.  This  provided  (4x3x3)  36  estimates  for  K*. 
Each  estimate  represented  a  particular  combination  of  Froude  number  (calculated 
in  the  flood  plain),  opening  ratio  and  abutment  shape. 

Table  V 1-2  presents  the  output  of  the  ANOVA  and  Table  VI-3  the  results 
of  testing  for  effects.  The  abutment  shape  variable  is  the  only  observed  sign- 
ificant effect  at  the  95  percent  confidence  level  (Table  VI-3).  This  means  that 
the  only  variable  whose  various  categories  measurably  affected  K*  at  the  stated 
level  was  abutment  shapes.  Removing  spur  dikes  from  the  data  base  eliminated 
this  effect.  This  indicates  the  effect  on  K*  was  significant  for  spur  dikes  but 
not  for  wingwall  or  spil 1  through.  The  testing  also  indicates  that  the  effect  of 
the  various  categories  of  M  or  F  on  K*  was  not  statistically  significant.  For 
example,  the  effect  on  K*  of  differing  opening  ratios  was  not  observed  in  the 
data  at  the  95  percent  significance  level.  The  effect  of  Froude  number  on  K* 
was  not  statistically  significant  at  the  stated  level.  This  does  not  deny  the 
existence  of  a  relationship  but  rather  states  that  the  data  analyzed  do  not 
confirm  its  existence  at  the  given  level  of  confidence. 

In  order  to  further  test  for  significance,  the  conservative  estimate 
of  the  error  mean  square, a2,, was  relaxed.  By  observing  the  ratio  of  mean  squares 
in  Table  VI-3,  it  can  be  noted  that  the  three-way  interaction  mean  square  is 
considerably  less  than  the  interaction  mean  square  FA.  FA  was  used  as  an 
estimate  of  a2.  However,  after  the  three-way  mean  square  was  substituted  and 
the  proper  test  ratios  calculated,  no  significant  differences  were  noted  other 
than  those  previously  indicated. 
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TABLE  VI-2 
ANOVA  OF  K*  WITH  AND  WITHOUT  SPUR  DIKES 


WITH 

SPUR  DIKES 

WITHOUT  SPUR  DIKES 

Level 

s  of  Variabl 

es 

Levels 

;  of  Variables 

F 
M 

A 

4 
3 
3 

F 
M 
A 

4 
3 
2 

Grand 

Mean   1 .70 

Grand 

Mean   1 .55 

Source  of 
Variance 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Froude  (F) 

0.39 

3 

0.13 

0.24 

3 

0.08 

Opening  Ratio 
(M) 

1.13 

2 

0.56 

0.44 

2 

0.22 

F  X  M 

0.77 

6 

0.12 

0.36 

6 

0.06 

Abutment  (A) 

1.68 

2 

0.84 

0.13 

1 

0.13 

F  X  A 

0.83 

6 

0.13 

0.75 

3 

0.25 

M  X  A 

0.45 

4 

0.11 

0.25 

2 

0.12 

F  X  M  X  A 

0.59 

12 

0.04 

0.38 

6 

0.06 

TOTAL 


5.i 


35 


2.59 


23 
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The  conclusion  drawn  from  this  analysis  is  that  the  opening  ratio, 
abutment  shape  and  Froude  number  do  not  show  a  strong  effect,  if  any, 
on  K*  using  data  collected  for  36  runs  with  roughness  configuration  C.  It  does 
not  state  there  is  no  relationship,  but  rather  these  data  do  not  demonstrate  a 
relationship  at  a  reasonable  level  of  significance.  The  error  is  considered 
to  lie  in  the  large  data  scatter  as  shown  in  the  following  graphical  analysis. 

Graphical  Analysis 

A  graphical  analysis  of  all  the  backwater  coefficient  data  divided  by 
abutment  shape,  spur  dike  and  roughness  configuration  is  presented  in  Figures 
VI-3  to  VI-16. 

For  comparative  purposes,  the  plot  of  K*  vs  M  for  Type  I,  subcritical 
flow,  taken  from  Hydraulic  Design  Series  Bulletin  No.  1,  is  superimposed  on 
the  flume  data.  Also  included  are  previous  curves  developed  by  Colorado  State 
University  flume  studies.  Plotting  position  M,  which  is  the  ratio  of  the 
length  of  the  opening  to  the  width  flood  plain,  is  the  same  for  both  sets  of 
data.  For  convenience,  the  five  patterns  and  their  roughness  densities  are 
repeated  below: 

Density 


Pattern 

Flood  PI 

ain 

Ma 

in  Channel 

P 

100% 

100% 

A 

100 

50 

B 

75 

75 

C 

75 

25 

D 

50 

50 

Numbers  shown  beside  points  on  Figures  VI-4to  VI-16  are  theoretical 
Froude  numbers  calculated  in  the  opening  at  normal  depth,  taken  from  the 
Consolidated  Data  Tables  under  the  heading  "FROUDE."  Note  these  are  not  the 


108 


if) 

*> 

_ 

^ 

d 

1 — 

c 

K 

2 

2 
Ld 

LU 

O 

_l 

3 
00 

LU 

(/) 

O 

LU 

ro 

CJ 

a: 

i 

CO 
XI 

cc 

LU 

> 

II 

2 

LU 
< 

o 

C/> 

LU 

or 

£ 
^ 

LU 
Q 

3 

o 

O 

Li_ 

< 

_l 
Z> 

00 

< 

2 

> 

Z> 

X 

^ 

H— 

X 

a: 

r 

o 


q 
to 


T 
Q 


o 


A 


109 


<  H 
5  Q 


-* 


OQ 
\ 

JO 


CO 

cr 

*> 

LU 

^£ 

i- 
i- 

h- 

< 

2 

Q_ 

Ld 

CO 

O 

CO 

sf- 

L_ 

LU 

i 

Li_ 

-z. 

LU 

X 

> 

o 
o 

o 
o 

LU 

cr 

cr 

Ld 

cr 

3 

H- 

X 

Lu 

^ 

e> 

£ 

ID 

^ 

O 

CJ 

cr 

< 

X 

CD 

b 

^ 

_l 

3 

CL 

^ 

CO 

X 

< 

cr 

o 


o 

to 


o 


_j_H L0 


110 


e> 


5P 

Q. 
Ul 

±.  a 

Ul  _i 
a.  < 
O  2 

5§ 
cz 

(C 

UJ 

ffl 

id 


ui 

Q 

Z> 

o 

q: 

u 

i 

x 

X 

X 


Ul 


CO 
Uj 

Q 

3 


o 

•J1 


CD 


QO 


to 

or 

* 

UJ 

h- 

h- 

fer 

^ 

Q_ 

UJ 

— 

cn 

U 

CO 

U_ 

LU 

U. 

p* 

LU 

T 

o 

CO 

u 

z> 

or 

o 
or 

h- 

T 

<. 

CO 

^ 

z> 

*: 

o 

o 

or: 

< 

X 

CD 

h- 

_i 

^ 

_i 

Z) 

Q 

;> 

if) 

X 

< 

or: 

in 
> 

LU 

or 

3 
CD 


111 


<  H 
Q. 
Ul 


O 


^  a 

a.  < 

S  o 
coz 

^  ° 

Ul  H 

oo 

u.  3 

i 

X 
X 
X 


o 


O 
to 


o 


o 


CD 
\ 

-Q 
II 

2 


in 

~Z. 
LU 

O 

u_ 

LU 

o 
u 

CT 
LU 
h- 

I 

u 

< 

CQ 


X 

< 


o 

z 
or 

LU 

l- 
co 

CO 

UJ 

X 

o 
o 

X 

o 

o 
or 

x 
H 


0_ 
CO 

or 
o 
u_ 


CO 

I 


or 

CD 


J 


112 


t  x 

<  K- 

is 

Q.  < 

O  2 

5§ 

c/>z 

<*  Q 


oo 
tr  z 

i 

X 
X 
X 


_j 

3 

e> 
o 

_t 
o 

UJ 

o 


p 


~\ — h 


o 
to 


p 

CM 


"T~ 

O 


CD 
\ 

-Q 


Q 


c/> 

cr 

*> 

Ld 

^ 

h- 

r- 

£ 

2 

cr 

LU 

o 

CO 

CO 

Ll_ 
U_ 

LU 
O 

u 

LU 

X 

1 

O 

3 

> 

cr 

O 

LU 

LU 

cr 

cr 

h- 

3 

< 

X 

CD 

o 
< 

CD 

(3 

3 
O 

cr 

X 

b 

li_ 

2 

_j 

3 
2 

CO 

X 

< 

cr 

o 


113 


o 
I) 

♦ 


<  t- 

IS 

Q.  < 

O  2 


Ul 


UJ 


i§ 

z  q: 

uj  H 

o  o 

u.  3 

X 
X 
X 


03 


UJ 


k> 

2 

^ 

cr 

H 
Ld 

UJ 

O 

LL 

U_ 

cn 

Li_ 
UJ 

CD 

hi 

(J 
U 

-z. 

I 

rr 

o 

UJ 

ID 

h- 

o 

< 

tr 

£ 

V 

UJ 

CJ 

^ 

< 

n 

CD 

q: 

^ 

Z> 

ID 

CL 

2 

CO 

X 

< 

or 

CO 

> 

UJ 

or 

Z> 
CD 

U. 


O 


o 
to 


o 


o 


114 


h 

X 

< 

1- 

Q. 
Ul 
Q 

7 

III 

_l 

Q. 

< 

O  5 

ir 

^ 

o 

;* 

</) 

£C 

n 

UJ 

CO 

2 
=> 
2 

UJ 

1- 
o 

q: 

hi 

h 

Q 

co 

m 
oo 

U.  =} 

x 

X 
X 


CD 


CD 


CO 

2 
or 

H 

1- 

~P 

h- 

III 

< 

n 

O 

U_ 

CO 

u_ 

CO 

LU 

UJ 

o 

7* 

u 

X 

rr 

(D 

ii  i 

_) 

h- 

O 

< 

tr 

£ 

Y 

LU 

U 

*: 

< 

Q 

m 

cr 

^ 

3 

Z> 

a. 

2 

co 

X 

< 

cr 

0> 


LU 

cr 

Z> 

o 

U_ 


-*x 


115 


5  Q 


en 
ce 
ui 

CD 

2 
ZD 
Z 

u 

Q 

go 

<£  Z 
U.  3 
I 

X 
X 
X 


-<*. 


CO 

\ 


O 


z 

*: 

or 

— ' 

Lu 

h- 

1- 

LU 

1- 
< 

O 

Q_ 

U_ 

to 

U_ 

CO 

LU 

LU 

O 

z 

U 

X 

Ld 

g 

CD 

z> 
o 
cr 

B£ 

UJ 

^. 

^ 

o 

< 

Q 

CD 

or 

ID 

2 

Q_ 

Z) 

CO 

^ 

or 

X 

< 

o 

Ll 

o 

> 

LU 

cr 

O 


ay 

o  u. 


o 


o 
to 


Q 

CM 


o 


-uo 


*X 


116 


?  o 

a.  < 
O  2 

a: 

UJ 

m 

5 

2 


a 
ui 

i- 
o 

3 
IK 
UJ  H 

=>  s 

o  o 

CE  Z 
U.  3 
• 

X 
X 
X 


T- 

o 


-r- 

o 

ro 


o 


-I j — 


ffl 


Q 


V) 

2 

* 

Cd 
UJ 

2 

1- 

5 

Ld 

CL 

O 

U_ 

CO 

U. 

CO 

UJ 

UJ 

O 

2 

U 

X 

O 

CC 

=> 

UJ 

o 

h- 

or 

< 

£ 

UJ 

:*: 

^ 

o 

< 

Q 

CD 

cr 

2 
3 
2 

Z> 
CL 
CO 

X 

(T 

< 

O 

> 

UJ 

a: 

Z> 

o 


T" 

o 


o 


117 


t  x 
<  I- 

is 

a.  < 
O  2 

S§ 

coz 
a: 

UJ 

CD 

5 


OO 

U.  3 
I 

X 
X 
X 


-<* 


m 
\ 

■i 

2 


in 

•z. 

v: 

ce. 

UJ 

h- 

h- 

^ 

H 

hi 

^ 

LI. 

o 

U_ 

CO 

U_ 

CO 

UJ 

UJ 

o 

CJ 

X 

cc 

UJ 

CD 
O 

< 

(T 

o 

< 
m 

_l 
_l 
< 

CD 

X 

< 


C\] 

> 

UJ 

a: 

Z> 
CD 


O 

U_ 


o 


o 
ro 


O 


"T" 

o 


118 


o 


o 


*  ... 

CM  U3 

C-J  — 


s^ 

a.  < 
O  2 

coz 

IK  Q 


-^; 


X 

< 


GO 


CO 

2 

a: 

^ 

UJ 

2 

UJ 

a: 

o 

CO 

Lt_ 

CO 

U- 

LU 

UJ 
O 

X 

o 

£   cr 

=> 

f     UJ 

o 

*   !5 

or 

£ 

_i 

^ 

_j 

!      CJ 

< 

< 

ca 

ro 

> 

LU 

a: 

=> 
O 

li- 


ar 
o 

Ll 


C> 


-T- 

O 
to 


Q 
CM 


I 

o 


-I— 


•o 


119 


o 


<  I- 
a. 

Id 

±  a 

a.  < 
O  2 

sS 

Sh 
o 

o: 

H 
co 
m 
o 

3 


UJ 

a 

O 
K 
U- 
1 

X 
X 
X 


o 


o 
to 


a 


o 


m 


X 
< 


o 


*> 

z 

^: 

tr 

UJ 

h- 

h- 

2 

H 

UJ 

< 

— 

Q_ 

O 

CO 
CO 

Ld 

UJ 

o 

2 

CJ 

X 

CD 

cr 

Z) 

LU 

h- 

< 

o 
a: 

H? 

_i 

^ 

_j 

o 

< 

< 

£ 

CD 

CD 

2; 

> 

UJ 

en 

Z> 
CD 


a: 
o 
u_ 


120 


DL   < 

O  5 


hj 


Ul 


1^ 

iij  I— 

oo 

CC  Z 

I 

X 
X 
X 


-<*. 


a: 

^ 

UJ 

LU 

cl 

) 
CO 

O 

u_ 

LU 

o 

CO 
CO 
UJ 

CJ 

X 

ii 

e> 

3 

cr 

Z> 

LU 

O 

h- 

or 

< 

j 

£ 

-j 

^ 

_i 

cj> 

< 

< 

£ 

CD 

o 

lO 

> 

UJ 

or 
S2 


x 

< 


or 
o 
u_ 


o 


-r- 
q 
to 


121 


CD 


<  H 
Q. 

UJ 
^  Q 

S^ 
Q.  < 

O  5 

sS 

coz 
<r  q 

&  UJ 

Si- 
^  o 


UJ  H 

o  o 
<r  z 
u.  ID 
I 

X 
X 
X 


•  • 


O  Ll 


T 
o 


m 


I- 

UJ 

o 
u_ 

Ll. 

UJ 
O 
O 

cr 

UJ 

I 

o 
< 

CD 


X 

< 


a: 

Ld 

t- 

< 

CO 
CO 

UJ 

z 

X 
CD 

z> 
o 
or 


< 

CD 


o 
u_ 


CD 

> 

UJ 

en 

Z> 
CD 
Ll. 


+  M 


122 


same  experimental  Froude  numbers  used  in  the  ANOVA.  The  lines  shown  on  each 
graph  are  from  top  to  bottom:  USGS  field  data,  and  two  lines  showing  Colorado 
State  University  flume  data.  These  are  shown  separately  on  Figure  VI-3.  Of  the 
K*  values  reported  in  the  Appendix,  only  those  less  than  0.5  and  greater  than 
5.0  are  now  shown.  All  runs  which  had  a  Froude  number  greater  than  1.0  were 
also  excluded. 

The  data  do  not  show  the  consistency  required  for  definitive 
conclusions  and  quantitative  assessments  of  major  variables.  As  indicated  by 
the  ANOVA  the  plotting  parameter (M=b/B)  as  shown  on  the  figures  does  not 
indicate  a  strong  relationship,  if  any,  to  K*  for  most  of  the  figures. 
Nevertheless,  certain  qualitative  assessments  can  be  made: 

1.  The  effect  of  varying  roughness  patterns  is 
indeterminate  in  regard  to  the  maximum  backwater 
coefficient.  This  might  be  alleviated  by  a 
subjective  sorting  of  the  plotted  data,  not  done 
here. 

2.  An  inverse  relationship  between  Froude  number  as 
calculated  in  the  opening  and  K*  appears  to  exist 

for  any  opening  ratio.  Note  here  that  due  to  the  resis- 
tance of  the  roughness  elements  the  higher  Froude 
numbers  correspond  to  lower  flows.  Effective  Manning 
roughness  in  the  flume  increased  with  flow. 

3.  As  expected,  scatter  tended  to  decrease  with  decreasing 
opening  ratio.  This  was  due  to  the  small  backwater 
produced  by  large  openings  and  therefore  the  difficulty 
in  accurately  measuring  this  parameter. 

4.  Finally,  as  also  indicated  by  the  ANOVA,  any  future 
analysis  of  K*  could  combine  wingwall  and  spillthrough 
data  into  one  set  of  data  due  to  the  similarity  of 
results  for  each  roughness  pattern. 
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BACKWATER  COEFFICIENT    (K*)   FOR  SKEWED  CROSSING 

Analysis  of  Variance 

All  skewed  crossing  runs  were  made  with  roughness  pattern  D  (flood 
plain  50  percent  and  main  channel  50  percent),  using  45°  wingwall  abutment 
faces.  Categorization  by  Froude  number  was  based  on  the  Froude  numbers 
calculated  using  the  entire  flume  at  the  given  flow  and  normal  depth.  Eccentric 
categories  were  divided  into  north,  central  or  south,  depending  on  the  opening 
orientation  to  the  north  or  south  sides  of  the  flums.  The  definition  of 
eccentricity  defined  earlier  in  this  report  is  not  used  here  as  it  does  not 
adequately  assess  the  direction  of  eccentricity  from  the  centerline.  Eccentricity 
as  used  in  this  part  of  the  analysis  is  defined  in  Figure  VI-17.  This  definition 
is  applied  in  order  to  take  into  account  the  expected  differences  in  the 
maximum  backwater  height  and  location  caused  by  energy  losses  due  to  eddies 
generated  by  the  opening  location.  Calculations  of  eccentricity  for  each  run 
are  contained  in  the  Appendix. 

Opening  ratios  were  categorized  according  to  that  portion  of  the 
opening  perpendicular  to  the  flow  for  both  sizes  of  opening.  Data"  were  further 
divided  into  15°  and  45°  levels.  Levels  for  each  variable  are: 

Variable 15£ 45f 

Froude  (F)  .077  .091  .129  .077  .091  .129 
Eccentricity  (E)  North  Central  South  North  Central  South 
Opening  Ratio  (M)  2.89     7.73  2.12     5.66 

The  angle  of  skewness  was  not  analyzed  because  the  opening  ratio 
categories  varied  with  the  angle  of  skew.  This  made  it  impossible  to  statisti- 
cally test  for  variances  in  K*  due  to  the  angle  of  skew.  An  examination  of 
Table  VI-4  shows  that  18  (3x3x2)  estimates  of  K*  are  needed  for  both  the  15° 
and  45°  skewed  crossings 
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NORTH   WALL 


SOUTH  WALL 
a)  NORTH  ECCENTRICITY 


NORTH    WALL 


SOUTH   WALL 
b)    CENTRAL    ECCENTRICITY 


NORTH   WALL 

SKEWNESS  \ 

^  FLOW 

SOUTH  WALL 
c)    SOUTH  ECCENTRICITY 

EXPERIMENTAL  ECCENTRICITY    DEFINITION 
FIGURE  VI-17 
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TABLE  VI-4 
ANOVA  OF  K*  FOR  15°  AND  45°  SKEW  CROSSINGS 


15°  SKEW  CROSSINGS 
Levels  of  Variables 


45°  SKEW  CROSSINGS 
Levels  of  Variabiles 


F 
E 
M 

3 
3 
2 

F 
E 

M 

3 
3 
2 

Grand 

Mean 

1.68 

Grand  Meai 

i  1 

.57 

Source  of  Sums  of 
Variance   Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Froude  (F)   1.93 

2 

0.96 

.93 

2 

0.46 

Excenctricity 
(E)       1.20 

2 

0.60 

1.48 

2 

0.74 

F  x  E  (FE)    .64 

4 

0.16 

0.93 

4 

0.23 

Opening  Ratio 
(M)       1.09 

1 

1.09 

0.29 

1 

0.29 

FxM  (FM)    .34 

2 

0.17 

0.32 

2 

0.16 

E  x  M  (EM)   1.59 

2 

0.79 

5.88 

2 

2.94 

F  x  E  x  M 

(FEM)      1.71 

4 

0.42 

0.43 

4 

0.10 

TOTAL 


8.54 


17 


10.30 


17 
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Table  VI-4  presents  the  sums  of  squares,  degrees  of  freedom  and 
mean  square  calculations  for  the  skew  runs.  Table  VI-5  presents  the  results  of 
testing  each  factor  at  the  95  percent  confidence  level.  No  significance  could  be 
attributed  to  any  of  the  main  effects  for  either  15°  or  45°  skew  runs.  The  one 
interactive  effect  of  eccentricity  and  opening  ratio  for  the  45°  skew  runs  is 
considered  an  anomaly  since  neither  of  their  main  effects  tested  as  significant. 
The  conclusion  drawn  from  the  ANOVA  is  that  the  data  show  no  statistically 
significant  effect  at  the  95  percent  confidence  level  which  may  be  attributable 
to  Froude  number  (as  calculated  in  the  entire  flume),  opening  ratio,  or 
eccentricity.  Effects  between  15°  and  45°  crossing  configuration  angles  were 
not  testable  because  of  insufficient  data. 

Graphical  Analysis 

Plots  of  K*  vs  M  are  given  in  Figures  VI-18  to  VI-20.  Numbers 
associated  with  each  point  refer  to  the  theoretical  Froude  number  calculated 
in  the  constriction  at  normal  depth.  Triangles  refer  to  a  south  eccentricity 
and  squares  to  a  north  eccentricity.  All  others  were  centered  in  the  flume  at 
the  angle  of  skew  shown. 

Graphically,  the  same  difficulties  regarding  measurement  error  are 
evidenced  in  the  skew  runs  as  were  found  in  the  centered  crossings,  perpendicular 
to  flow.  Measurement  error  decreases  with  increasing  backwater  as  noted  by  the 
clustering  of  points  at  smaller  M  ratios.  As  indicated  by  the  ANOVA,  the 
effects  of  opening  ratio  cannot  be  graphically  quantified  at  a  reasonably 
significant  level.  The  clustering  of  points  shows  that  the  Froude  value 
calculated  in  the  opening  perpendicular  to  flow  using  the  normal  unobstructed 
depth  has  an  inverse  relationship  to  the  backwater  coefficients.  However,  this 
relationship  has  not  been  quantified  by  any  further  regression  analysis. 
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TABLE  VI-5 
TESTING  OF  K*  FOR  SKEW  CROSSING 
AT  95  PERCENT  CONFIDENCE  LEVEL 


Test 

F  Statistic 

Siani 
15s 

Source  of  Variance 

Req.  for 
Signif . 

Calci 
15° 

jlated 
45° 

ficant 
45° 

Froude      (F) 

MSF/MSFEM 

6.94 

2.26 

4.24 

No 

No 

Eccentricity  (E) 

MSE/MSFE 

6.94 

3.72 

3.18 

No 

No 

FxE 

MSFE/MSFEH 

6.39 

0.38 

2.13 

No 

No 

Opening  Ratio  (M) 

MSM/MSFM 

18.5 

6.30 

1.83 

No 

No 

FxM 

MVMSFEM 

6.94 

0.40 

1.48 

No 

No 

ExM 

MSEM/MSFEM 

6.94 

1.86 

26.7 

No 

Yes 

FxExM 

No  Test 

-- 

-- 

— 

— 

-- 
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The  parameter  h*  used  in  the  calculation  of  K*  represents  the  height 
above  normal  water  surface  elevation  at  the  location  of  the  maximum  backwater. 
This  involves  three  sources  of  error:  estimating  the  location  for  measurement, 
the  water  surface  elevation  and  the  height  of  the  normal  water  surface  elevation, 
In  collecting  and  analyzing  the  data,  emphasis  has  been  placed  on  finding  the 
maximum  backwater  position  and  estimating  the  most  accurate  normal  water  surface 
elevation  possible.  These  errors  have  been  minimized  to  the  extent  possible. 
The  remaining  error  is  considered  attributable  to  the  small  calculated  values  of 
h*  due  to  the  combination  of  flows  and  opening  ratios. 
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ANALYSIS  OF  MAXIMUM  BACKWATER  LOCATION  U*/b) 

The  location  of  the  maximum  backwater,  l*,   was  defined  as  the  distance 
from  the  upstream  crossing  embankment  to  the  location  of  the  maximum  water 
surface  elevation  above  normal  surface  elevation.  The  location  of  this  point 
was  found  experimentally  using  the  apparatus  and  methodology  described  in 
Chapter  IV. 

ANALYSIS  OF  VARIANCE    (l*/b) 

As  has  been  done  previously,  an  analysis  of  variance  was  performed  on 
the  36  runs  comprising  roughness  pattern  C,  75  percent  flood  plain  and  25 
percent  main  channel.  The  purpose  of  this  test  was  again  to  provide  some 
statistical  screening  of  the  data  in  order  to  establish  a  level  of  significance 
for  the  important  variables.  The  dependent  variable  tested  was  i*/b.     An  ANOVA 
was  performed  with  and  without  the  spur  dike  data.  The  variables  tested  were 
Froude  number  as  calculated  in  the  flood  plain  (F),  opening  ratio  (M),  and 
abutment  shape  (A)  (Figure  V-3).  The  levels  of  each  variable  are  given  below: 

Froude  Number  (F) 
Opening  Ratio  (M) 
Abutment  Shape  (A) 

Output  of  the  ANOVA  and  significance  testing  of  the  major  variables  and 
interactions  is  given  in  Tables  VI-6  and  VI-7.  Note  the  estimate  of  the 
residual  mean  square,  sigma,  is  the  interaction  mean  square,  FA.  This  was  done 
in  order  to  give  a  slight  conservative  bias  over  the  three-way  interaction 
mean  square,  FMA.  No  differences  in  significance  results  from  the  choice  of 
either  as  an  estimate  of  sigma. 

Results  of  significance  testing  demonstrate  that  all  three  major  variables 
affect  the  dependent  parameter,  a*/b.     The  effect  noted  for  the  variable  abutment 
shape  appears  attributable  to  spur  dike  alone.  Conversely,  no  differing  effect 
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.059 

.072    .083 

.134 

.145 

.242    .387 

ingwall 

spill  through 

spur  dike 

TABLE  VI-6 
ANOVA  OF  l*/b   WITH  AND  WITHOUT  SPUR  DIKES 


WITH 

SPUR  DIKES 

WITHOUT 

SPUR  DIKES 

Leve' 

Is  of  Factors 

Levels  i 

of  Factors 

F 
M 

A 

4 
3 
3 

F 

M 
A 

4 
3 
2 

Gram 

i   Mean    1 . 

31 

Grand  Mean    1 .17 

Source  of 
Variance 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Froude  (F) 

1.64 

3 

0.54 

1.44 

3 

0.48 

Opening  Ratio 

(M)       9.56 

2 

4.78 

4.94 

2 

2.47 

F  x  M 

0.43 

6 

0.07 

0.43 

6 

0.07 

Abutment  (A) 

1.40 

2 

0.70 

0.01 

1 

0.01 

F  x  A 

0.20 

6 

0.03 

0,12 

3 

0.04 

M  x  A 

1.05 

4 

0.26 

0.27 

2 

0.13 

F  x  M  x  A 

0.12 

12 

0.01 

0.05 

6 

0.01 

TOTAL 


14.43 


35 


7.29 


23 
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on  £*/b  at  the  95  percent  confidence  level  is  noted  for  wingwall  and 
spillthrough  alone.  The  interactive  effect  of  opening  ratio  and  the  variable 
abutment  shape  may  be  due  to  their  interactive  effect  on  turbulence. 

GRAPHICAL  ANALYSIS    (l*/b) 

Figures  VI-21  to  VI-27  present  plots  of  all  the  i*/b   data  generated 
by  the  experiment.  Plotting  numbers  are  flood  plain  Froude  numbers.  The 
graphs  give  further  evidence  to  that  indicated  by  ANOVA  for  roughness  element 
patterns  C.  Opening  ratio  and  Froude  number  do  appear  functionally  related 
to  £*/b  and  the  addition  of  spur  dikes  to  the  embankment  faces  does  affect 
the  magnitude  of  the  £*/b  parameter.  These  facts  appear  to  be  significant 
for  all  the  data  taken,  though  except  for  roughness  pattern  C  they  have  not  been 
proven  statistically  significant.  Again  any  future  analysis  of  £*/b  could 
combine  wingwall  and  spillthrough  data  with  only  a  slight  increase  in  data 
scatter. 

Due  to  the  few  points  taken  for  any  set  of  combinations  and  the 
difficulty  in  pooling  the  data  onto  one  or  two  graphs,  a  regression  analysis 
is  presented  below  using  the  same  data  shown  in  the  figures.  This  difficulty 
arises  because  of  the  impossibility  of  keeping  Froude  numbers  in  the  experimental 
setup  constant  between  roughness  patterns  and  the  incomplete  testing  due  to  time 
exigencies. 

REGRESSION  ANALYSIS    (l*/b) 

The  notation  describing  the  regression,  taken  from  the  dimensional 
analysis,  is  given  below: 


rVb  =  f(M,  A,  F,  Cf)  (54) 


where 


M  =  Opening  ratio,  b/B 
A  =  Abutment  shape 
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F  =  Froude  number  (flood  plain) 
Cf  =  Roughness  element  pattern  (A  through  D) 

Regression  requires  a  general  form  of  the  equation  to  be  defined  before 
proceeding.  In  this  analysis  only  linear  and  log-linear  were  considered.  These 
were  chosen  since  the  data  appeared  to  follow  a  log-linear  to  linear  pattern 
as  shown  in  the  preceding  section.  Though  both  were  analyzed  only  the  log-linear 
is  presented  here  due  to  its  generally  superior  performance.  The  form  of  the 
log-linear  equation  is  given  below: 

^  -  b,  Hb2  Ab2  Fb3  Rb"  (55) 

Table  VI-  8  shows  multiple  correlation  coefficients  (R)  for  the  four 
independent  variables  in  the  experiment.  The  data  in  this  table  have  been 
log  transformed.  Note  that  abutment  shape  is  treated  as  a  categorical  variable. 
In  the  regressions  that  follow  this  means  that  abutment  shapes  (wingwall,  spill- 
through  and  spur  dike)  are  coded.  These  are  given  below: 

Abutment  Variable 


X, 

x5 

Wingwall 

1 

1 

Spill  through 

1 

10 

Spur  dike 

10 

1 

Roughness  patterns  are  also  coded  by  using  A  values  defined  In 
Chapter  IV.  These  are  repeated  here  for  reference. 
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TABLE  VI- 8 
CORRELATION  COEFFICIENTS  FOR  INDEPENDENT  VARIABLES 


Variable 
Number 


1  1.000    -0.049    -0.691  0.536  -0.281  -0.151 

2  1.000         -0.205  0.071  0.081  -0.553 

3  1.000         -0.177  0.048  0.010 

4  1.000  -0.495  -0.116 

5  1.000  -0.112 

6  1.000 


where  Variable  1  =  £*/b 

2  =  Froude  Number  (flood  plain)  (F) 

3  =  Opening  Ratio  (M) 

4  =  Abutment  Variable  (categorical) 

5  =  Abutment  Variable  (categorical) 

6  =  Roughness  Element  Pattern 
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Roughness  Element  Pattern        Lambda  (x) 


p 

0.0242 

A 

0.0198 

B 

0.0160 

C 

0.0144 

D 

0.0107 

The  effect  of  coding  is  to  identify  categories  into  which  the 
dependent  variable  can  be  classified  in  an  attempt  to  "explain"  variance. 
The  numerical  coding  in  no  way  affects  the  tests  of  significance  or  conclusions 
drawn.  However,  it  does  modify  coefficients  and/or  exponent  values. 

Table  V I -9  presents  results  of  four  log-linear  regressions  completed 
on  the  i*/b   data.  Regression  equation  1  is  the  best  least  squares  fit  to  the 
data  obtainable  from  the  data  presented  previously  for  i*/b.     Approximately 
77  percent  of  the  total  variance  is  "explained"  by  these  four  independent 
variables.  All  the  exponents  test  significant  at  the  95  percent  confidence 
level.  Figure  VI-28  is  a  plot  of  the  log  residuals  vs  log  of  i*/b.     Note  that 
no  unusual  distribution  of  the  residuals  is  obtained.  Approximately  one-half 
the  residuals  are  positive  and  one-half  negative. 

Analysis  2,  Table  VI -9,  is  the  functional  relationship  between  i*/b   on 
M  and  F  alone.  Fifty-two  percent  of  the  total  variance  is  accounted  for.  This 
would  indicate  the  abutment  shape  variable  and  roughness  element  pattern  accounts 
for  25  percent  more  of  the  variance  when  added  to  the  equation  at  this  time. 
Note  that  the  amount  of  variance  explained  by  a  variable  is  dependent  on  the 
order  in  which  independent  variables  enter  the  regression. 

The  problem  of  estimating  the  relative  importance  of  the  independent 
variables  to  the  regression  equation  is  a  problem  that  cannot  be  directly 
assessed.  This  is  due  to  the  interrelationships  (correlations)  between  the 
independent  variables  (Table  VI-8  ).  If  all  cross  correlations  were  zero,  the 
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TABLE  VI-  9 
LOG-LINEAR  REGRESSIONS  FOR  a*/b 


Analys 

is 

Constant 

M 

F 

Ai 

A2     R 

stat 

R2 

1 

7.2*io"3 

0.193 

0.275 

1.41 

0.939  0.263 

82.9 

0.77 

2 

0.173 

0.179 

0.521 

-- 

-- 

68.0 

0.52 

3a 

0.100 

0.196 

0.288 

1.35 

0.943&  0.932b 

72.2 

0.87 

e 

0.057 

0.138 

0.221 

1.31 

0.977^  0.980^ 

65.0 

0.84 

a0nly  Roughness  Element  Patterns  B  and  D. 

Does  not  contribute  significantly  to  the  regression, 
e0nly  Roughness  Element  Patterns  B  and  C. 
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multiple  correlation  coefficient  (R2)  would  be  the  sum  of  the  correlations 
shown  in  Table  VI-10  under  R2.  In  regression  equation  1  this  sum  is  .867  as 
compared  with  the  actual  R2  for  the  regression  given  in  Table  VI-9  of  0.77. 
The  amount  of  variance  "explained"  by  the  addition  of  a  variable  in  this  problem 
is  directly  linked  to  the  order  of  the  addition. 

This  problem  can  be  partially  alleviated  by  the  calculation  of  Beta 
weights  (Table  VI-10).  Beta's  are  standardized  partial  correlation  coefficients 
that  can  be  used  in  a  regression  if  all  variables  are  in  standard  form. 
Standard  form  means  each  variable  is  subtracted  from  the  mean  of  that  variable 
and  divided  by  its  deviation. 

Standard  Form:  Z..  = 


x    Sx 


Beta  variables  can  be  viewed  as  a  standardized  form  of  b. .  This  removes  the 
problem  of  various  scales  and  units  between  independent  variables.  The 
disclaimer  "partial"  refers  to  the  fact  that  the  effects  of  variables  other 
than  the  one  to  which  the  weight  applies  are  held  constant.  With  these  facts 
in  mind,  it  can  be  said  qualitatively  that  opening  ratio  is  the  most  important 
variable,  abutment  variable  1  (spur  dike),  Froude  number,  and  the  roughness 
pattern  are  of  lesser  importance,  and  abutment  variable  2  is  of  least  impor- 
tance for  the  prediction  of  z*/b. 

Regression  equation  3  is  an  attempt  to  test  for  significant  differences 
between  roughness  patterns  B  and  D.  These  patterns  represent  uniform 
distributions  of  decreasing  density  (75  percent  to  50  percent).  More 
specifically  the  question  that  can  be  answered  is  "Does  the  knowledge  of 
which  element  roughness  pattern  (B  or  D)  explain  any  additional  variance  after 
the  other  three  independent  variables  have  been  included."  In  this  case  the 
additional  variance  is  an  increase  in  R2  of  .005.  For  any  practical  purpose 
this  increase  is  negligible. 
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TABLE  VI -10 

BETA  WEIGHTS 

bi 

Bi 

Ri2 

Froude  Number 

-0.56 

-0.393 

0.002 

Opening  Ratio 

-0.713 

-0.702 

0.477 

Abutment  Variable 

1 

0.150 

0.366 

0.287 

Abutment  Variable 

2 

-0.027 

-0.069 

0.079 

Roughness  Pattern 

-0.580 

-0.33 

0.022 
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Regression  equation  4  was  used  to  examine  the  effect  of  element 
roughness  patterns  B  and  C  on  the  prediction  of  a*/b.     These  have  flood  plain 
and  main  channel  differences  of  75-75  and  75-25  percent,  respectively.  The 
results  were  essentially  the  same  as  for  equation  3.  The  knowledge  of  the 
roughness  pattern  after  the  other  variables  are  taken  into  account  is  of 
marginal  practical  use. 

The  results  of  equation  3  and  4  can  be  accounted  for  in  two  ways: 

1.  The  effect  of  changes  in  pattern  densities  which 
changes  effective  Manning's  roughness,  is  accounted  for 
by  Froude  number,  or 

2.  The  parameter  i*/b   is  not  sufficiently  sensitive 
to  show  the  change. 

However,  due  to  the  correlation  (R)  of  Froude  number  to  roughness 
pattern  of  -0.553  (Table  VI-8)  the  former  is  considered  more  likely. 

A  summary  of  conclusions  from  the  regression  analysis  of  i*/b   on 
F,  M,  A,  Cf  follows: 

1.  Regression  equation  1  provides  an  accurate  representation 
of  the  i*/b   data. 

2.  The  use  of  Froude  number  and  opening  ratio  alone  explain 
a  little  over  50  percent  of  the  variance. 

3.  Froude  number  accounts  for  the  effect  of  changing  roughness 
densities  for  patterns  B  and  D. 

4.  The  effect  of  a  main  channel  of  lesser  density  after 
the  variables  F,  M,  and  A  are  included  is  marginal 
for  predicting  z*/b.     The  knowledge  of  the  main 
channel  differences  adds  little  to  a  regression 
relationship  beyond  that  added  by  the  other  variables. 
This  can  definitively  be  said  for  only  patterns  B  and  D 
and  has  not  been  extened  to  the  entire  data  base. 

5.  The  inclusion  of  spur  dikes  has  a  significant  effect  on 
£*/b. 
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The  regression  analysis  has  attempted  to  assess  the  significance 
of  varying  roughness  patterns  on  the  parameter  i*/b.     This  approach  is 
restricted  due  to  the  difficulty  of  assessing  absolute  importance  of 
independent  variables  and  the  lack  of  sensitivity  inherent  in  the  parameter 
2*/b.  Possible  future  analysis  of  different  variables  might  eliminate  the 
latter  problem  while  the  former  is  solved  only  by  further  experimentation. 
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ANALYSIS  OF  MAXIMUM  BACKWATER  HEIGHT  (h*/yQ) 

Maximum  backwater  height,  defined  as  the  maximum  difference  in 
backwater  surface  elevations  over  normal  elevation  (h*),  is  calculated  as 
described  in  the  Appendix.  The  normal  depth  y  is  the  average  of  the  three 
normal  depths  taken  over  the  north,  central  and  south  sides  of  the  flume. 
The  values  of  the  parameters  h*  and  y  used  in  this  part  of  the  analysis  can 
be  found  in  the  Appendix  under  headings  "HSTAR"  and  "YTILD,"  respectively. 

ANALYSIS  OF  VARIANCE   (h*/y   ) 

B o 

An  analysis  of  variance  was  performed  on  the  36  runs  comprising  the 
complete  set  for  roughness  pattern  C  (flood  plain  75  and  main  channel  25 
percents).  The  major  variables  and  levels  are  the  same  as  given  in  the 
previous  section  on  backwater  location.  The  dependent  variable  is,  of  course, 

h*/y„. 
J  o 

Table  VI-11  is  output  of  the  ANOVA  for  the  data  both  with 
and  without  spur  dike  results,  respectively.  Testing  (Table  VI-12) 
shows  that  Froude  number,  opening  ratio,  and  one  interaction  FxM  are  significant 
for  both  sets  of  data  at  the  95  percent  confidence  level.  The  significance  of 
the  abutment  variable  is  due  to  the  spur  dike  data. 

GRAPHICAL  ANALYSIS    (h*/y   ) 

Figures  VI-29  to  VI -35  are  plots  of  the  dependent  variable  h*/y  and 
opening  ratio  (M).  Numbers  included  in  the  plot  are  Froude  numbers  calculated 
over  the  entire  width  of  the  flume. 

A  review  of  the  plots  indicates  that  for  the  entire  data  base: 


153 


TABLE  VI-11 
ANOVA  OF  h*/y0  FOR  ROUGHNESS  PATTERN  C 
WITH  AND  WITHOUT  SPUR  DIKE  DATA 


WITH 

SPUR  DIKES 

WITHOUT  SPUR  DIKES 

Analy 
h* 

sis  of  Variance 
433 

Analys 
h* 

is  of  Variance 
432 

Level 

s  of  Factors 

Levels 

of  Factors 

F 
M 

A 

4 
3 
3 

F 

M 
A 

4 
3 
2 

Grand 

Mean    .12 

Grand 

Mean    .  1 1 

Source 
Variance 

Sums  of 
Sauarps 

Degrees  of 
Freedom 

Mean 
Squares 

Sums  of 
Squares 

Degrees  of 
Freedom 

Mean 
Squares 

Froude 

[F) 

0.11 

3 

0.03 

0.06 

3 

0.02 

Opening 
(M) 

Rati  o 

0.22 

2 

0.11 

0.14 

2 

0.07 

FxH 

0.04 

6 

0.01 

0.02 

6 

0.00 

Abutment  (A) 

0.00 

2 

0.00 

0.00 

1 

0.00 

F  x  A 

0.00 

6 

0.00 

0.00 

3 

0.00 

M  x  A 

0.00 

4 

0.00 

0.00 

2 

0.00 

F  x  M  x 

A 

0.00 

12 

0.00 

0.00 

6 

0.00 

TOTAL 


0.39 


35 


0.24 


23 
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XXX-FROUDE  NUMBERS  IN  FLOODPLAIN 
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1.  H*/y0  is  strongly  related  to  opening  ratio  and 
Froude  number. 

2.  Spur  dikes  appreciably  alter  the  relationship 
noted  above. 

3.  The  effect  of  various  roughness  element  patterns 
appears  slight  as  can  be  expected,  since  the 
value  of  y  varies  with  roughness. 

In  regard  to  any  further  analysis  this  would  indicate  that  differing 
wingwall  and  spill  through  abutment  shapes  do  not  appreciably  reflect  in  h*/y 
parameter  nor  do  the  various  roughness  element  patterns.  The  presence  or 
absence  of  spur  dikes  must  be  noted  due  to  its  strong  effect  as  shown  on 
Figures  VI-31  and  VI-34. 
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VII.  CONCLUSIONS 


In  relation  to  the  volume  of  data  collected,  the  information  thus  far 
obtained  can  be  viewed  only  as  a  partial  analysis.  Topics  bearing  on  the  two- 
dimensionality  of  the  data  have  not  been  given  the  attention  possible  in 
relation  to  the  data  collected.  Important  energy  transfer  mechanisms  due  to 
the  roughness  elements  have  not  been  investigated  nor  has  the  effect  of 
differing  roughness  patterns  on  the  basic  energy  relationships  been  adequately 
characterized.  These  are  topics  beyond  the  scope  of  this  contract  that  can  be 
pursued  in  later  analysis. 

Conclusions  reached  in  the  present  analysis  are: 

1.  The  primary  parameter  recorded,  water  surface  elevation, 
was  sensitive  to  every   level  of  the  independent  variables. 
The  independent  variables  were  Froude  number,  opening  ratio, 
and  abutment  shape. 

2.  For  any  use  of  the  data  other  than  an  analysis  of  water 
surface  elevations,  data  describing  wingwall  and  spill  - 
through  abutment  shapes  can  be  pooled  with  little  increase 
in  data  scatter. 

3.  The  presence  of  spur  dikes  constitutes  a  strong  influence 
on  the  location  and  height  of  maximum  backwater.  Any 
general  purpose  prediction  equation  for  these  and  other 
parameters  must  make  allowances  for  their  effect. 

4.  A  return  to  normal  depth  was  shown  to  occur  approximately 
5  to  10  ft  below  the  constriction  for  the  centered,  normal 
openings. 

5.  The  effects  of  the  constriction  on  depth  extended  upstream 
the  entire  length  of  the  flume  (60-foot  reach). 
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6.  The  wide  scatter  evidenced  in  the  backwater  coefficient 
data,  K*,  makes  impossible  any  direct  contribution  to 
modifying  the  calculation  procedure  in  regard  to  K*  for  the 
FHWA  backwater  method.  This  scatter  may  be  attributable  to 
the  small  h*  values  obtained  in  the  experiment. 

7.  The  parameter  l*/b   is  statistically  related  to  the  opening 
ratio,  Froude  number,  abutment  shape  and/or  spur  dike  and 
roughness  pattern.  Approximately  77  percent  of  the  total 
variance  is  accounted  for  by  these  variables. 

8.  A  graphical  analysis  of  the  parameter  h*/y  showed  a  similar 
relationship  as  demonstrated  above  for  the°dimensionless 
parameter  a*/b. 

9.  A  strong  relationship  exists  between  roughness  pattern  and 
Froude  number.  When  Froude  number  is  included  in  a  predictive 
equation  for  the  parameter  i*/b   the  relative  additional 
prediction  capability  provided  by  the  addition  of  the 
variable  for  roughness  pattern  is  marginal.  This  was  noted 

in  the  analysis  of  the  two  roughness  patterns  B  and  C  (Flood 
Plain  75  percent,  Main  Channel  75  percent  and  Flood  Plain 
75  percent  and  Main  Channel  25  percent),  respectively.  This 
may  be  partly  attributable  to  the  relative  insensitivity  of 
the  parameter  i*/b. 

10.  A  similar  comparison  between  roughness  patterns  B  and  D 
(Flood  Plain  75  percent,  Main  Channel  75  percent  and  Flood 
Plain  50  percent  and  Main  Channel  50  percent)  using  the 
parameter  i*/b   indicated  that  differences  noted  in  the 
parameter  solely  attributable  to  the  roughness  pattern 
were  negligible. 

11.  The  data  collected  in  this  experiment  do  not  provide  a 
sufficient  basis,  when  taken  alone  without  additional 
experimentation,  to  justify  major  changes  in  Hydraulics 
Design  Bulletin  No.  1  -  Hydraulics  of  Bridge  Waterways. 
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VIII.  RECOMMENDATIONS 


A  backwater  method  which  is  useable  to  the  highway  design  engineer  must 
have  two  characteristics: 

1.  Easily  understood  and  implemented  and 

2.  Reasonably  accurate  results. 

The  methods  currently  in  general  use  adequately  satisfy  the  first  criterion  but 
in  many  instances  not  the  second.  Future  research  must  be  directed  to  improving 
the  accuracy  of  the  various  methodologies  without  seriously  affecting  their 
level  of  difficulty  by  a  careful  blend  of  physical  model  and  prototype  informa- 
tion coupled  with  computer  simulation  using  advanced  numerical  techniques  such 
as  the  Finite  Element  Model. 

To  this  end,  the  recommendations  listed  below  deal  principally  with  two 
areas  of  interest:  improvement  of  the  data  collection  effort,  and  further 
analysis  and  interpretation  of  the  data,  particularly  where  research  efforts 
could  most  rapidly  advance  the  state  of  the  art. 


DATA  COLLECTION 

1.      The  Soundness  of  the  Large  Scale  Roughness  Elements  Must  be 
Scrupulously  Maintained. 

A  major  and  important  distinction  between  this  and  previous  experiments 
was  the  use  of  large  .scale  roughness.  The  simulation  of  prototype  conditions 
proved  sensitive  to  the  movement  of  fairly  small  portions  of  the  elements. 
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This  sensitivity  was  further  documented  by  including  varying  numbers  of  roughness 
elements  in  the  constriction.  In  order  to  maintain  reproducibility  of  results, 
a  program  of  constant  repair  and  attention  must  be  run  simultaneously  with 
experimentation. 

2.  The  Roughness  Elements  Must  Be  Spaced  to  Facilitate 
Varying  the  Length  of  the  Opening  and  Installing  Spur 
Dikes. 

Due  to  the  sensitivity  of  backwater  to  the  roughness  elements  the 
experimental  matrix  should  be  scoped  to  avoid,  to  the  extent  possible,  movement 
of  the  roughness  elements.  This  is  particularly  the  case  in  regard  to  the 
testing  of  spur  dikes. 

3.  Automated  Velocity  and  Depth  Measurement  Equipment 
Should  be  Used  to  Gather  Information. 

The  use  of  automated  sensing  equipment  requires  a  large  initial 
investment  which  can  often  appear  prohibitive.  However,  this  investment 
decreases  operating  costs,  shortens  data  collection  time  and  guarantees  a  known 
uniformity  of  measurement  error.  Automated  sensing  equipment  particularly 
lessens  human  measurement  and  recording  bias  and  is  highly  recommended  as  a 
worthwhile  additional  initial  expense. 

4.  The  Data  Base  Should  Be  Further  Expanded  to  Include 
Other  Slopes. 

The  present  configuration  of  the  flume  allowed  only  one  slope  to  be 
used  in  the  experiment.  The  simulation  of  other  slopes  would  allow  the 
expansion  of  the  data  base  to  cover  a  larger  set  of  prototype  conditions. 

5.  Undistorted  Models  Are  More  Suitable  for  Study  of 
Point  Effects  of  Flood  Plain  Roughness  and  Bridge 
Contraction  on  Backwater. 
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DATA  ANALYSIS  AND  INTERPRETATION 

1.  The  Data  Base  Should  Be  Used  to  Calibrate  the  Two -Dimensional 
Finite  Element  Model 3   Developed  as  Part  of  this  Study. 

In  many  instances,  useful  comparisons  between  particular  independent 
variables  were  not  possible  due  to  missing  runs  caused  by  time  and  money 
constraints.  Often  this  severely  limited  the  data  analysis  possible  and  the 
strength  of  the  conclusions  drawn.  This  is  particularly  evident  when  comparing 
the  effects  of  varying  roughness  element  patterns.  An  effort  directed  along 
such  a  plan  of  study  would  allow  "handbook"  information  to  be  developed  describ- 
ing turbulence  exchange  coefficients,  ranges  of  coefficients  for  varying  situa- 
tions, and  roughness  estimates  for  field  conditions. 

2.  The  Federal  Highway  Administration  Design  Bulletin  No.    1 
on  Hydraulics  of  Bridge  Waterways  Should  Be  Updated. 

The  information  developed  in  Recommendation  1  immediately  above  could 
be  used  to  update  Bulletin  No.  1  to  reflect  the  additional  findings.  The 
emphasis  in  this  task  is  on  the  needs  of  the  design  engineer,  which  are 
distinctly  different  from  those  of  field  data  collectors  or  theoretical 
hydraulicians . 

3.  Other  Data  Analyses  Beyond  the  Scope  of  this  Effort  Should 
Be  Performed. 

Other  recommended  data  analyses  include: 

a.  Investigate  the  shape  of  the  maximum  backwater  as  a 
function  of  the  various  main  channel  and  flood  plain 
roughness  element  densities.  Since  a  large  portion 
of  the  data  collected  includes  difficult-to-obtain 
flood  plain  information,  this  provides  a  logical 
extension  of  the  current  study. 

b.  Compare  selected  portions  of  the  collected  data  with 
field  behavior. 
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Evaluate  the  revised  Geological  Survey  method  for 
computing  average  approach  flow  line  length  and 
exit  losses  using  the  flume  data. 
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